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7 THE DIFFRACTION OF X-RAYS AND ELECTRONS 


BY FREE MOLECULES 


By M. H. PIRENNE, Dr.Sc. (Liége) 


2 plates, 82 text-figures. 


I2s. 6d. net 


This new volume in the Cambridge Series of Physical Chemistry is intended to give an account 
of the theoretical basis of the study of x-ray diffraction by gases, and of the information it 


has yielded about the structure of atoms and molecules. 
attention has been paid to the hypotheses and 


ideas and experimental results and special 


The book deals with fundamental 


principles underlying theories as well as to their limits of validity. 
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UNIVERSITY 


PRESS 
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Binding Cases for volume 57 (1945) and 
previous volumes may be obtained for 
6s., inclusive of postage. For 9s. the six 
parts of a volume will be bound in the 
publisher’s binding cases and returned. 


THE PHYSICAL SOCIETY 


1 Lowther Gardens, Prince Consort Road, 
London S.W.7 


——_ 


SMALL TRANSFORMERS AND CHOKES 
to meet any practical specified requirements for 
Research and Experimental Purposes 


We have considerable experience in the design and manufac- 
ture of special transformers and chokes for a wide variety of 
purposes and include amongst our clients many well-known 
Research Laboratories, Government Institutions, ; Technical 
Colleges and Industrial Research Departments. 

Dn receipt of data giving the desired electrical performance 
ind circuit conditions we will quote for prototypes or 
juantities without obligation and are in a position to give 
yery quick delivery. Our prices are reasonable, our goods of 
uigh quality and our technical service unequalled. We also 
tock most of the conventional types of transformers and 
chokes. 


SOWTER TRANSFORMERS, 
27 Homefield Road, Wimbledon, S.W.19. 
*Phone : WIMBLEDON 0244, 


THE SIR JOHN CASS TECHNICAL 
INSTITUTE, 
Jewry Street, Aldgate, E.C.3 


DEPARTMENT OF PHYSICS 
AND MATHEMATICS 


Full-time and part-time courses for B.Sc. 
(London) degrees and for post-graduate 
studies 


Special courses in— 
X-RAY SPECTROSCOPY 
ILLUMINATION 
MATHEMATICAL PHYSICS 
ACTUARIAL MATHEMATICS 


Session begins September 23rd, 1946 


Full particulars may be obtained on 
application to the Principal 


OptTIcAL BENCH 
PHOTOGRAPHIC ENLARGER 


PHOTOGRAPHIC COLOUR 
FILTERS 


OPTICAL FLATS 


THE 
LEECH OPTICAL Co. Ltd. 
WHETSTONE LEICS. 


Telephone: NARBOROUGH 2213 


PROJECTION EQUIPMENT 
MICROSCOPY 


APLANATIC MAGNIFIERS 
PRISMS TO ANY SPECIFICATION 


y, 
LEECH 


BLOOMED SURFACES 
ALUMINISED SURFACES 


RHODIUMISED SURFACES 
MIRRORS OF ALL TYPES 


London Office : 
94, HATTON GARDEN, E.C.1 
Telephone: HOLBORN 1752 
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SELENIUM PHOTOCELLS 


for use with direct reading instruments or valve amplifiers. | 


TEST TUBE COLORIMETERS | 
for Haemoglobin and other determinations. l 
| 


REFLECTOMETERS 


for comparing reflectance from all surfaces. 1 


i} 


EVANS ELECTROSELENIUY 
LIMITED | 
BISHOP’S STORTFORD. HERTS. 


ELECTRICAL AND INSTRUMENT SP 


Greater tensile, elastic and fatigue strengths 
than any other non-ferrous alloy, a higher conductivity than | BERYLLIUM COPPER 
any of the bronzes and excellent resistance to corrosion and 
wear—these characteristics of Mallory 73 Beryllium Copper 
have made it first choice for instrument springs, diaphragms 
and bellows, current-carrying springs, snap action switch 


Properties of 


blades, contact blades and clips. MALLORY 73 BERYLLIUM COPPER | 

Supplied annealed or lightly cold worked, it has good | aicersheat cesineie 

forming properties and is readily fabricated into springs : . 

and parts of complicated shape. A simple heat treatment ee net 

then develops its remarkable properties. Ae cy . mn 
Limit of proportionality 


Available as sheet, strip and wire, in a range of tempers to tons per sq. inch 47.50 | 
suit users’ requirements, and as rod, tube, precision rolled 


: : : Fatigue limit 
hair-spring strip, and silver-faced contact bi-metal strip. g 


tons per sq. inch + ]9-29 
pf: Vickers pyramid hardness | 

Electrical Conductivit 
per cent 1.A.C.S. 23.25 


J OH NSON, MATTHEY & CoO. LIMITED, 
Controlling MALLORY METALLURGICAL PRODUCTS LTD. | 
73-83 Hatton Garden, London, E.C.i 


Full details are given in our beok!st, | 
Telephone: HOLborn 9277 


which will be sent on request. 
G.D.40 
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DORAN. 


VITREOSIL 
POTENTIOMETERS APPARATUS 


BRIDGES A.C. & D.C. 
GALVAN OMETERS 
STANDARD CELLS 


ELECTRONIC TEMPERATURE 
CONTROLLERS 


are a few of the wide range of Instruments 
of modern design and high accuracy now 
in course of manufacture for all types of 
laboratory measurement and_ industrial 
process control This photograph—by the way, it won 
International Exhibition in 1939—shows 
by some pre-war apparatus made to speci- 
fication in VITREOSIL pure fused silica, © 
but you need hardly be reminded that 
this work can always be done for you. 


Sere els, THE THERMAL SYNDICATE LTD. 


Tel. : Stroud 15. Head Office: Wallsend, Northumberland. 
London Depot: 12-14, Old Pye Street, 
Westminster, S.W.1. 


Doran Instrument Co. Ltd. 


ACN 


InST 


FOR INDICATING - RECORDING: CONTROLLING é 


i 
TEMPERATURE ee 


We shall be glad to supply full particulars 


BOWEN INSTRUMENT Co. Ltd. 
9/13, NEWTON ROAD ~- LEEDS, 7 


Telephones: 41036/7 


‘ 


Write for characteristics 


BASICALLY BETTER - 
AIR -S PACED 


TRANSRADIO LTD. t6 THE HicHwa 


CO=AX /W/0S CABLES 


iv 


VISUAL TWO-DIMENSIONAL 
delineation of any recurrent law. 


RELATIVE TIMING OF EVENTS 
and other comparative measure- 
ments with extreme accuracy. 


PHOTOGRAPHIC RECORDING 
of transient phenomena. 


SIMULTANEOUS INDICATION 
of two variables on a common 
time axis. 


Completely embracing all the above func- 
tions, the unique Cossor DOUBLE BEAM 
Cathode Ray Tube, as embodied in Model 
339 Oscillograph, is inherently applicable 
to all problems arising in the 


RECORDING, 
INDICATING and 
MONITORING 


of effects and events in Electrical engineer- 
ing ; and in Mechanical engineering when 
the effects can be made available as a vol- 
tage. Recurrent traces are studied visually 
and transients may be photographically 
recorded with Model 427 Camera. 


A. C. COSSOR LTD., 


Instruments Dept. 
Cossor House, London, N.5 


*Phone : CANonbury 1234 (30 lines). 
’Grams : Amplifiers Phone London. 


SERRA RG ORIEN ERT 


CER] 
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CARPENTER HIGH SPEE} 

POLARIZED RELAY 

Main features of Standard model : Hig} 

speed. Short transit time. Negligible di] 

tortion up to 300 c/s. Contact gap | 
function of input power, hence small d 

No contact chatte} 


tortion almost down to failure point. His 
contact pressures. 


High sensitivity—robust operation 
5 mv.a. at 100 c/s. or 0-2 mw. D 
Great ease of adjustment. Magnetic bil 


adjustment giving absolutely smooth co} 
trol. Balanced armature—hence immunij 
to considerable vibration and no positio 
error. 


Full details and Test data available on request. 


TELEPHON 
MANUFACTURING 
co. LTD. 


Contractors to British - Dominion = Colonial and 
Foreign Governments. ' 


HOLLINGSWORTH WORKS, | 
DULWICH, LONDON, S.E.2! | 


Telephone 
GIPsy Hill 2211 (10 lines). 
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ANSON 


\ 


CBA: 


STASBILISTOR 


he ONLY constant voltage equipment giving a 
onstant UNDISTORTED output with changes in 


ains voltage and load. Available as a complete unit 


) cover, as a wired-up Chassis as illustrated, or as 
set of components. 


Undistorted 
output wave 


Constant out- 
put with mains 
changes of 15% 


Constant out- 
put from zero 
to full load 


Rapid response 
to mains or load 
changes 


Special voltage and frequency models now available 
which, at the designed load, will hold the voltage 
between the extreme limits of + 0°14% of rated out- 


put voltage for simultaneous variations of + 6% in 
voltage and + 2% in frequency. 


Write for Publication E.E.2 to Dept. P.S 


ESTINGHOUSE BRAKE AND SIGNAL CO. LTD. 


YORK WAY, KING’S CROSS, LONDON, N.I. 
SS EE LT EET DE SE TE OCTET SIE IS 


DRAYTON 
F lay” 
Metal 
Bellows | 


i 


WG 


SS 
sss 


eS 


WEES 


These bellows are formed from the initial 
tube in one gradual continuous operation, 
resulting in a uniformity of wall-thick- 
ness unattainable by any other method. 
The tough resilient product is tested to 


many times the maximum rated working 


pressure during formation. Customers’ 


end plates can be fitted prior to forming 


so that the soldered joint is also pressure 
tested. 


For Gland Seals ; Refrigeration 
Control; Thermostatic and 


Pressure Operated Devices, etc. 


Every Gland uniform in life and performance. 
Pretested during forming. 


Absolutely reliable in operation. 


@ 
@ 
@ 
@ Provision of closed end eliminates one joint. 


Root Diam. 3/8” to3” Outside Diam. 9/16” to 44” 


Write for the ‘‘ Hydroflex ’’ Brochure. 


DRAYTON REGULATOR & INSTRUMENT CO. LTD, 
WEST DRAYTON West Drayton 2611. MIDDX. 
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THE CMG8 


The CMG 8 is a typical example of the OSRAM range of emission- 
type photocells which constitute the essential means of converting 
light changes into electric current. Widely used in sound 
projectors and industrial apparatus they are non-microphonic 
and of convenient size. Outstanding features include :— 


# Linear response for sound reproduction giving undistorted output. 


High sensitivity to artificial light. 


A detailed technical data sheet is available on request. 


Osram &6.C Osram 


PHOTO CELLS CATHODE RAY TUBES VALVES 


a i 7 


way, London, W.C.2. 
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Telecommunication Measuring 
instruments. 


Telecommunication Components, 
Loading Coils for Telephone 
Cables. 


Toroid Cores, Coils and 
Transformers. 


Magnetic Dust Cores for all 
frequencies. 


Magnetic Powders. 


Quartz Crystal Units for 
frequencies up to 20 Mc/s. 


Rectifiers for instruments. 


Selenium Photocells etc. etc. 
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RESEARCH & 
ENGINEERING 


BTH Products are available 
as instruments and apparatus 
for research, development, and 
engineering, including all branches 
of electronic control. 


THYRATRONS MAGNETRONS 
CRYSTAL RECTIFIERS 
SPECIAL VACUUM ELECTRONIC DEVICES 
GLASS, AND GLASS METAL SEALS 
RECTIFIERS 
FRACTIONAL HORSEPOWER MOTORS 
‘ETC. 


RUGBY 


THE BRITISH THOMSON-HOUSTON COMPANY LIMITED, RUGBY, ENGLAND. 


i 
| 


A3595N 
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HILIPS GEIGER-COUNTER X-RAY SPECTROMETER 


. . The new Philips Geiger-Counter X-Ray 
Spectrometer is a revolutionary development 
based on newly applied principles of diffraction 
technique. Using a direct-reading method, no 
film is required. The speed of analysis is increased 
and greater accuracy obtained. The Spectrometer 
is simple, rapid and easy to operate. It provides a 
fast medium for the qualitative and quantitative 
analysis of crystalline and certain amorphous 
materials, quickly identifying chemical substances 
present and their state of chemical combination. 
Under optimum conditions of resolution, diff- 
raction angles are obtained with an accuracy 
of 0:03 degree. The reflected beams are located 
by means of a highly sensitive Geiger-Counter 
tube mounted on a movable arm that traverses 
the scale. Radiation intensity is integrated and 
totalized by an electro-mechanical counter, while 
average intensity is indicated by a micro-ammeter, 
both connected to the Geiger-Counter through 
electronic circuits. 


HILIPS X-RAY DIFFRACTION APPARATUS 


PHILIPS 


LIPS LAMPS LTD, CENTURY HOUSE, 


. . . This apparatus constitutes the ideal flexible 
diffraction machine. With a power unit rated at 
60 kvp. 50 ma. any of the modern techniques can 
be readily employed. The full wave rectification 
means a longer tube life because the tube filament 
runs cooler, and for equal target loading the 
filament emission is only about half of that 
compared with self-rectified or half-wave recti- 
fication. Furthermore, contamination of the 
tube target, due to deposits of tungsten, is thereby 
greatly reduced. This results in a clean tube 
target, no contaminated radiation to complicate 
the pattern on the film, and longer, more efficient 
tube life. Among the many applications for this 
machine are :—Chemical Indentification . Phase 
Analysis . Solid Solubility Study . Single Crystal 
Studies . Particle Size Measurements . Metallo- 
graphic Studies . Molecular Weight Determination 
Structure Identification . Stress Analysis . Fibre 
Studies . Process Control . Orientation Studies. 


METALIX 


SHAFTESBURY AVENUE, LONDON, W.C.2, 
a 
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RAY TUBES 


ILFORD RECORDING FILMS AND PAPE| 


Ilford Limited market a comprehensive range of sensitised materials for us 
recording instruments of all types. The following brief details cover a seleg 
of films and papers specially prepared for Cathode Ray Tube record 


FILMS 
5R101 Panchromatic, very high speed, medium contrast. Specially suitable for recording 
8R101 = fluorescent (phosphate) screens. Sometimes used for green and blue screens. 


5G91 Orthochromatic, very high speed, high contrast. The fastest material for recording | 
and blue fluorescent screens. Suitable for high-speed transients. 


5B52 Blue sensitive, high speed, medium contrast. Specially for recording blue fluor@ 
(orange screens. Minimum image spread enables successful recording in one trace of 
base) variations in writing speed. Suitable for high-speed processing. 


Films 5/1000 in. thick (Code No. beginning with 5) supplied in rolls. 
Films 8/1000 in. thick (Code No. beginning with 8) supplied only in flat sheets. 


PAPERS | 
BP1 Blue sensitive, high speed, high contrast. Glossy surface. Medium thickness | 


Widely used for recording blue fluorescent screens and for many other forms of reco} 

where a paper base is required. | 
BK1 Blue sensitive, high speed, high contrast. Similar to BP1, but on thin base with 4 

matt surface. 
Rolls of film and paper are supplied in 35, 60, 70, 120 and 300 mm. widths and in 25, 50 and i 
lengths. 35, 60 and 70 mm. widths can be supplied with standard 35 mm. positive perforatio| 
order. Film is also available in 16 mm. width (perforated) — 50 ft. and 100 ft. rolls only, 


limited range of special sizes and packings can be supplied subject to certain conditions which 
be notified on request. | 


It is highly desirable in the interests of uniformity that designers of new apparatus 
should work within the limitations of British Standard Specification No. 1193/1945 
entitled “ Light- sensitive Film and Paper for Recording Instruments.” | 


ILFORD LIMITED - ILFORD - LONDC 


Zea 
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TYPE 
AF 200 
FREQUENCY RANGE HARMONIC CONTENT 
50 to 20,000 c/s. Less than 3% over main portion of the 
ACCURACY frequency range at full rated output. 
+ 1% + 2 qs. R.F, CONTENT 
OUTPUT Less than 0:25% at full output. 
2 watts into 600 or 10 ohms. DRIFT 
HUM Less than 10 c/s a day after first half- 
Less than 0°5% at full output. hour of operation. 


PRICE £50 


ROR EHAM WOOD-HERTS 


TELEPHONE : ELSTREE 11357 | 
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Tue world-wide use of “AVO” Electrical 
Testing Instruments is striking testimony to their 
outstanding versatility, precision and reliability. 
In every sphere of electrical test work they 
are maintaining the ‘‘AVO’’ reputation for 
dependable accuracy. which is often used as a 


standard by which other instruments are 


XIV THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


LLARD-B.Ti 


ELECTROMETRIC 
TITRATION APPARATUS | 
TYPE-E920 


Ss ey 7 700 
| 
| | 
| | i} 
Reais han ane 650 
WITH W/20 FERROUD AMMONIUM SULPHATE | 
ELECTRODES: BRIGHT PLAT(NUM AND | 
aa SATURATED CALOUEL | 


ii mV 
500}-- — | 500 
450 
boi) ry 2 3 * al 400 
«Co Fea 
CCS Fe $O4 


Titration of Vanadium in Steel with N/20 Ferrous 
Ammonium Sulphate Electrodes: Bright Platinum and 
Saturated Calomel. 


Tuis apparatus has been designed to supply the need 
for a robust, self-contained potentiometric titration 
apparatus which is sufficiently simple to be used in routine testing by unskilled 
operators and is yet capable of meeting the requirements of the industrial research 
chemist. 

The apparatus is operated from the 50-cycle supply mains and the end point is 
detected by a “‘ Magic Eye” indicator. Thus there are no batteries to be replaced 
and no delicate galvanometer to be damaged by mechanical shocks or electrical 
overloads. A special circuit eliminates all possibility of drift during a titration, 
and changes of mains supply voltage do not give rise to any inaccuracies. 


OTHER TYPES OF MULLARD APPARATUS 


HIGH AND LOW FREQUENCY A.C. BRIDGE EQUIPMENT. GRAIN MOISTURE METE 


THES ROGCEE DINGS OF DHE PEYSIGAL SOCIE DY. XV 


From magnetic pole to magnetic pole, wherever magnetic circuits 
have to be designed for high efficiency, more and more 
manufacturers are specifying %* “ Ticonal” permanent magnets. 


* Mullard “ Ticonal”’ magnets are anisotropic, i.e., they possess greatly 


increased properties along the preferred axis. They were the first 
anisotropic magnets commercially available : they are still the best. 


‘TICONA 


ERMANENT MAGNETS MADE BY 


ULLARD 


) MULLARD WIRELESS SERVICE CO. LTD., CENTURY HOUSE, SHAFTESBURY AVE., LONDON, W.C.2z. (8D) 


" 


(REGD. TRADE MARK) 
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Apparatus 


An Improved Lightweight 
Instrument for the Psychiatrist 


The advantages of this new apparatus are 
that it is simple and quick to operate, it is 
absolutely safe both for patient and doctor, 
and it is light in weight and handy in size. 
It is robust so that it stands up to the strain 
of routine practical use, and it has been 
designed in such a way that the risks of 
technical breakdown have been reduced to 
a negligible minimum. 


Its improved design is based on the rejection |]} 
of the idea that the convulsant dose can be 
calculated in advance by measuring the 
resistance of the patient. Circuits for the 
measurement of this resistance have there- 
fore been eliminated, thus — simplifying 
not only the construction but also the 
working of the apparatus. (Gerald Caplan, 
J.Ment.Sci., 1945, 91,200). 


9 A new type of electronic time switch is 

incorporated, which gives an entirely safe 
andreliable performance under aliconditions. 
An arrangement of push button switches 


prevents accidents to the patient or the 
nursing staff while the electrodes are being 


STANDARD MODEL. 200/250 volts, 50 
cycles A.C. The output is 40-150 volts for 
0.1-1.0 second. It is housed in a strong metal 


cabinet measuring 9 in. X 9 in. x 8 in., and adjusted, and at the same time reduces the 
weighs 16 Ibs. — ; actual procedure of administering the shock to 
PRC Ft SBE CS en es a its simplest terms. For this reason, and also 
Models are avaiiable for immediate delivery on x 

14 days trial. because the apparatus is so small and easy to 


handle, the treatment timeis greatly reduced. 


MULTITONE ELECTRIC GCOMPANY, LIMITED. 


223+7, St. John Street, Clerkenme! [eel omic nee ore 


C.R.C. 103 
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POTENTIOMETER 
TYPE D-72-A 


For all measurements based on potential difference 


XVit 


sasures from — 0.01 volt to +-1.92 volts on three ranges 
Pal, <O.L and <0.0r. 


allest subdivision on X 0.01 range, IO microvolts. 
all negative reading ensures absolute zero volts 


sition. ACCESSORIES FOR USE 
curate current control by stud rheostat and slide wire. WITH THIS INSTRUMENT 
arate standard cell balancing circuit. VOLT RATIO BOX .._ .. TYPE A-202-A 
rminals for two external circuits selected by switch RTE OD ens oan any Sones ce 
: GALVANOMETER _... .. TYPE D-41-A 


LAMPSTAND .. oe .. TYPE D-74-A 


Muirhead & Company Ltd., Elmers End, Beckenham, Kent. Tel. : Beckenham 0041-2. 
OVER 60 YEARS DESIGNERS & MAKERS OF PRECISION INSTRUMENTS 


C.B.G. 35 
) 
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THE “SHIRLEY” 
MOISTURE METER 


Originally developed for the rapid determination 
of the moisture regain of cotton, the instrument is now 
available for a wide range of materials, including wool, 
rayon, flax, jute, hemp, etc. 


XVili 


Information concerning Meters suitable for 


testing other materials will be given on request. 


oe ee aes 
The instrument is both quick and accurate, giving 
instantaneous direct readings of Moisture Regain. 
Worked from A.C. Supply (consumption is only 25 
watts), the meter is so simple that an unskilled person 


can obtain correct readings. 


THE RECORD ELECTRICAL CO. LTD. 
BROADHEATH : ALTRINCHAM : CHESHIRE 
Phone : Altrincham 3221/2 ’Grams : “ Infusion’’ Altrincham 
London Office : 28 Victoria Street, Westminster, S.W.!. 


LARGEAR 
CRVARIOGE 
EN > 


Electrical Specn. 2 € 


Variable amplitude, out- 
put 2 watts max. Self- i 
contained power supply Sy 
for 200-250 v. 50 c.p.s. SS 
input. Output imped- 
ance 500 & or 5,000 Q as 
selected by switch. Pre- 
set contrcl for matching 
replacement oscillator 
valves. 


SPA beregr OY 


A very useful Sine Wave Oscillator supplying a 

pure waveform of fixed frequency, equal to 

expensive B.F.O.’s. Suitable as excitation 

source for bridges or as laboratory standards. 

Standard oscillators 500, 900 or 1,000 c.p.s. are 

available. Oscillators of other frequency can 
be supplied to order. 


LABGEAR 


DESIGNERS AND MANU-ACTURERS OF 
HIGH CLASs 
INSTRUMENTS & ELECTRONIC APPARATUS 
WILLOW PLACE, CAMBRIDGE 


Telephone : 
Cambridge 245 4 (2 lines) 


Telegrams : 
Labgear, Cambridge 


Phone: Abbey 5148. 


BGEAR AUDIO OSCILLATOR 


easy valve replaceme 


a \ 


| 


Mechanical Specn 


Solidly — construct 
metal cabinet of ple 
ing modern style. Sil} 
li x 8x9 ins. high. 

movable back cover f: 


Finish: Black — wrinl 
paint with nickel plat 
fittings. 


Weight: 20 los. appr 


Full specification of th 

and other Labgear La’ 
oratory Instruments a! 
available on _ reques}j 
Please mention th 
journal when writin\}) 


| 
f 

| 

i 

| 
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Progress hy Guality 


AVIMO 


35 mm. CAMERA 
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ABSTRACT. The estimation of the best combination of lens and film for aerial 
reconnaissance photography is not a simple problem, owing to the number of variables, 
such as focal length and aperture of the lens, the size and speed of the film, forward motion 
and speed of the aircraft which are involved. For quantitative considerations it is necessary 
to adopt a measure of the detail-revealing capacity of the negatives, and for this purpose 
the resolving power was chosen, since it was not only obviously related to the capacity 
of the negative to reveal detail, but also capable of justification as a possible means of 
estimating the performance of the lens in optical image formation. In order that the 
resolving power figures may correspond as closely as possible to practical performance, 
careful choice of test-object was necessary, for experiment showed that different test- 
objects gave different results. A test-object consisting of two rectangles separated by their 
width, with a difference in density of 0:2 between them and the background, was adopted. 
Apparatus with such test-objects at the focus of a collimator of long focal length was set up, 
and with it measurements were made of the photographic resolving power all over the field, 
and in a set of planes embracing the focus, for a representative collection of air camera 
lenses. The results showed that the photographic resolving power improved as the lenses 
were stopped down from their maximum apertures until very small apertures were reached, 
and that as a general rule the resolving power diminished towards the outer parts of the field. 
The variation with aperture was accounted for, in general terms, by theoretical considera- 
tions, and certain regularities in the change of resolving power with angular distance from 
the optical axis were also explained. A figure for the performance of any lens was obtained 
by averaging the resolving power over the field of the lens, having regard to the greater 
area covered at the larger angles. It was found that the performance of any of the lenses 
could be represented with very fair accuracy by the formula 


ade ae : E No. |%8 
LG tan? 0 


where R is the mean resolving power in lines per mm. over the field up to an angular 
distance @ from the axis, f is the focal length of the lens in inches, and F’ No. its relative 
aperture. G is the granularity of the film in density-microns at a density of 1-0. 

This formula has no substantial theoretical basis, although it can be accounted for 
to some extent. It has proved useful, for it assists in the choice of the focal length, aperture 
and angle of view most appropriate for any given purpose in aerial reconnaissance photo- 
graphy. It may also be used, at least until the progress of lens design renders it inaccurate 
in form, for estimating the merit of a given lens design by comparison of its performance 


with that given by the formula. 


> 


§1. INTRODUCTION 
N order to make clear the reason for examining the problem with which this 
paper deals, a brief sketch is given in this Introduction of the work on aerial 
photography of which it forms part. 
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Soon after the war began it seemed likely that photography from aircraft would } 
be an important item in it. In spite of the considerable development of aerial) 
photography in the war of 1914-1918, the problem of estimating what were the}, 
best possible results, and how they could be secured, still seemed to be largel | 
unsolved, especially on the theoretical side. Consequently, very early in the war, | 
the Kodak Research Laboratories took much interest in this problem in order tal} 
discover what sort of film was best for aerial photography. Preliminary discussions}) 
soon showed that it was one of some complexity. The main use of aerial photo- | 
graphy was likely to be in reconnaissance, in which it is important to record they 
utmost possible detail. The other use, for map making, requires distortion of | 
the image introduced both by the lens and by the deformation of the film base tc}f 


I) 


be reduced to the minimum, and presumably some loss of definition would bat | 
tolerated if that were the necessary price to pay for absence of distortion. Aerial) 
survey for map making, however, did not become an urgent problem till later in 
the war, and this paper is concerned with one broad aspect of the general problem 
of securing the best possible definition in aerial photographs. | 

The nature of the general problems may be gathered from the following| 
considerations. It is necessary to give at least a certain minimum exposure if thef 
negative is to be of value. Thus the time of exposure T, the focal ratio F. No. 
the brightness of the subject B, and the speed of the film S must be in a ve 
definite relation. This relation, in fact, is 


Soc(F. No.)?/BT. 


To see that this is so, it is only necessary to note that it constitutes the practical] 
definition of S. For a given value of B, therefore, a wide variation in values of S'jf] 
F.No. and T is possible, which will, nevertheless, satisfy the condition that thef 
exposure must be adequate. But out of the possible sets of values of S, F.No. an 
T some will lead to sharper photographs than others, as may be seen in the following 
way. Itis amatter of common experience that sharper photographs are obtained] 
with slow films than with fast ones. But if the speed of the film is decreased therjf 
either J must be increased or F'. No. decreased, or both. An increase in 7} 
however, necessarily causes an increase in the blurring of the photograph by the} 
forward motion and vibration of the aircraft. Also it is a matter of common! 
experience that a moderate reduction in the aperture of photographic objectives o} : 
normal construction improves the sharpness of the photographs taken with themi] 
so that a decrease in F'. No. may be expected to diminish the sharpness. Thal 
increase in sharpness expected from the use of slow film is therefore offset, more : 
or less, by the increase in time of exposure or aperture, or both. Hence, unless sf 
remarkable exactness of balance occurs between the various effects involved, web 
must expect that some particular set of values of S, F.No. and T will lead to the 
sharpest photographs. | 

The foregoing argument merely demonstrates that there is probably ai 
optimum as regards sharpness. In order to understand how this optimum ij} 


definition and how the optical properties of the lens are affected by a change is} 
aperture. 

A little consideration will also show that we must know more about the lens tha ' 
the variation of its properties with aperture. On general grounds it is clear that 


Hf 
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the larger the image on the film can be made, even if the faults in the image given 

by the lens are magnified in proportion, the less serious will be the relative degra- 
_ dation of the image due to the film. We are thus led to expect that better results 
_ might be obtained with lenses of long focal length than with short. Fora given size 
of film, however, the area of ground covered in a single photograph is less the 
| longer the focal length of the lens. From an operational point of view it might be 
(necessary to cover a large area of ground not only for estimating, for example, 
‘ the damage suffered by a town from aerial bombardment, but also because it is not 
easy to get a given target into the centre of the picture under operational conditions. 
. Thus it is also important to know whether the angle of view for which a lens is 
|designed has any important effect on the definition. 
_ Two considerable gaps in our knowledge thus existed. One was the way in 
which the film influences the definition, and the other the way in which the 
definition obtained with a lens of given design depends upon its aperture, focal 
length and the size of the film with which it is used which determines its operative 
angle of view. Until these gaps were filled it could not be asserted with certainty 
, that the best was being achieved with available equipment, nor was it possible to 
indicate the most profitable directions in which improvement could be made. 
It might be suggested that this difficulty could be overcome, in empirical fashion, 
by making aerial photographs with every possible combination of lens and film. 
It is, however, neither easy nor inexpensive to make reliable tests in an aircraft, and 
the amount of experiment to arrive at an incontestable conclusion, when so many 
variables are involved, would be prohibitive. 

To begin with, a purely theoretical attack was made. It was hoped that it 
might be possible so to relate the resolving power of a photographic lens, as 
.observed visually, for preference, and the resolving power of the film, as found with 
a lens of extremely high resolving power, or with a test-object in contact with the 
film, that the resolving power of the combination of photographic lens and 
film could be estimated. The attempt proved abortive. Experiments carried 
Jout in the Kodak Research Laboratories, and a further more extensive series 
carried out by Mr. R. J. Hercock of the Ilford Laboratories, revealed no relation 
f the sort expected. Mr. Hercock, in fact, did discover an empirical relation 
which fitted his results, but the value of this for the desired use seemed doubtful. 
However, further theoretical examination of the problem, on similar lines to the 
riginal attempt, but this time bringing in the grain of the film explicitly as a 
(factor disturbing the recognition of detail, suggested a relationship between the 
resolving power obtained in the negative and the quality of the image formed by the 
jlens, the contrast, turbidity and granularity of the film, and the contrast of the 
object photographed, which was not in serious disagreement with any of the known 
facts. Previous to this work no satisfactory reason had been found, for instance, 
for the fact that photographic materials display a maximum in the curve of resolving 
ower against exposure or density. ‘The turbidity of the emulsion layer had been 
recognized as a feature reducing the contrast of the image formed in the layer, and 
it had been recognized that on development this contrast would be increased or 
decreased according to the slope of the (density: log. exposure) curve of the 
material, at the point of exposure. ‘Thus the resolving power would be expected 
to increase with density and then to remain constant throughout the region where 
the slope of the (density : log. exposure) curve remains constant (the straight-line 
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region of the characteristic curve). In fact, the resolving-power curve exhibits 4) 
maximum near the commencement of the linear region of the characteristic curve) 
This was accounted for by the continual increase of granularity with density | 
because an increase in granularity reduces the visibility of any given density) 
difference in the developed image. Advantage was taken of the recognition of thi 
property to construct an instrument for measuring granularity, which proved olf] 
great value later on in the work on aerial photography. An account of this 
instrument has already been published (Romer and Selwyn, 1943). | 

By means of this theory it was possible to account, qualitatively at least, for the} 
very great differences frequently observed between the resolving powers obtained] 
visually and photographically, for the fact that the resolving power in prints was|| 
not greatly affected by the contrast of the printing paper, and for the considerable} 
increase of graininess of prints made from negatives taken under hazy atmospheriq}} 
conditions compared with those taken under clear conditions. It was realised i 
however, that more exact establishment of this theory would require much timeif] 
For the time being, therefore, the theory was used only as a guiding principle,}) 
and attention was concentrated on objectives of greater practical importance. 

The initial attempts to estimate the effect of varying the aperture, angle of views) 
and focal length were speculative and gave only a fragmentary idea of the nature off 
the variation. Although, therefore, they suggested the general direction in whic 
the definition might be expected to vary, as the aperture, angle of view and focall, 
length were changed, it was clear that many tests on lenses would be necessary} 
before it could be said that the performance of camera objectives in aerial phot 
graphy was understood. Some of the early arguments are used, in a for 
modified by later experience, in the body of this paper. 

During the time these attémpts at the problem had been in progress, numerou 
discussions had taken place with the staff of the Royal Aircraft Establishment andi 
some degree of agreement had been reached on the type of work which could be} 
carried out in the laboratory, in contrast with that which must necessarily be done 
in aircraft. At this stage, the Kodak Research Laboratories were placed unde 
contract to the Ministry of Aircraft Production, with a fairly definite programme o 
investigations in view. Foremost among these investigations was the establisa| 
ment of a technique of photographic testing of air camera lenses by which theij 
performance in aerial photography could be estimated. It is with this problenifj 
and the results of systematic tests of all the lenses available that this paper if 
concerned. ‘The results obtained enabled a test to be made of the belief, fo: 
which only slight and doubtful evidence was originally available, that thef 
differences between lenses could be accounted for, for the most part, by differenced 
in angle of view, aperture and focal length. In other words, given the foca 
length, aperture and angle of view, it should be possible to estimate with moderatd, 
accuracy the definition given by the lens. It will be seen that this was so for thd 
lenses available for test. ‘This general formula has not the validity of a law of 
nature ; it is a concise expression of the state to which lens design had arrived by 
about 1943, and no doubt departures from the general approximate rule which ha 
been found will be more and more frequent in the future as the art of designing 
lenses progresses. Already some notable advances had been made in this direction] 

The general formula connecting the resolving power with the aperture, angl 
of view, and focal length which was arrived at by this investigation would enabl4 
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an estimate to be made of the best lens to be used under any given conditions, 
provided that the exposure time was sufficiently short to allow the effect of move- 
ment of the camera during flight to be neglected. But as soon as it was seen that 
the general formula in question was a fair representation of the performance of 
lenses on the ground, a laboratory investigation was begun on the effect of the 
forward motion of the aircraft upon the definition of negatives.* By adopting 
the general formula arrived at in the survey of lenses (with some limitations), 
using certain results obtained in subsidiary work concerning the variation of 
granularity with speed of films, and the information embodied in exposure tables 
worked out by the Royal Aircraft Establishment, it was found possible to express 
the definition obtained in aerial photographs as a function of all the “‘ practical ”’ 
variables entering into the problem, such as speed and height of aircraft, speed of 
film, aperture of lens (itself determined by the time of exposure and intensity of 
illumination of the ground), angle of view, and focal length of the lens. This 
formula is capable of use in so many different ways, according to the information 
required, that it has not been subjected yet to exhaustive practical test. One 
interesting feature of this particular work, however, was that it suggested that the 
best time of exposure, for anything approximating to normal conditions, was that 
already in use by the Royal Air Force. It is not certain whether this result should 
be regarded with most satisfaction by the laboratory workers or by those who 
arrived at the exposure by severely practical considerations. 

In addition to the above work, several investigations of interest have also been 
carried out. ‘These relate to the spread of light in the emulsion layer caused by the 
scattering from grain to grain, acheck onthe theory of resolving power of emulsions 
for the type of test-object used in the lens tests, and the nature of the correlation 
between speed and granularity of films for air photography. An investigation has 
also been made on the effect of obliquity of the incident light on the resolving 
power of an emulsion. Finally, the work on the interaction between the image 
given by the lens and the film, which leads to the observed resolving power, has 
been pursued to a state at which we can say with some confidence how the inter- 
action comes about and what determines the resolving power in the negative. 


§2. RESOLVING POWER AS A MEASURE OF PERFORMANCE 


In the early stages of the work doubts were expressed whether measurements 
of resolving power could be used satisfactorily as a measure of performance. 
These doubts appear to have arisen, in the main, from acquaintance with the 
results of resolving-power tests using test-objects of high contrast. It is known, 
for instance, that the visual resolving power of a telescope, for objects of high 
contrast, may be a misleading guide to the state of correction of the lens. The 
same criticism can also apply to photographic tests, especially if, in addition to the 
contrast of the test-object being high, a fine-grain high-contrast photographic 
material is used. In such cases also the exposure is very important. This may 
be demonstrated by a modification of Jones’s tone-reproduction diagram (1921). 
The curve A in figure 1 (a) represents the intensity distribution in the image of a 
point source of light, such as might be obtained from an imperfectly corrected 
lens. In figure 1(b) this has been plotted on a logarithmic scale of intensity in 


* This investigation is not described in this paper ; it was carried out by Mr. J. M. Gregory, 
and will be published later. 
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the lower right-hand quadrant. By projection on to the characteristic curve of th} 
photographic material (shown in the upper right-hand quadrant) the rou 
density distribution may be obtained as shown in the upper left-hand quadran} 
The full radius of the disk of light produced by the lens is 5 units, whereas und) 
the particular conditions of exposure chosen for illustration, the full radius of | 
image is only 2 units. Such a possibility is not unknown in practice. Sharpne | 
beyond expectations, based only on the state of correction of a lens, may be obtaine}} 
on suitable material by careful attention to exposure. | 

In order to explain how it is that resolving-power tests can be an adequate gul 
to performance in many cases, and in aerial photography in particular, it 1s cor) 
venient to use the theory of resolving power mentioned in §1. That theory, te h | 
stated precisely, requires that the test-object be in the particular form of lon} 
parallel lines across which the brightness varies as a sine-curve. ‘The main effed} 
of lack of sharpness is to decrease the contrast between light and dark parts of th 
image, in such a way that the smaller the scale of the image the greater is t 
reduction of contrast. If the contrast in certain details is above a certain minimur}} 
then those details will be seen, but if below, they will not. In more definite terms} 
if the detail is in the form of the above sine-curve test-object and the contrast i | 
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Figure 1. 


(2) Schematic representation of intensity across image of point source as given by imperfed§ 
lens (Curve A), Curve B represents the intensity distribution when the point source 
surrounded by a uniformly illuminated field. 

(6) Tone-reproduction diagram showing detivation of density distribution in the develope} 
photographic material. 
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measured in suitable fashion, the curve of contrast against the logarithm of t | 
reciprocal of the size of detail will follow a course such as that shown in figure 2 
The appropriate measure of contrast is log P, where (1+ P)/(1— P) is the ratio o 
maximum to minimum brightness in the image. The value of log P is approx 
mately equal, for values of log P up to 1-8, to logy + loga+log P,+ log P,, wher 
y is the slope of the curve of density against log exposure for the photographi 
material, (1 +a)/(1—a) is the ratio of maximum to miminum brightness in t 
test-object, log P, is a function of d (the separation between neighbouring pairs aj, 
dark lines in the image of the test-object) and is determined by the image-forming 
properties of the lens, and log P, is a similar function of d which describes how thy 
contrast is reduced by the scattering of light in the emulsion layer of the photaj) 
graphic material. | 

If the aerial image is being viewed, P, is unity and y is unity, there being n 
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photographic material involved, and the threshold contrast is represented by a 
straight line parallel to the axis of log 1/d at a value of log P of approximately 73: 
The threshold for the photographic case is represented approximately by a straight 
line at a slope of 3/4,* the ordinate of which at log (1/d) =0 is equal to log Gt — 2-0, 
where G is the granularity of the photographic material in density microns. 
The reason for the increase in the threshold contrast for smaller detail is that the 
granularity becomes relatively more obtrusive and disturbing to recognition of 
detail the finer the detail. 

Where a log P curve intersects a threshold line, the contrast of the image is 
equal to the contrast threshold of the receptor (either the eye or the photographic 
material) and the value of 1/d at this point is equal to the resolving power. Thus in 
figure 2 the intersection of the log P, curve with the threshold contrast line for the 
eye gives the visual resolving power of the lens ; the intersection of the logy + log P, 
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Figure 2. Contrast in visual and photographic images, and thresholdsfof contrast plotted 
against reciprocal of size of image (logarithmic scales). 


= with the threshold contrast line for the film gives the film resolving power, 
_and the intersection of the log y + log P, + log P, curve with the threshold contrast 
line for the film gives the resolving power of the lens and film together. Fora set 
of test-objects of different contrast, the log P curves will consist of a set all of the 
same shape (namely that of the curve log P,+log P,+logy), but displaced from 
each other vertically by the differences in log a for the test-objects. If, therefore, 
we are prepared to make resolving-power tests on a number of test-objects of 
different contrast, we are able to reconstruct the curve of log P,+log P,+logy. If 
the observations are visual the reconstructed curve will be log P,; if the obser- 
vations are photographic a knowledge of the y-value and of the light-scattering 
properties of the photographic material will enable log P; to be obtained from 

* These figures are derived from those given by Romer and Selwyn (1943) and apply therefore 


to Cobb test-objects. It is known, however, that approximately the same values apply for sine} 
Curve or other test-objects. 
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log P,+log P,+logy. Itis assumed also, of course, that the granularity of the fil i 
is known, otherwise the appropriate threshold line cannot be drawn. Log P, i | 
determined only by the optical properties of the lens; it can be calculated from th¢ : 
light distribution in the image given by the lens of a point source. 

Although the method of presentation may be unfamiliar, a set of measurements} 
of resolving power for test -objects of different contrast is equivalent, therefore, to 4} 
description of the image-forming properties of the lens, and in that respect is no 
markedly inferior to any other description. Indeed, it is difficult to see what othe1 
measurable index of performance could be used: at any rate, no other has beer 
suggested. 

The diagram in figure 2 also indicates how far above the threshold contras} 
the actual contrast is, and to that extent is a measure of the visibility of objects 
which are greater in size than at the resolution limit. Now in aerial photography} 
for reconnaissance, the quality of a lens is determined, very largely, by the Ho 
smallest details visible. If details of a certain size and contrast can be recognized) 
larger details of the same contrast, or details of the same size and greater contrast} 
will certainly be recognized, and it is not of much importance whether they can A | 
recognized with a certain facility or whether they can be recognized with a littl 
greater facility. In other words, a resolving-power measurement is likely to b 
an adequate guide to the performance of an aerial camera lens. Before leaving{ 
this discussion it is worth interpolating several conclusions to be drawn from a . 
not directly connected with the immediate problem of estimating lens performance} 
It will be seen that it is not possible to derive much information about the photo 
graphic resolving power from the visual resolving power for a test-object of give 
contrast, or vice versa, even if the properties of the photographic material are} 
known in detail. If we know only the resolving power of the film, measured, foa} 
instance, by a test-object in contact with the emulsion layer, the situation is eve 
worse. Nevertheless, from a knowledge of the general properties of the curves 0 
log P, and log P, against log (1/d) rough estimates or limits of resolving power ca 
be obtained from limited data. 


Choice of test-object 


So far the discussion of the possibility of using a resolving-power measuremenj] 


statements are concerned, on the assumption that the test-object is of the sine-curvaf!_ 
type. But, in fact, detail of such a character must be rare—the only exampl 
approximating to it which comes to mind is a corrugated iron roof upon which the!) 
sun is shining not too obliquely. ‘The question therefore arises whether such 4). 
test-object can be used for testing the capacity of the lens to resolve details of othea}_ 
types, or whether we must seek other test-objects approximating more closely taf| 
those which arise in practice. In the latter case the above theory must be used only}! 
as an approximate indication of general trends. There is already some evidence ii 
the literature that the value of resolving power obtained for a photographic materia} 
varies with the form of test-object employed. Ross (1924) found that a test-objeci 
consisting of a number of parallel lines separated by a distance equal to their width 
gave the same values of resolving power as a ‘“ double-star”’ test-object consisting})) 
of pairs of disks, but higher values than those given by a fan-shaped test-objec 
Sandvik (1927) investigated the variation of resolving power with the ratio of widt 
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of line to width of space, and found that the resolving power fora given photographic 
material increased uniformly with the logarithm of the line/space ratio. In 1928 
he showed experimentally that the resolving power increases with the contrast of 
the test-object. The theory to explain these two results put forward by Frieser 
(1935, 1938) was based on the sine-curve test-object, and it was through his work 
that the simplification introduced by using a sine-curve test-object was first 
appreciated. 

Some early experiments were made with a test-object (shown in figure 3) which 
is fairly typical of the parallel-line type. This was photographed by means of an 
8-inch F'/2-9 lens stopped down to F/8, using different films, at different exposure 
levels and different focus settings. A study of the negatives showed that the 
length of the lines of the test-object influenced the observer’s estimate of the 
smallest group which could be resolved. For example, there was usually a gradual 
transition from groups which could be recognized clearly to groups which could 
not be recognized, so long as the resolution limit occurred at some group within a 
set of five of constant line length. When, however, the transition from resolved to 
unresolved groups occurred at a group near X in figure 3, where there is a sudden 
change in line length, it was found that the larger group on one side of X was 
clearly defined, whilst the smaller group on the other side of X showed no trace of 
line structure. The number of lines also influenced the result, for if only a small] 
portion of the group at the resolution limit were exposed to view, the observer did 
not recognize it as resolved, but when more of the area was exposed to view the 
observer had no hesitation in deciding that the pattern was recognizable. The 
results obtained with the test-chart described by Sayce (1940), in which the lines 
are all of the same length but with the separation between them diminishing from 
end to end of the chart, are not likely, therefore, to bear any simple relation to those 
obtained with other charts. 

Experiments were also carried out using letter charts of the type used by 
ophthalmic opticians, for the reason that the letters would presumably be unknown 
to the observer, and thus provide a better test of recognition than test objects the 
shape of which was already known to the observer. Such letter charts were, 
however, not very satisfactory, for they did not provide a very uniform scale of 

detail, and in any case the observer soon memorized the sequence of letters. 
There seemed little doubt from the somewhat scattered evidence that different 
resolving powers would be found with different test-objects, but this would be of 
no importance if the resolving powers yielded by different test-objects always bore 
the same ratio to each other. To check this point, comparisons were made 
between the parallel-line test-object shown in figure 3 and a test-object which was 
designed to eliminate the defects of the parallel-line type. This new object 
contained the minimum number of lines necessary for a determination of resolving 
power (namely, two) and in which the ratio line-length/line-width was constant at 
all sizes. Moreover, the length of line was no greater than was necessary to 
permit the recognition of the longer dimension, and therefore the direction, of the 
“test-object. A test-object described by Cobb and Moss (1928), consisting of two 
bars of length three times their width, separated by a distance equal to their width, 
had hitherto been used for visual acuity measurements; its dimensions fitted the 
above desiderata sufficiently closely to suggest the name “ Cobb test-object”’ 
for this type. The experiment was carried out with a composite chart consisting 
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of groups of Cobb test-objects and groups of parallel line test-objects. Phot 

graphs were made at different settings of focus of an 8-inch F/2-9 “ Pentac”’ lens 
and the limit of resolution was assessed by a number of observers. ‘The resolving}, 
power limit as determined by the Cobb test-object was plotted against that deted} 
inined by the parallel-line test-object, as in figure 4. Each point represents th]] 
mean estimate of a number of observers. (The departure of some of the point 

from the smooth curve is probably due to the fact that in the parallel line test-obje¢ 
the line frequencies of consecutive groups were in the ratio of 1-2 to 1, whereas if} 
the Cobb test-object the ratio was 1-1 to 1.) It will be noted that although 4 
the position of focus the two test-objects give substantially the same value 
resolving power, there is an appreciable departure from equality on each side d 
the focus ; at positions inside the focus, the parallel-line test-object leads to highef} 
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Figure 4. Comparison of parallel-line and Cobb test-objects. 


values of resolving power than the Cobb test-object, whilst beyond the focust 
reverse is true. ‘I'he explanation for these results must be that the distributioa} 
of light in the image given by the lens are not the same on both sides of the foc 
position. ‘The effect of this change must then be different in the image of ‘tl 
multi-line test-object from that in the image of the two-line techie 
Obviously the resolving powers for different test-objects do not, in general), 
bear the same ratio to each other, and if resolving-power tests are to be of value th } 
test-object must be chosen to simulate, as near as may be, the type of objects to t 
photographed, It was therefore decided to adopt the Cobb type of test-objed} 
which is probably as nearly similar to the detail in air photographs as can 
obtained in a simple geometrical pattern. Although no comparisons with oth 
test-objects for accuracy have been made, it is possible that this type of test-obj 
does not yield results of quite so high precision as others. By using a number. 
Cobb test-objects disposed at random, and plotting the number of a given sii 
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Figure 3. Parallel-line test-object. 
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Figure 5. Form of test-object adopted for resolution tests. 
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recognized against the size an estimate may be made of the size at which half of 
each group are recognized, and the precision may be thereby somewhat improved 
(Romer and Selwyn, 1943). This method could not be applied to lens examination 
mainly because the complete test-object must be kept reasonably small. In the 
form used for the present work the test-chart consists of groups in aspiral arrange- 
ment, as shown in figure 5. Half the test-objects are vertical and half horizontal, 
so that measurements can be made along radial lines and along tangential lines. 
The definition of the lens used in preparing the test-chart and of the collimator 
used for the tests is likely to be worse for points off axis, so the smaller groups, 
the sharpness of which will be relatively more affected by poor definition in these 
two lenses, are placed near the centre. 

Aerial photographs are different from ordinary photographs in that the 
atmosphere itself appears luminous and throws a veiling light over the surface of 
the earth, as a result of the scattering of the light by small suspended particles, 
Thus while an average subject photographed from the ground has a brightness 
range from the brightest part of the subject to the darkest of 160: 1 (Jones, 1941), 
aerial photographs show that the ratio for subjects seen from moderate and high 
altitudes is only 10:1 or so. Detail can range in ratio of maximum to minimum 
brightness from 1:1 to 160: 1 for photographs taken on the ground, and to 10:1 
for aerial photographs. A rough mean for the detail in the terrestrial photograph 
may therefore be taken as \/160:1, and for the aerial photograph as 1/10: 1, 
say 3:1. Examination of details in aerial photographs larger than that just 
resolved, so as to avoid the reduction in contrast due to lack of sharpness, suggested 
that a much more nearly typical value of ratio of brightness was 1-6: 1. In terms 
of the usual photographic units, therefore, the appropriate density difference in 
the test-object was estimated at 0-2 (i.e. log1-6). Figure 5 has been reproduced 
to approximately this density. It is possible that for testing lenses for use in 
ordinary photography on the ground a somewhat higher contrast than this is 
desirable. 

Suppose we now return to the example shown in figure 1, and consider what 
happens when the contrast is decreased. ‘To decrease the contrast of a point 
source we must place around the point source an extended source of uniform 
brightness (as occurs for instance in astronomy where the stars are surrounded by 
the general illumination in the night sky). The image given by the lens then 
becomes as shown in curve B of figure 1 (a) and the photographic image as in curve 
B of figure 1(b). It will be seen that when the contrast is low the photographic 
image of the point source extends to the full diameter of the optical image. ‘Thus 
we may conclude that changes in the exposure are likely to affect the image less if 
the contrast is low and that the outer and weaker parts of the optical image are 
important when the contrast is low, but may be unimportant if the contrast is high, 
The latter point is, in fact, usually true whether the image is observed visually or 
photographically. Suppose that the images of a point source given by two lenses 
have intensity distributions as shown in figure 6 at A and B. A represents an 
image of moderate diameter, into which all the light from the lens comes to a focus. 
B represents a central bright image surrounded by a flare of 1/100 of the intensity, 
covering an area 100 times that of the bright core; the amount of light in the core 


is thus equal to the amount of light in the flare. The log P, curves corresponding 


to these images are shown at A’ B’ of figure 6. The visual and photographic 
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threshold lines are annotated ‘“‘ Visual High Contrast” and “‘ Phot. High Contrast | 
If the contrast of the point object is now reduced by surrounding it by a uniformly} 
illuminated field, the log P, curves A’ and B’ should be shifted downwards by a | 
appropriate distance but, for convenience, the threshold lines have been shifted} 


upwards by the same amount. They are annotated “ Visual, Low Contrast’’ | 


and “ Phot. Low Contrast” respectively. In order not to complicate the argu 


ment, the effect of the spread of light in the emulsion layer is neglected. It will bel 
seen that the resolving power when the contrast is high is determined very largely |} 
by the bright central core of the image which suffers from flare, but that as the}j, 


contrast is decreased the flare becomes much more important. We shall be 
referring to this point later in connection with the effect of aperture on the resolving 
power. 


Choice of photographic materials 


It has already been mentioned that the prospects of simplifying the problem} 
by using the theory of resolution were remote. |The work therefore had to be} 
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Figure 6. Effect of flare round image on log P; curves, and effect of this on resolving power 
for high and low contrast test-objects. 


proceeded with on an empirical basis, and for this reason the choice of photographic 
materials was governed by two practical considerations: firstly, it was required ta 
know what standard of performance was obtained with existing equipment, 
and secondly, it was required to know how far this standard might be expected tos 
rise, if the properties of the film could be improved, or if by some means the 
exposure time could be lengthened to permit the use of slower fine-grain materials 
The first consideration naturally led to the choice of the existing Service type films, 
the Kodak versions of which were ‘“‘ Super-XX Panchromatic Aero Film ’? 
and “‘ Panatomic X Aero Film”. The former is the faster of the two and has a 
lower resolving power. Partly in order that the effects of chromatic aberration o 
the lenses should remain the same in all the photographic tests, and also because 
panchromatic sensitization was not likely to be superseded by any other, such as} 
infra-red, it was desirable that the film chosen to satisfy the second reed 
1] 


SETS 


should also be panchromatic. It should also have high resolving power and 
moderate speed. No existing high-resolution film has high speed, but the 


emulsion manufacturer is continually striving to increase resolving power Bee: | 
retaining a useful speed. Allowing for possible advances in this directio | 


| | 
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and also in devices (for reducing vibration and compensating for the forward 
motion of the aircraft) which would permit of longer exposure times, it was 


_ decided to choose Kodak ‘“ Microfile’” Panchromatic film. This has a resolving 
_ power about three times that of Super-XX, whilst its speed is about 1/80 that 


of the latter. If the speed of this film could be increased, and the exposure time 


_ increased, there was some possibility of it being brought into use. 


Development conditions 


Since the resolving power of the film is dependent on the contrast, not only of 
the test-object, but also of the emulsion, it was important that development 
conditions should be standardized. The Service type films were developed to 
give about the same contrast as obtained in the R.A.F. Photographic Recon- 
naissance Units, namely y*=1-4 for Super-XX and y=1-6 for Panatomic-X, 
both films being developed in Kodak developer D.19b. Microfile film, in its 
normal use for microcopying documents, etc., is developed to a high contrast, but 


_ for pictorial photography a high contrast in the negative would lead to a density 


\ 


range too great to be printed on existing positive materials. Development was 
therefore standardized to make y=1-6. Kodak developer DK.20 was used 
initially, but this was later replaced by Kodak developer D.76d, which was more 
conveniently available, and conditions were adjusted to produce the same negative 
characteristics as formerly. 


Description of apparatus 


The examination of lenses of focal lengths from 5 inches to 20 inches was 
carried out with the apparatus shown diagrammatically in figure 7. The test-object 
(B) was mounted at the focus of a collimator objective (A), and illuminated by 
light from a tungsten filament lamp. An opal diffusing screen (C) to ensure 
uniform illumination, and a minus-blue filter (as used in R.A.F. practice) were 
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Figure 7. Arrangement of apparatus. 


mounted between the source and the test-object. The test-object was exposed 
by means of a roller blind shutter (E) attached to the lamp-house (F) and a slide 
(D) carrying neutral filters, graded so that the brightness of the test-object could 
be altered in steps by a factor of 1/2, was fitted between the lamp-house and the 
shutter. The collimator objective was a Wray achromatic doublet of focal length 
60 inches and aperture F/15, permitting lenses of aperture 4 inches to be tested. 


* In this connection y is the slope of the straight portion of the characteristic curve of the 
material. 
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The camera incorporated a lens holder (K), capable of carrying lenses of || 
external diameters up to 6 inches, and a film holder (G) sliding on rails, guided by | 
two V-grooves and a knife-edge. The film (in strips 35 mm. wide) was held | } 
against a register glass of thickness 0-15 inch, in conformity with the standard || 
register glass used in the smaller R.A.F. cameras. ‘The length of the film holder 
was such that a total angle of view of 56° could be covered with a 20-inch lens. | 
The camera could be rotated about a vertical axis M the position of which was) 
adjusted to coincide with the front surface of the lens to be tested. ‘This was to | 
ensure that the lens did not move out of the beam of light from the collimator. 

A larger collimator and camera was constructed later to accommodate lenses of | 
focal lengths up to 50 inches. The collimator lens was a Cooke “ Photo Visual ” | 
objective of focal length 134 inches and aperture 7 inches. | 

Both collimator objectives were well corrected, and their influence on the 
results obtained with different aircraft lenses was negligible. The lateral aber- 
rations of a collimator objective of focal length F are reduced by a factor f/F in | 
the image plane of a test lens of focal length f, and the longitudinal aberrations 
by a factor f2/F2. A number of 20-inch lenses, examined with the aid of both 
collimators in turn, showed the same values of resolving power, within the limits | 
of experimental error. 


§3. EXPERIMENTAL PROCEDURE 
Exposure conditions 


The illumination in the focal plane of a lens is not uniform, but decreases as ]]} 
the angular separation from the axis is increased. In a thin lens the effective 
aperture is reduced as the angle of incidence is increased, and the principal ray 
strikes the focal plane at a greater angle the greater the angle of incidence. Thus 
the intensity of illumination J at any angle is given by J=/, cos! 6, where J, is the |] 
intensity of illumination at the centre of the field. If the lens is of finite edge 
thickness, or, more particularly, consists of separated components, an additional 
decrease in illumination occurs, due to “ vignetting’. Absorption and reflection 
losses due to the register glass, as used in R.A.F. cameras, also increase with 
angle of incidence. Minor variations with angle arise from reflection and 
absorption losses in the different components. 

The exposure given to the film therefore decreases with increasing angular | 
separation from the axis. Since the resolving power of the film varies with the ) 
exposure it receives, the resolving power of the lens and film used together will 
vary with angle due to film properties as well as to lens properties. In addition, 
the subject may cover a considerable range of brightnesses: in aerial scenes this 
range is comparatively small, as stated above. In order to make a complete | 
survey of the problem with regard to these practical considerations, it would be 
necessary to collect data on lens-film resolution for a number of different exposure 
levels. It was thought, however, that such a procedure would unnecessarily | 
lengthen the programme of experimental work and complicate analysis of the data | 
and all variations in resolving power for a given lens-film combination eord | 
restricted to those due to the lens alone, by adopting one exposure only. The | 
exposure chosen was that which would result in the maximum film resolution — 


which occurs, for the most negative materials, when the density of the negative is 
about 1:0. | 
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Measurements of the relative illumination at different positions in the field of 
_ each lens were made in order to discover how serious the variation was, and whether 
| the technique of giving constant exposure to all parts of the film was justified. 


Measurement of relative illumination in the focal plane 


The relative illumination in the focal plane of each lens was measured by 
, photographic photometry in the following way. Using an illuminated disk in the 
_ place of the test-object, exposures were made, of constant duration, at 0° and at 5° 
intervals on each side of the axis, up to the full angular extent of the test-lens, on a 
. film of moderate contrast. This was repeated for different apertures of the lens, 
| and the exposure level for each aperture was adjusted so that the densities obtained 
on development lay on the linear portion of the sensitometric curve of the material.* 
. Sensitometric film-strips were developed together with the test-strips exposed by 
_the lens, and the values of log exposure corresponding to the densities produced 
_by the lens were read from the sensitometric curve thus obtained. Since we 
_ require to know only the relative values of the intensity at different points in the 
field, these log exposure values may be considered to be values of log intensity. 
In order to allow for possible variations in the shutter speed, and the intensity 
of the light source, two exposures were made (by moving the film laterally inside 
the film-holder) at each position, and, therefore, at all except the axial position there 
were four exposures for each angle. It was found that there was little variation in 
_the densities so obtained, and the values of log intensity were usually reproducible 
to within about +0-02. Greater precision than this could have been obtained by 
controlling lamp-voltage and shutter speed, but was considered unnecessary. 


Measurement of visual resolving power 


Although the research was to be primarily for the purpose of assessing the 
photographic performance of lenses, preliminary experiments had shown that 
some useful information might be gained if the resolving power of the aerial image 
were measured. The design of the apparatus did not permit visual measurements 
to be made easily for extra-axial images, and these measurements were confined 
to the axial image. Test-objects were made on Kodak “‘ Maximum Resolution”’ 
Plates, to the same design as for the photographic test-object. The plates 
referred to were the only type available for recording faithfully the number of 
lines per mm. required for the visual measurements. _‘T'wo test-objects were made, 
one having clear lines on an opaque background, and the other having clear lines 
ona background of density 0-2. The high contrast test-object was used in order 
‘to measure the resolving power of a lens in the way in which this measurement is 
usually carried out, and the low contrast test-object to provide measurements of 
resolving power approaching more closely those obtained photographically. 
| The test-objects were placed at the focus of the collimator objective, and a 

microscope was used to observe the image produced by the lens under test. ‘The 
focus of the microscope was adjusted until the maximum number of groups of the 
test-object was resolved, and the number of the smallest resolvable group was 
recorded. The relative aperture of the microscope objective was always greater 
than that of the lens under test. 


* It is not essential to the method that all the densities on the test negatives lie on the linear 
portion of the characteristic curve of the material, but this ensures the same degree of accuracy 


i 


for all values. 
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Photographic resolution tests 


After preliminary trials to ascertain what region of focus was to be explored} 
and to find the exposures necessary to give a density of about 1-0 at each positio iy 
in the field, the film-holder was set and exposures were made on the axis and at 5°} 
intervals up to the full angular extent of the lens, on both sides of the axis. This| 
procedure was repeated at different settings of the film-holder, until sever) 
planes, passing through the region of focus, had been explored. This series was 
then repeated at all apertures of interest, and for the other films. | 

Since most of the lenses were examined at about 5 apertures, this procedure 
resulted in about 100 negatives per lens, with 20 or more images on each negative. | 
The amount of work involved in measuring the resolving power in the negatives, 
not to speak of the complexity of the problem of analysing the data, imposed 
limitations on the thoroughness with which the inspection could be carried out. 
It was therefore necessary to limit the measurements to one negative at each 
aperture, for each of the three films. 

For this reason, from each set of 7 negatives, that one was chosen in which the 
best axial definition was obtained. The only other plane of “focus” of any 
significance would correspond to the negative in which the best average definitio 
was obtained, having regard to the entire area of film which the lens was intended ta}f 
cover. The selection of this plane wouid involve the inspection of all the images 
however, and would not solve the problem. ‘This method was used later ir] 
connection with the focusing of lenses in cameras. It is worth pointing out tha 
in a lens afflicted with spherical aberration both the plane in which the best axial} 
definition occurs and the optimum focal-plane shift away from the lens as thef| 
aperture is reduced. ‘Thus in setting cameras to give the best performance at 
different apertures without re-focusing, it is necessary to compromise on a focal} 
setting which will be the optimum with respect both to the range of apertures and 
to the angle of view to be covered. Having chosen the negative, the number off 
the smallest group of the test-object recognizable at each angle was recorded and|p 
converted into resolving power in lines per mm., applying a correction discussed 
below. The resolving powers for radial and tangential lines were recorded} 
separately. | 


between the two methods. ‘The National Bureau of Standards apparatus has a| 
number of fixed collimators accurately positioned so that measurements o 
distortion may be made. A high contrast test-object is used in which the number 
of lines and the ratio of length to separation is not the same for all line frequencies, 1 | 
The photographic material used is the Eastman “Type VB Spectrographic’”| | 
plate, which has the same order of resolving power as “ Microfile” film. The 
problem in hand at the National Bureau of Standards is to calibrate lenses in terms 
of the optimum focal length having regard both to resolution and distortion, and} 
is, therefore, different from the problem with which this paper is concerned. 


Measurement of resolving power in negatives 


In determining the limit of resolution from the negative, the conditions of 
viewing were first standardized. ‘The illumination of the field of the microscop 
was kept fairly constant from negative to negative, and the degree of magnification: | 


t 
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was kept constant for a given photographic material within the obvious limits that 
the magnification should not be so small that the size of test pattern is near the 
limit of the acuity of the eye, nor so large that the obscuring power of the graininess 
obliterates the image. It was found that a convenient magnification was 20 
for the Service type films, and this had to be increased to 75 when observing good 
images on “ Microfile” film. . Since all the negatives were of approximately the 
same density, the conditions of illumination were kept constant, once a satisfactory 
level of illumination had been found. It is an advantage of the low contrast test- 
object that glare is absent when viewing the negative. The “ spurious” resolution 
effect referred to by Gardner and Case (1937) has not been observed by the authors 
with the low contrast Cobb test-object. This effect, which is observable with 
multi-line high contrast test-objects, is that lines of a particular size appear to be 
resolved, although closer inspection reveals more lines than appear in the original 
test-object. With the Cobb test-object one should observe a strong line in the 
middle and possibly weaker lines on the outside if this effect occurred, but in fact, 
when the contrast of the test-object is low, the effect cannot occur. 

The measurement of resolving power is to some extent subjective and is liable 
to errors which may appear large by physical standards. It is usually necessary 
for the observer to set himself some criterion of resolution, with the result that 
some observers judge as resolved a pattern which other observers would consider 
inadequately defined. ‘These differences occur largely amongst inexperienced 
observers, and it has been found that when a number of workers have gained 
experience the agreement between them is good. ‘The resolving-power determin- 
ations were carried out in the present work by one observer, who adopted the 
criterion that a group of a given size was resolved if the Cobb element appeared to 
consist of two separate patches, each quasi rectangular in form, and if the direction 
in which the “lines” were oriented agreed with that on the original test-object. 
It is sometimes found that the pattern appears to consist of two lines pointing in a 
direction at right angles to that in the original test-object. This is due to fortuitous 
“clumps”? of grain. 

L. C. Martin and C. A. Padgham* have investigated the precision of photo- 
graphic resolution measurements under conditions similar to those described in 
this paper, and their results agree with estimates made by the present authors 
which suggest a standard deviation of about 10% for repeated measurements on 
the same negative. It was found, in addition, that the percentage precision of 
measurement was the same for both high and low contrast test-objects. 

In converting the number of the smallest groups which can be resolved into re- 
solving power in lines/mm., a correction must be applied for the increase in magni- 
fication for extra-axial images. In practice the resolving power was calculated from 
the dimensions of the test-object rather than by measurements on the negative. 
If 1/R, is the separation between lines in the object, f’ the focal length of the 
collimator, f the focal length of the test lens, and 6 the angular separation from the 
axis at which the image falls, then the line frequency R in the negative is given by 


Rp=(RpF cos )/f 
and Ry =(RoF cos’ #)/f, 
where the suffixes R and T indicate radial and tangential lines respectively. 
* Private communication. 
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As a result of this distortion of the test-object the length of radial lines becomes||_ 
greater than 3 times the width of line whilst the length of tangential lines becomes) 
less than 3 times the width of line by a factor cos. This change of shape is}| 


however, small enough to be of little importance. | 
i 


§4, EXPERIMENTAL RESULTS 


‘The ieuses examined in this survey were mainly of British manufacture | 
although a few American and German objectives were included. The focal 
lengths ranged from 5 to 49 inches and their maximum apertures from F/2-9 to F/7 | 
‘The types of construction included wide and narrow angle, triplet and telephote 
types but in spite of this variety there was sufficient similarity in their performance 


to enable results to be calculated as typical of the majority for presentation below. 
| 
| 


Variation of resolving power with aperture 


In figure 8 resolving power is plotted against /. No., both on logarithm 
scales, for six of the lenses tested. The results are quite typical. The measure 


ments refer to the axis only, but include visual measurements with high and low) 
contrast test-objects, as well as photographic measurements on the three films|}, 


with low contrast test-objects. The broken line indicates the value of resolving 


power calculated from the formula 1/A(F. No.), A being taken as 0-55. This is aj 


customary value to adopt for a perfectly corrected lens (Conrady, 1929). Strictly), 


it refers to visual observation of two close point-sources seen against a non 
luminous background. It will be seen that except at the highest apertures the 
observed resolving powers for high contrast Cobb test-objects approximate t | 
those given by the formula, being usually some 10°% lower. The greater departures 
at the higher apertures must be due to uncorrected aberrations. With lo 
contrast test-objects the visual resolving powers are inevitably less than for high} 
contrast test-objects, but the deviations from the formula at large apertures are} 
also so considerable that the curves of resolving power against F. No, display 4 
definite maximum. The photographic resolving powers are much lower than thd 
visual, and diminish in the order Microfile, Panatomic-X, Super-XX. They alsa 
display maxima in the curves of resolving power ‘against F’. No. and the maxima 
tend, though perhaps not in strongly marked fashion, to larger values of F’. No. ix 
the order Microfile, Panatomic-X, Super-XX. 

It is not difficult to account for these results in approximate fashion. To dd 


this we shall use the theory already mentioned, which, although strictly applica 


only to sine-curve test-objects, may be expected to give a qualitative account of the} 
results to be expected from other types of test-object. It is first necessary td 
make some estimate of the curves of log P, against log (1/d) for different apertures 
and for this purpose it is convenient to note that if the lens were perfectly correcte 4 


(and if the light from the test-object were monochromatic—though this does not}) 


make an important difference) the curves for log P, for different apertures would 
be all of the same shape but separated from each other in the log (1/d) direction by 
an amount equal to the logarithm of the ratio of successive F-numbers. A set of 
such curves is shown by the light lines of figure 9. 

When the path differences introduced by imperfect correction are less. thar 
half a wave-length, the diameter of the central nucleus in the image of a point| 
source is no greater than when correction is perfect. There is, however, some| 


| \ 


| 
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Figure 8. 


Typical curves showing variation of visual and photographic resolving power 


with aperture (for axial region only). 


(a) Visual, (b) Visual, (c) Microfile Film, (d) Panatomic X film, (e) Super XX film. 
Curve (a) refers to high contrast test-object, and curves (b) to (e) to low contrast test-object. 
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Figure 9. Log P; curves for perfectly and 
imperfectly corrected lens at different 
apertures. ‘Thin lines are curves for per- 
fectly corrected lens and thick lines for 
imperfectly corrected. 
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Figure 10, Estimation of log Pe curves from 
resolving power of films. 
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loss of light, and the light which is lost from the nucleus appears in the rings| 
surrounding the nucleus, usually as a general luminous halo. (Martin, 1926 ;| 
Conrady, 1929). An illustration of the effect of a halo of this sort is provided by 
figure 6, which may thus be used as a guide for estimating how the curves in light ul 
lines of figure 9 must be modified when a small amount of aberration is present, | : 
When the amount of aberration is large the central nucleus, if it still persists, 
becomes larger than in an aberration-free image at the same aperture. T he | 
nucleus may not be evident in some cases, but there also the effective size of the 
image is greater than that of the corresponding aberration-free image. ‘Thus when, | 
the amount of aberration is considerable the log P, curve necessarily lies sub-}f} 
stantially to the left of the curve for a perfectly corrected lens of the same aperture, . 
even at low values of log P,. When the amount of aberration is small, on the|}} 
other hand, as indicated above, the log P,; curve remains of approximately the same! 
shape as that for a perfectly corrected lens, but is shifted downwards slightly. | 
Now the general result of decreasing the aperture of a lens is to reduce the path ! 
differences at the focus. Thus at high apertures the log P, curves will, if the} 
aberrations are considerable, lie substantially to the left of the curves shown by} 
light lines in figure 9, while for small apertures they will not differ much from these} 
latter curves, but will tend to lie a little below them. ‘To add more precision it has} 
been assumed that at F/45 the lens is practically free from aberration, while at F 4 
a guess at the shape of the curve has been based on the values of the visual resolving}} 
power obtained on lens No. 137006 for test-objects of high and low contrast. The} 


curves for the apertures intermediate between F/4 and F/45 were drawn by eye toll 
form a consistent family. ‘The whole family of curves thus obtained is shown by| 
heavy lines. It should be noted that a family of real log P; curves might not 
display such an obviously systematic variation as those shown in figure 9. 
We now require to estimate the influence of the diffusion of light in the emulsio 1 : 
layer. Fairly satisfactory results have been obtained, in other connections, by| 1 
assuming that the intensity of illumination in the emulsion decreases mars | ; 
tially with the distance from a long, narrow, brightly illuminated strip of the} 
emulsion. If this is true, then the curve of log P, against log (1/d) has a definitelf| : 
shape which is the same for all photographic materials but which lies toward | 
higher values of log (1/d) the less the diffusion of light in the emulsion layer.}} 
The actual position of the P, curve can be located on the log (1/d) scale, once the} 
threshold line (determined eS the granularity and y-value of the material) has been: H| 
fixed, and the resolving power measured. ‘This is illustrated in figure 10. The] 
resolving powers of the three films were measured using a contact test-object,|]] 
with a density difference of 0-2 between the light and dark parts. It should not be 
concluded that the diffusion of light is exactly the same in Panatomic-X as in 
Super-X-X—the calculation is not a very exact one. | 
Adding together the log P curves of figures 10 and 11 gives the curves off 
log P,+log P,, and these are shown in their upper range in figure 11 for Super- | 
XX and Panatomic-X, for which they are conveniently identical. Figure 11 also 
shows the curves of log P; over the lower range: for low values of log P only visual 
resolution is possible, and it is not necessary to include the log P, + log P, curves., 
The resolving power at any aperture may be read from these curves. When 
plotted, the variation of resolving power with aperture is as shown in figure 12,, 


which also includes values for Microfile film. The qualitative resemblance to the: 
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curves shown in figure 9 will be obvious, and in the main phenomena shown by 
photographic lenses at the centre of the field may be regarded as accounted for. It 
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Figure 11. Log P; and log P; +log P, curves and threshold curves for imperfect lens 
shown in figure 9. 
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ANGULAR SEPARATION FROM AXIS (DEGREFS) 


Figure 13. Typical [resolving power, angular separation from axis] 
curves. Full curves: radial lines. Dotted curves: tangential 


lines. 
should be particularly noted that the results of the tests are inconsistent with any 
assumption except that at full aperture important amounts of aberration exist at 
the centre of the field. . 


| 
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Although there is no published evidence that the variation of photographi«| 
resolving power with aperture had been investigated, either theoretically 03 
experimentally, before the work described inthis paper was commenced, there had} 
been some discussion in photographic literature of the fact that the resolving powej}) 
of a lens when used with photographic materials is much less than the figure usual ! 

ccepted as the theoretical limit. As the discrepancy is particularly great at high} 

apertures, it has been suggested on several occasions (for example, Zeigler, 1936) 
that the reduction in resolving power is due to the wide angle of the cone of ligh | 
incident on the film and the finite thickness of the emulsion layer. Ina paper see | i 
late in the progress of this work Roeder (1941) also subscribes to this view, althoug 
he does, in addition, suggest that at high apertures the lower contrast of the imageé 
has a greater effect upon the photographic than on the visual resolving power: 
The proposed explanation is, however, not correct, as has been shown by Gregory 
(1946). A similar explanation for the poor resolving power observed off-axis has . 
also been shown by Gregory to be incorrect. Pinoir (1944) also discusses the above . 
points, and shows that there is serious quantitative disagreement between the 
proposed explanation and the observed resolving power. It is interesting that in 
the paper quoted, Pinoir discusses the phenomenon of resolving power b | 
arguments verv similar to those we have used. His work was, of course, quite 
independent of ours. 


Variation of resolving power with angle 


Curves of resolving power against angular separation from the axis (hereafter) 
referred to as (R, 6) curves) were obtained for these lenses in the three films] 
In figure 13 is a set of such curves for the 14-inch lens No. EY 118K. The (R, 9}} 
curves for each aperture refer to the plane in which the best axial resolving power) 
occurs at that aperture. The general character of the curves in figure 15 is typica 
of all the lenses so far examined. Individual lenses show minor divergencies off 
shape, but the similarity in the (R, @) curves over a wide range of types of lens isH)) 
surprising. At the particular focus setting adopted, the resolving power is highes i| 
at the centre of the field. At points off-axis there is, except at special points, aif} 
difference between the resolving powers for radial and tangential lines, and the} 
resolving-power curve for radial lines usually has a secondary maximum at someff 
point in the filed. i 

The same general shape of curve is obtained for all three films, the average leve | 
of resolving power being lower when the Service films are used than when Microfile}} 
film is used. ‘This is to be expected because of the greater granularity and the 
greater diffusion of light in the Service films. But it has already been shown tha | 
aberration of an important amount exists in the centre of the field, where the}! 
resolving power is greatest. Consequently it must be concluded that, away from} 
the centre, aberration is present to an even more important degree. Stopping {| 
down the lens is for that reason more effective in improving photographic resolutio 
at points off-axis than on-axis, as may be seen from figure 13. Some curious 
changes, such as in the curves for Microfile film, where the minimum at F/8 is) 
lower than at F/6-3 but increases as the aperture is decreased below F/8, are almost'| 
certainly due to the re-focusing at each aperture. | 

One curious feature of the (R, 8) curves for radial lines deserves some mention. | 
It was found that the secondary maximum and the minimum in the (R, 6) curves: 
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occurred at more or less constant positions in the field. Table 1 gives the tangents 
of the angles corresponding to these positions, and also the tangents of the semi- 
_ angle of view, defined as the angle at which the resolving power falls to zero. It 
will be seen that, with few exceptions, the maximum occurs at about 0-8 of the 
linear field and the minimum at about 0-5 of the linear field. 
In most photographic objectives, the field curvature is such that the sagittal 
and tangential surfaces depart from the Gaussian image plane, and then re-cross 
it at some point towards the edge of the field. It is to be expected, therefore, that 
if the shape of the (R, 0) curve is determined by field curvature, the resolving power 
' should first decrease as the angular separation from the axis is increased from zero, 
and then increase as the image surface swings back to re-cross the “‘ focal plane”. 


Table 1. Positions of secondary maximum, minimum and zero in (R, 9) curves 
for radial lines 


tan 0 ’ : 
L Ratio Ratio 
oe At max. At min. At zero (1) : (3) (2. : (3) 


(1) (2) 


156449 0-18 0:09 


(3) 


73051 0:36 0-25 

154107 0:27 0-18 

233408 0:47 0-32 

227628 0-38 0-28 

EY118K 0:47 Or52 

104-inch 263403 0°47 0:28 
84-inch 137006 0-58 0-36 
230863 0-47 0-28 

174351 0:47 0-32 

143352 0:58 0-36 


Means...... 0°77 0:50 


Suppose that the surface upon which radial lines are imaged (the sagittal 
surface) is generated by the rotation about the x-axis of the curve 
NenOVsE DV eee M8 ait eae (1) 
The secondary maximum in the (R, @) curve will then occur where the above curve 
crosses the focal plane, that is at 


CHO ee ene Wee cea: (2) 
and the secondary minimum where x is a maximum, that is at 
¥1=(«/2B)*. ae) 


The “‘edge”’ of the field cannot be determined easily. However, there-is a fairly 
close limit to the ‘“‘edge”’ since the curve of x falls away from the focal plane 
extremely rapidly after crossing it. It may be assumed that worse definition is 
tolerated at the “edge” of the field than at the secondary minimum in the (R, 6) 
curve. ‘hat is to say, a greater value of x may be tolerated at the “edge” of the 
field than at y,, and the limit of the field may be determined by making » at the 
“edge”? —7 times as great as at y,, where 7 is greater than unity. ‘Thus 


na? /4B = ay," — Bys*, 
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1+7)?}]# ‘ol 

whence p= pee | : | 
H| | 


Thus the ratio yy/yz, which corresponds with the ratio (1) : (3) of table 1, is given byl 


2 4 | iY 
Dd Reser 


and the ratio y,/y3, which corresponds with the ratio (2): (3) of table 1, is given b ; 


nbo= La] 


| 
The first ratio is equal to 0-75 and the second to 0-53 when 7 is equal to 6, which is}} 
probably not an unreasonable value. | | 
The ratio ,/y, is independent of any assumption as to the value of y, and if 
equal to \/2. The experimental value, as determined by the ratio of the 7 values}, 
for the secondary maximum and the minimum, is 1-56, i.e. within 10% of the 
calculated value. Considering that no special efforts were taken to ascertain the} 
precise positions of these features, the agreement is close enough to give some] 
support for the assumption that for these lenses the sagittal field curvature is 01 
the form described by equation (1). It may be observed that the relative positions}; 
of the minimum, secondary maximum, and zero in the (R, #) curves are independent 
of the coefficients « and 8. Equation (1) represents the third and fifth orde 
field-curvature for an infinitely narrow pencil of rays. The positions of the} 
secondary maximum and minimum do not change when the aperture is varied. 
however, and it may be inferred that the above discussion is equally valid fos 
lenses of high aperture. . ' 
In the published field-curves for a number of lenses (for example, Hay and 
von Rohr, 1932) the tangential field-curves display the characteristics of equation 
(1) as often as do the sagittal field-curves, and yet the (R, 6) curves for tangential 
lines rarely exhibit the secondary maximum. ‘The reason for this may be that thd 
aberrations of aperture, which do not apparently affect the validity of the about! 
discussion for radial lines, are more serious in their effect on the resolving powet 
for tangential lines. Such defects as coma and transverse chromatic aberrations] 
as well as the effects of vignetting, might be expected to exert a greater influence 0 : 
the resolution of tangential lines than of radial lines. | 
It has already been mentioned that the drop in resolving power away from the 
axis cannot be explained away as a result of the obliquity of the light. Washe { } 
(1945) has recently shown that even if a lens is completely free from aberrationdf| 
the resolving power will be less for points off-axis than on-axis, owing to the smaller 
angular aperture of the cone of light reaching the image. The effect is, however,}| 
too small to be of any importance in real objectives of the type tested, and in any, 
case cannot be put forward to account for the present results since reducing the} 
aperture improves the definition. 


Variation of illumination across field } 


Typical curves of log relative illumination plotted against angular separation} 
from the axis are shown in figure 14. In most cases the measurements were made} 
up to angles much greater than the lenses were intended to cover. The decrease 
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in illumination from centre to corners of the film size for which the lenses were 
designed was in most cases fairly small. The resolving power of most photo- 
graphic materials increases rapidly as the exposure is increased from zero, and 
reaches a maximum value when the mean density on the film is about 1-0. At 
greater exposures the resolving power falls slowly. If the exposure at the centre 
of the field is such as to give a density of 1-0, at other points in the field the film will 
have a lower effective resolving power if no compensation is introduced. In 
practice the variation of film-resolving power may be reduced to a minimum 
by arranging that the optimum exposure is reached at some intermediate zone 
of the film. 
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Figure 14. Decrease of intensity towards corners of field. 
Broken curve represents [=I cos* @. 


From a curve of film-resolving power against log exposure, the effective film- 
resolving power at the corners of the field was calculated for the standard R.A.F. 
film sizes, and it was found that in only a few lenses was the vignetting serious 
enough to introduce any appreciable change in the effective film-resolving power. 
Moreover, since at the corners of the field the definition is determined by the 
lens rather than the film, the resolving power of the lens-film combination is very 
insensitive to changes in film characteristics. 

Although vignetting tends to reduce the effective film-resolving power, it may 
have a beneficial use in removing some of the aberrant light from outer zones of the 
lens at the higher obliquities. Examples of this effect were encountered in studying 
the results of the present work. 
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§5. ANALYSIS OF RESULTS 


| 
An exact but simple comparison between the performances of different lenses | 
is almost impossible. Not only is the problem complicated by differences in focal | | 

length, aperture, angle of view, and the differences in resolving power for radial i| 
and tangential lines, but the variation of resolving power with angle from the axis. || 
is different for different lenses. If, however, we are prepared to forego detailed |} | 
comparison, and to take a more general view, the problem is not so difficult. The i| 
interpreter of an aerial reconnaissance photograph requires to be able to recognize 
details of objects on the ground. In practice this means, at present standards of | 
performance, that the finer the detail resolved the better the photograph. The | 
value of a photograph is also dependent upon the number of details visible in it, | 
that is upon its area. If, therefore, a system of calculation can be devised which | 
will give a figure for the total number of resolved details in an aerial photograph, | 
that figure may be used as an indication of the performance of the camera in| 
question. Some considerable discussion has taken place on this point and the 
best system of calculation is not yet agreed upon. We shall in this account keep | 
to the original method. 
If the resolving powers for radial and tangential lines are Rp and Ry respectively, 
we can imagine that this represents the resolution of a rectangular unit of area 


1/RpRy. Consider the plane of the film to be divided up into a series of concentric | 
zones by means of circles of radii 79, 79/2, 791/3, etc., so that each zone has area 
mrg2. The number of resolvable elements in a zone is evidently 77,2Rp Rp, | 
if Rp and Ry are the mean values of the resolving powers over the zone. | 
The values of 7) were chosen somewhat arbitrarily to be 10 mm. for a lens of 
5 inches focal length and proportionately for other focal lengths, so that corres- 
ponding zones for all lenses fall at the same angular separation from the axis. It 
follows that the number of resolvable units in any zone is proportional to (Rf)?, 
where R is the geometric mean of Rg and Ry. The figure is thus proportional to. 
the square of the mean angular resolving power. 

The calculation of the averages was facilitated by re-plotting the (R, #) curves 
as (R,) curves, where n is the number of zones enclosed by an angle @ to the axis. | 
and is equal to 161-3 tan? 9. | 

Having thus arrived at a method of estimating the performance of a lens, on a_ 
given film, we are now in a position to examine whether there is any general rule | 
connecting the performance of a lens with its aperture, focal length, angle of view 4} 
and type of film with which it is used. 

There were indications of such a general rule before the full data were obtained. | 
For instance Hay and von Rohr (1932) give figures for spherical aberration, |) 
coma and field curvatures for a large number of lens specifications, taken from | 
patents. ‘These are all for a standard focal length of 100 mm. The curves given |} 
for spherical aberration extend to a value of aperture which is presumably the full | 
aperture at which the lens is intended to work, and the curves for curvature of | 
field extend to an angle which is presumably the angle covered by the lens. When | 
these values are plotted one against the other and a mean curve is drawn through 
them it is found that there is an approximately constant ratio between the F’. No. 
and the square of the angle of view. This is a natural consequence of the difficulty | 
of reducing aberrations to tolerable values at high apertures and wide angles. On 
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the assumption, therefore, that all these lenses were intended to give approxi- 
mately the same level of performance, it may be concluded that the relation between 
the average number of resolvable units per zone, N, the F. No., the number of 
zones, n,,, which the lens covers before the resolving power falls to zero, and the 
focal length, f, will be approximately of the form 


F.. No. — 
eS 


Ny 


since 
62, =tan*6, =n,,/161-3. 

In an attempt to determine the nature of ¢, a three-dimensional model was con- 
structed, in which the average number of resolvable units per zone (obtained on 
Super-XX film) was plotted against (Ff. No./n,,) and f. This model confirmed that 
(Pf. No./n,,) was an appropriate argument in the functional relation, and further- 
more, showed that the #’. No. need not be restricted to the maximum aperture. 
The same relation held, approximately, for all apertures. To test whether this 
was also true if m was varied, the average number of resolvable units per zone 
up to the mth zone was plotted against F. No./n for different apertures. With one 
or two exceptions, a single curve was obtained for each lens ; in the exceptional 
cases, different values of N were obtained for the same value of F. No./n, according 
to whether the F. No. were high or low. In general, however, the agreement 
was good enough to warrant further consideration, and on re-plotting the curves on 
double-logarithmic paper, the points were found to approximate to a straight line 
foreachlens. At this stage the data for the other films were included and a number 
of graphs, such as that shown in figure 15, were obtained. 

Assuming that the points for different apertures lie in the same straight line, 
the average number of resolvable units per zone (N) may be expressed in the form 


= BeNoAe 
en x. 


n 
where A is a function of the focal length and the film properties and 4 is the slope 
of the straight lines and varies, though slightly, from lens to lens and film to film. 
The value of 6 is approximately 0-6. If, on the other hand, it is necessary to 
account for the fact that in some cases the points for different apertures lie on 
different straight lines, equation (4) must be written 


N=A a re (5) 


For the sake of simplicity, we shall adopt the form in equation (4), for we can thus 
set F’. No./n=1 and find a single value of A for each lens-film combination. This 
value of A was found to be directly proportional to focal-length, in the following 


manner :— 


for “ Super-XX ” film A— 3kf, 
for “‘ Panatomic-X”’ film A= 4Rf, 
for ‘‘ Microfile ” film Apa lok: 


The factor multiplying the focal length contains a constant k the value of which 
depends upon the units in which the resolving power is measured and on the size 
of the zones into which the film is divided. 
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The following considerations may help to explain why N is proportional to fi | 
Suppose we have two series of lenses of the same F’. No. but of different focal 


lengths, such that in one series there are no aberrations whatever, and the other 


different focal lengths, and with considerable amounts of aberration. For th 


lenses without aberrations the resolving power in lines/mm. will be the same for alll). 


1] 


focal lengths, because the image size is independent of the focal length. The 
angular resolving power will therefore be proportional to the focal length an i 


B-INCH F/56 I4-INCH_ F/4-5 
264095 227628 


eke) 


Z| 


2O1NCH F/63 : 3eINCH F/63 
5337/LP ' 275668 


100 


F NUMBER 
n 


Figure 15. Variation of N with F. No./n. 
(a) Microfilm film, (6) Panatomic-X film, (c) Super-XX film. 


N will be proportional to the square of the focal length. When the lenses have 


considerable amounts of aberration, the size of the image is proportional to the} 


focal length, so that the angular resolving-power, and WN also, will be independent 
of focal length. In fact, even in this case, N will increase a little with focal length 
because the diffusion of light in the emulsion layer will be of less importance 


when the optical image is large, and also because the threshold of resolution is | 
lower for large images than for small. Now the path-length errors are propor-|| 
tional to the focal length, and, therefore, however wellor ill the lens may be designed, | 


1M 
series is of lenses all of the same design, but made to different scales, so as to have} 
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there is some short focal length at which the lens is nearly perfect and some long 
focal length at which the correction is bad. In other words, for short focal 
lengths N will vary as f2 and for long focal lengths will be constant, or hearly so. 
For moderate focal lengths, therefore, it is not unreasonable to find N varying 
approximately as f. What is perhaps more surprising is that the proportionality 
between N and f is not seriously disturbed by changes in aperture. Figure 16, 
which is illustrative only, will perhaps make the argument clear. 

The general level of the resolving power varies according to the film used for 
taking the photographs. ‘Two quantities have arisen during the theoretical 
discussions to which these differences of level might be attributed. They are the 
granularity and the diffusion of light in the emulsion layer. The dimensions of 
granularity are [Density: Length] and the dimension appropriate to the diffusion 
of light is [Length], since it may be conveniently measured by the distance in 
which the light flux in the emulsion layer coming from a long illuminated narrow 
strip drops by the ratioe. Either quantity, therefore, will satisfy the dimensional 


UMIT FOR LONG if 


2| 


VA 
| APPROXIMATE LIMITS OF —————>| 


VALIDITY 


t 
Figure 16. To show that for short focal iengths Ne f?, for long focal lengths N is constant, 
and that intermediately Nx f, approximately. 


requirements of the formula. _ But in fact the diffusion of light is not likely to be 
of any importance. It will be seen from figure 2 that log P, does not drop much 
below 0 for values of 1/d less than 20. ‘Values as high as this are observed only 
onthe axis. Assoonas any appreciable angle is covered, the resolving power drops 
considerably below 20 lines/mm. and the diffusion of light in the emulsion layer is 
an unimportant item in the degradation of the image. 

The granularity values, measured by the method described by Romer and 
Selwyn (1943), at a density of 1-0 are 1-6, 1-2, and 0-3 density-microns for Super- 
XX Aero, Panatomic-X Aero and Microfile films respectively. The reciprocals 
of these figures are in the ratio 3:4: 16, from which it may be concluded that the 


complete formula is 


eee el tN Olas 

pies n ] ‘ 
in which b has the approximate values of 0-55, 0-6 and 0-75 for Super-XX Aero, 
Panatomic-X Aero and Microfile films respectively. Substituting these values 
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is found to be approximately 5-5 x 104. 
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and in table 3 the frequency with which these differences occur is shown. 


‘Table 2. Mean percentage differences of experimental values of N from those| 


Film 
Lens 


5-inch F/4 143352 
8-inch F/4 230863 
8-inch F/2:9 174351 
8-inch F/5-6 264095 
81-inch F/4 137006 
104-inch F/6 263403 
14-inch F/4:5 227628 
14-inch F/6-3 EY118K 
14-inch F/5-6 72754 
20-inch F/5-6 73051 
20-inch F/5-6 277490 
20-inch F/6-3 154107 
20-inch F/6-3 174856 
20-inch F/6:3 5337/LP 
50 cm. F/5 2451477 
24-inch F/6 7887 
24-inch F/6 EA504 
30-inch F/6-3 156449 
36-inch F/6-3 275667 
49-inch F/7 402490 


given by equation (6) 


Super-XX Panatomic-X 
(b=0°55) (6=0-6) 
—21 —10 
—14 —13 
—12 — 8 
—15 — 6 
— 7 +21 
+ 4 + 3 
— 24 — 24 
— 5 — 8 
+25 +27 
+16 +23 
—19 a 
— 3 — 8 

+10 
+10 — 1 
—10 
— 3 
ee 23 
a 34 Sele 
— 9 — 1 
—23 —28 


Microfile 
(6==0°75)r 


Wall al Fs aes a teal 
Bmw 
NEN SUN OYwW 


aul) 
+15 


Table 3. Distribution of differences in N. 
Percentage No. of cases | 
difference % j 

(Super-XX) |(Panatomic-X)} (Microfile) 
10°% or less 10 8 ‘i 
11%-20% 5 2 6 
20° or more 5 5 6 


] 


i) 


i 


| 


in (6) and comparing the calculated value of N with the experimental values, Kil), 
In table 2 the departure of the perfor 


mance of individual lenses from that of the “average” lens as given by (6) is shown 


1 2 


The standard deviation for repeated measurements of the resolving power}! 


of the same image on the same negative is of the order of 10°, but in view of the} 


i] 


large number of readings made for a given lens, at different angles to the axis,| 
and since all the readings were made by the same observer, the averaged values of i 
resolving power are probably reliable to within much less than 10%. (It should |} 
be noted that since N is proportional to the square of the angular resolving power, 
the differences in terms of angular resolving power are approximately half of | 
the values listed in table 2.) The differences shown in table 2 may therefore 


be regarded as real differences, and due to differences in performance between 


different lens types. 
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A single value of 6 may be adopted without serious change in the reliability of 
the empirical relation. ‘Thus, eliminating the parameters introduced originally, 


= Ff’. No. |" 
N=5:5x10!x Al 7 =] resolvable units per zone 
F. No. |° 
=2:6 x 108 x abe since n=161-3 tan? 4 


ye 10? cline NGa| ee resolvable units per sq. mm. 
= * PR tan? 6 on the film 


_ 207 F. No. 
~ fG | tan? 6 


Ihe mean resolving power on the film is thus 


a NOM co ee 
Fame lines per mm. 


ind the mean angular aes power, in radians“, is 


207f 7]? [F.No.7°3 
25-4 =| [ars | ne a i er (7) 


ince f is measured throughout in inches. 


§6. DISCUSSION OF EMPIRICAL RELATIONSHIP 


The formula which has been arrived at to represent a general view of the per- 
ormance of aircraft camera lenses has no fundamental or theoretical basis to 
nsure that it is generally valid. ‘This has been made plain during the latter part 
f this investigation, as a result of measurements on new designs. These have 
hown some improvement in performance, which, while it does not so far make 
he above formulae seriously in error, undoubtedly shows that with continued 
mprovement a new formula, approximating more closely to the theoretical limit 
f perfection, will be necessary. For the time being, however, the above formulae 
re a satisfactory statement of the level of performance to be expected from a lens 
f given dimensional properties, provided that these dimensional properties are 
ot outside the range of those of the lenses tested. 

To illustrate this last point it is convenient to use again the argument of §5. 
t will be obvious that the rule Nec f cannot be obeyed for very short focal lengths 
or which Nef? ,so that the empirical formula predicts higher resolving powers 
han are theoretically possible. And it will also be obvious that Nec f cannot be 
rue beyond focal lengths at which the lens may be regarded as badly corrected. 

The correction of aberrations is better in the centre of the field than towards its 
uter parts, with the result that the diffusion of light in the emulsion layer is 
elatively more serious in the centre of the field than elsewhere. ‘Thus the 
nportance of the diffusion of light in the film is dependent upon the angle of view. 
n the formula the influence of diffusion does not appear explicitly, but may be 
iwvolved implicitly in the term G, since there is an approximate connection between 
* and the diffusion. But the influence of granularity also varies over the field of 
ne lens, as a result of its relatively greater obscuring effect for fine detail than for 
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large. ‘To represent the influence of granularity by a single term unconnected) 
with the angle of view is, therefore, theoretically incorrect. It is more tha | 
probable that the different values of 6 found for the three materials arise because off 
this effect. Unfortunately the data are quite inadequate for investigating this 
question more fully, and any attempt at more precise formulation of the experi-+ 
mental data than is represented by equation (5) would be somewhat artificial. 

There are a few minor limitations to the generality of the results. The average; 
resolving power was computed for circular films, since the analysis would have; 
been rather complex if rectangular films had been considered. In addition, 
in the collimator method of examining lenses, only a small fraction of the total 
field of the lens is illuminated, and the effect of light scattered by the lens is les 


than under practical conditions. Another departure from strictly practica 
conditions was that the brightness of the test-object was increased for extra-axiall 
exposures to compensate for the decrease in illumination due to vignetting, etc. 
In spite of the foregoing discussion, it is somewhat surprising, in view of the 
known differences in construction of the lenses under examination, that their 
performances may be defined, as well as they are, solely by the dimensional) 
properties of focal length, aperture and angle of view. ‘The empirical formula) 
implies the existence of a standard of performance, and this can only mean 
that the designer endeavours to reduce the aberrations to a level which 
experience has taught him will lead to a satisfactory result. Another designer 
may aim for a different type of correction, but the extent of the correction upo 
which he decides will probably be such as to lead to a photographic performance} 
approximately equal to that accepted by the first designer. This is not to say tha 
the empirical expression is inflexible; the lenses upon which it was based were} 
designed during and after the 1914-18 war and up to the early years of thet 


second war, and all that is claimed is that the formula represents the standardl 
prevailing up to about 1943, when this work was completed. Indeed, since that} 
time advances have been made which have raised the performance of individuallf} 
lenses to such an extent that a repetition of the present research with a wide enoug! 
range of the latest products of the British optical industry would result, if not in e 
change in the form of the relationship, then certainly in a change in the constants.|} 
For the time being, however, the formula is useful, not only for its original purpose, 
but also as a standard of reference by which the merit of a given lens for air photo- 
graphy may be estimated by comparing its performance with that calculated from! 
the formula, and for estimating the effect of moderate changes in the focal length, 
aperture and angle of view of a lens of given basic design. | 
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ABSTRACT. A critical examination is made of the Birge-Sponer method of deter- 
mining dissociation energies of diatomic molecules. Dissociation energies of certain 
classes of molecules as determined by linear extrapolation are compared with values 
determined independently. For the ground states of ordinary molecules the linear extra- 
polation tends to give too high results, often by around 20%. For molecules in which 
ionic binding is important the extrapolation often comes too low, sometimes much too low. 
For the molecules in which one atom has a 4S, ground state the linear extrapolation is 
‘much too high. The method appears to work satisfactorily for ionized molecules. Excited 
states of molecules cannot, in general, be relied on to give good extrapolations. These 
empirical observations are discussed in terms of the potential-energy curves and the 
binding forces in the molecules at large inter-nuclear distances. 


§1. INTRODUCTION 


by extrapolation of the vibrational energy levels to their convergence limit 

was first put forward by Birge and Sponer (1926). Since then some other 

methods have been used as well to ascertain the dissociation energies of certain 

‘molecules, but for the large majority of diatomic molecules we are still dependent 
on the Birge-Sponer extrapolation. 

When data are available for a large number of vibrational levels, the graphical 


| HE method of determining the dissociation energies of diatomic molecules 
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extrapolation, made by plotting the vibrational energy intervals AG(v+ 3) against) 
the vibrational quantum number 2, or sometimes against the vibrational energy; 
term G(v), may be expected to lead to a fairly reliable result. In the majority) 
of cases, however, only a few levels are known and the interval between these: 
decreases approximately linearly, and it is then customary to make a linear’ 
extrapolation. It has been generally realized that there is a tendency for this} 
linear extrapolation to come too high. Birge (1929) discussed the form of the} 
curve of AG(v+4) against v for a few molecules and showed that its shape and) 
curvature changed rather unexpectedly in some cases ; the sudden inflections} 
which he believed he had established for the excited states of O, and some of the} 
halogens, have not, however, received much support from later, more precise, | 
spectroscopic observations. Rydberg (1931) suggested that an alternative form) 
of extrapolation, in which [AG(v+4)]? is plotted against v, might be better in} 
some cases ; and since this method usually gives a lower value of the dissociation. 
energy than that obtained by linear extrapolation, it has sometimes been used.) 
Nevertheless, the linear extrapolation is still in most general use and forms the} 
basis for most of the values of dissociation energies now in the literature. | 

For a fair number of molecules accurate values of the dissociation energy have) 
been obtained by other methods, such as by observation of band structure up to 
the convergence limit at which continuous absorption begins (O, and the halogens), 
of atomic fluorescence (halides of alkali-metals, Tl, etc.), of predissociation 
(HgH, N,, CO, SO, etc.), or by thermochemical methods (alkali-metal halides, | 
etc.) or by combinations of these methods (HCI, HBr, HI, NO). ‘These methods} 
have been reviewed in the author’s recent book (Gaydon, 1946) and the values for} 
the dissociation energies have been listed. It is the purpose of this paper to make } 
a systematic comparison between these accurate values of dissociation energies | 
with those obtained by linear extrapolation, with a view to making better use of 
the linearly extrapolated values for the many molecules for which this is the only | 


| 


method available. 

In their original paper Birge and Sponer pointed out that their method of| 
extrapolating was only justified if the vibrational energy levels converged to a| 
limit at a finite value of the vibrational quantum number v. Kratzer (1924) 
showed that if the force between the nuclei at very large internuclear distance, 7, | 
is expressed as a power series in (1/r), then, if the dominant term at large 7 is 1/7? 
or 1/r’, there will be an infinite number of vibrational energy levels, so that the} 
curve of AG(v+ 3) against v would approach the v axis asymptotically. If, 
however, the dominant term in the expression for the force is 1/r* or higher | 
power of (1/r), then the maximum value of v will be finite. For non-polar molecules 
the force at large internuclear distance varies as a high power of (1/r) and the| 
extrapolation is therefore permissible. For polar molecules dissociating to| 
ions, however, the force must vary as 1/r? at large distances and the extrapolation 
should not be made by the method of Birge and Sponer. It is clear that the rate 
of convergence of the vibrational levels near the limit will depend on the form of 
the law of force at large distances and we may, therefore, expect the curve of 
AG(w + 4) against v to show a different trend for classes of molecules with different 
forms of binding. It therefore appears best to treat separately certain general 
classes of molecules. 
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§2, EXTRAPOLATION FOR THE GROUND STATES OF NEUTRAL 
MOLECULES WITH NON-IONIC BONDING 

First, we shall compare the accurately determined dissociation energy D, and 
the linearly extrapolated value Dj, for the ground electronic states of those 
neutral molecules for which ionic binding forces are not unduly important. 
The values of the dissociation energies are taken from the author’s book, where 
full references may be found. For N,, NO, CO and CN the values have been 
discussed by Gaydon and Penney (1945) and differ from those in the older literature 
(e.g. Herzberg, 1939). For SO and Sy, values higher than those given by Herzberg 
have been adopted. 

It is possible to make the linear extrapolation either graphically from the first 
part of the curve of AG(v+4) against v or analytically. The latter method is 
more usual and has been adopted here. 

If the energies of the vibrational levels are expressed as a power series 


G(v)=a,(u-+4)—xw (Ot) tec., «nave es (i) 


where wv takes integral values 0, 1, 2, ...., it can be shown that, neglecting higher 
terms, the energy interval from the lowest level, v=0, to the limit, which is the 
required dissociation energy , 1s 

wee 


Snore salen fy ree pee) Wilt eis (2) 


in > 4x0, 2 


It is often given as D,,,, = w?/4x,w, where wy =w,—x,w, ; for all practical purposes 
this is satisfactory, but it may be pointed out that it is not strictly correct as it 
differs from (2) by 4x,«,. 

The values of D,,, have been calculated with the aid of equation (2) from 
collected values of w, and x,w, given by Herzberg (1939) ; for CIF the value is 
from a more recent paper by Wahrhaftig (1942). Table 1 shows the values of 
D>; Din» and the ratio Dy/D 


lin? lin* 


Table 1. Values of Do, Dy, (in ev.) and Dy/D,,, for the ground states. of 
molecules with non-ionic bonding 


Molecule Dy Din = Do) Onn Molecule De Die De) Dis 
Hy, 4-478 4°55 0-99 Oz 5-084 6°35 0-81 
KH 1-86 2-2 ue 0-92 sO 5-148 6:3 0-82 
AgH Ia) 27, 0°85 Se 4-4 57) 0-77 
OH 4-4 553 0-83 Ses 2°8 4-5 =» 10-62 

TeO 27725 D5 0-48 
HCl 4-431 S15 0-86 
HBr | 3-60 4-78 0-75 Nz 9-764 11:8 0-83 
HI 27.9) 4-0 0-69 P, 5-033 6°66 0-76 
As» 3-94 Bo O77 
Li, dst? 1-46 0:77 Bi, 1-72 227. 0-76 
Naz 0-77 1-07 Oey 
K, 0:51 0:74 0-69 CIF 2-616 1:93 d-S5 
(0 2-476 2°44 1-02 
CO plea (1 0-99 Tet 2°153 3:08 0-70 
NO 6°49 7°6 0-85 Br, 1-971 3-02 0-65 
CN 7°6 9:9 OF IBr 1-818 2°85 0-64 


528 A. G. Gaydon 


The table includes practically all the non-ionic molecules for which adequate | 
data for D, and x,w, are available. The exceptions are certain molecules con- || 
taining elements of Group II of the Periodic Table, which have *S, ground states |} 
and which, for reasons given presently, are best considered separately. It will be |} 
noted that OH, HCl, HBr and HI are included in table 1 ; ionic forces probably |} 
contribute to a certain extent to the bonding in these molecules, but since the 
electron affinity of the non-metal is in each case small (around 3 to 4 ev.) compared 
with the ionization potential of hydrogen (13-5 ev.), these molecules will dissociate 
to neutral atoms and the force between the atoms at large distances is unlikely to be 
appreciably affected by ionic forces, 

It will be seen from table i that there is a general tendency for the true dis- 
sociation energy to lie below the linearly extrapolated value. The average 
value of Dy/D,;, is 0°79, and for thirteen of the 28 molecules the values lie between 
0:75 and 0:85, with six above 0:85 and nine below 0:75. ‘Thus we see that for 
this class of molecules the most probable value of the dissociation energy is about 
20% below that given by the linear extrapolation. Clearly the scatter in the 
values of D,/D,;, is rather great, but apart from the general tendency for D,,, 
to be too high it is possible to distinguish at least one systematic effect. Do/Dyy 
tends to be highest for molecules composed of light elements and lowest when 
heavy elements are involved. This effect is well shown both by the group 
containing the halogen molecules and by that consisting of molecules built from 
Group VI elements. It is not easy to express this effect in quantitative form, 
but it would seem that one might do worse than to adopt a dissociation energy 
10% below the linearly extrapolated value for hydrides and molecules composed 
of light elements, while for molecules containing a very heavy element a value 
only 70% of that given by linear extrapolation seems reasonable. 

It may be pointed out that if the experimental data are not of high accuracy 
there may be a tendency for the next term, y,w,(v +4)’, in the series expansion 
for G(v) (equation (1)), to be masked by the experimental error. Since this term | 
is usually negative, the result will be to obtain a slightly higher value of the 
apparent x,w, than the true one, thus bringing D,,,, down slightly. 


§3. MOLECULES CONTAINING AN ATOM WITH 1S, GROUND STATE 


The form of the graphical Birge-Sponer extrapolation or the accuracy of the 
linear extrapolation must clearly be related to the potential-energy curve for the 
molecular state, especially its shape at fairly large internuclear distances. Now 
an atom with a closed shell or closed sub-shell of electrons, i.e. an atom in a 1S, 
state, has no resultant electric or magnetic moment and is, therefore, likely to exert 
a smaller force on another atom at a distance than is an atom with a resultant 
moment. This is borne out by the fact that the rare gases do not, with the 
exception of some excited states for Hes, form molecules at all. For the elements 
of Group II of the Periodic Table, however, although their ground states are 
‘So, many stable diatomic molecules are known. _ In these cases the binding forces 
are believed to be mainly due to excited states of the atom, although as a result _ 
of interaction between the various potential-energy curves and the rule that 
curves of the same species may not cross, the ground states of these diatomic 
molecules do probably correlate with ground atomic states. Clearly, the form 
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of the extrapolation is likely to be affected by the fact that the binding is largely 
due to excited states, and that at large internuclear distances the force between 
the atoms is likely to be abnormally small. It thus seems desirable to treat 
molecules of this type separately. 

Although a large number of diatomic molecules formed from elements of 
Group II are known, for relatively few of them is it possible to compare D, and Dy, ; 
apart from the scarcity of accurate values of the dissociation energy, for many 
of these molecules, especially the halides, the band systems are not strongly 
degraded, so that the values of x,w, derived from the analyses of band heads are 
far from accurate. ‘The data are most complete for the hydrides; for many 
of these predissociation has been observed, which sets at any rate an upper limit 
for Dy. The values are given in table 2. 


Table 2. Do, Dy, (in ev.) and D,/D,;,, for hydrides of Group II elements 


Molecule Ds Dyin DolDun 
ZnH 0:84 1°35 0:62 
CdH 0-678 1-29 0-53 
HgH 0-372 0:64 0:58 
CaH < 1-70 2°6 < 0°65 
SrH = 1:66 2°6 = 0:64 
BaH z 1°82 2°6 Z 0:69 


Thus for these hydrides the true value of the dissociation energy is much less 
than the linearly extrapolated value, often being only about 60% of that value. 
The unusually rapid rate of convergence of the vibrational levels for these mole- 
cules is shown by the high values of x,w, and the strong negative curvature of 
the graphs of AG(vw++4) against v in many cases. For molecules other than 
hydrides the evidence is scanty. For CaCl predissociation shows that D)<2:76, 
while a poor linear extrapolation gives 5-1 ev., so that Do/D,;,<0°54. For HgCl 
recent measurements (Wieland, 1941) give a linear extrapolation to 1-6 ev., and 
a good graphical extrapolation to 1:13, the ratio of which is 0°70. While the 
data are clearly inadequate for sweeping generalizations, we may reasonably 
suspect that molecular states correlated to an atom in a 1S) state tend to give 
an abnormally high linear Birge-Sponer extrapolation. ‘The alkaline-earth 
metals mostly have low ionization potentials, and many of their compounds, 
especially their halides, may show partially ionic bonding which may confuse the 
relation between their true dissociation energies and the linearly extrapolated 
values. 


§4. MOLECULES WITH PARTIALLY IONIC BONDING 


For molecules composed of a metal with low ionization potential and a non- 
metal with high electron affinity the bonding force tends to be largely of an ionic 
nature. For all molecules except CsF the ionization potential of the metal exceeds 
the electron affinity of the non-metal, and the ground states are believed to 
dissociate to neutral atoms. Nevertheless, since the ionic forces in some of these 
cases are believed to be relatively important, it is probable that the potential 
energy curve will assume a different form from that for molecules in which ionic 
bonding is unimportant. ‘hus, it seems likely that the relation between the 
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true value of the dissociation energy and the linearly extrapolated value may | 
be abnormal. In making the comparison we are again handicapped by lack of 
data. The most ionic molecules are the halides of the alkali metals. Good | 
values of the dissociation energy are available for most of them, thermochemical || 
values agreeing closely with those obtained by observation of atomic fluorescence. | | 
These molecules do not, unfortunately, readily give discrete band spectra from | 
which values of w, and x,w, can be obtained ; for some of them, however, Beutler | 
and Levi (1934) have succeeded in observing some rather diffuse bands from which 
approximate values of w, and x,w, have been derived. The only other molecules | 
of this type for which data are satisfactory are TICI, TlBr, and Gal ; the measure- || 
ments are most precise in the case of TICI, for which the dissociation energy has : 
been derived both thermochemically and by atomic fluorescence, while the | 
values of w, and x,w, have been obtained by Howell and Coulson (1938) from | 
origin measurements. | 

The values of Do, Din, Do/Dyn» and also the difference between the metallic | 
ionization potential (/) and the electron affinity (Z.A.) of the non-metal are shown 
intable3. — It will be seen that for most of the molecules listed the true dissociation | 


Table3. Ground states of molecules with predominantly ionic bonding 


Molecule Ds; Din DolDun I-E.A. 
Lil 3-6 4-2 0-86 1-6 
NaCl 4-24 4-4 0-96 1-4 
NaBr SSIS 2-65 1-43 1°5 
Nal raat 3-4 0-91 Dit 
HCH 4-40 DEOT 1-63 0-6 
KBr BE9S USD 1:78 0:7 
KI Ses 2-0) 1:65 iheil 
TIC! 3-80 OE pe tbersy D8 
T1Br B 2-9 1:10 25 
Gal 2°85 2-9 0-98 2°8 


* From Howell and Coulson (1941). 


energy is much higher than that obtained by linear extrapolation ; this is particularly 
true for the halides of potassium, for which the values of x,w, are likely to be more 
accurate than for the halides of Na and Li. There is some doubt about the values | 
of Dy for GaCl, GaBr, InCl, InBr and InI, but here again it is quite likely that the 
true dissociation energy is much higher than the linearly extrapolated value. | 
For HF, although the ionization potential of hydrogen exceeds the electron 
affinity of fluorine by over 9 ev., itis believed that ionic bonding is fairly important; 
if we accept Herzberg’s (1939) value for Dy we obtain D,/D,,,=1-07; a rather 
lower value of Dy appears probable to the author, but it seems fairly certain that | 
the ratio is not less than 1. 

Thus we see that values for the dissociation energies of molecules with 
partially ionic bonding which have been derived from a linear Birge-Sponer 
extrapolation should be treated with caution. In some cases such values may be 
too low, instead of, as is usually thought to be the case, on the high side. It is not 
always easy to say definitely which molecules are subject to appreciable ionic 
binding forces; the rough rule, suggested by Mulliken, that ionic bonding is 
likely to be predominant when the metallic-ionization potential exceeds the electron 


The determination of dissociation energies 531 


affinity of the non-metal by less than 3 ev., is helpful, but there are cases when 
such forces are still important even when the ionization potential is a good deal 
higher than the electron affinity, as in HF. 


§5. THE EXTRAPOLATION FOR EXCITED STATES 


So far we have considered the extrapolation only for the ground states of 
molecules. For these, the perturbations of the potential-energy curve by those of 
other states are likely to cause only a general deepening of the potential minimum 
without any violent alteration in its general form. (The special cases where the 
ground state correlates with an atom is a 4S, state, and where the perturbations 
are due to ionic terms, have been discussed already.) For excited electronic 
states, however, mutual interactions are more likely to cause perturbations which 
will modify the shape of the potential curves in an irregular way, and deviations 
from a smooth Birge-Sponer extrapolation of the vibrational energy levels may be 
expected to be more frequent. 

It should be especially noted that interaction from ionic terms is likely to be 
more frequent and more severe for excited states. Thus, for the halides of Al, 
Ga, In, and TI, the dissociation limit to ions (e.g. to Alt + CI, etc.) lies below the 
limit going to a normal halogen atom and a metal atom in its lowest excited state* 
(e.g. to Al(2S)+ Cl(??P)). Many of the excited states of molecules formed from 
one strongly electro-positive and one strongly electro-negative atom are influenced 
by ionic forces, and may even dissociate to ions; in this case the Birge-Sponer 
extrapolation would take quite a different form, and it would be quite incorrect to 
assume a linear convergence of the vibrational levels, of which there would in this 
case be an infinite number. 

There are a fair number of examples of abnormalities for excited states. 
The most obvious group calling for comment is that of the upper states of the 
main systems of the alkali-metal hydrides. For all these the spacing of the 
vibrational energy levels at first increases so that x,w, has a fairly large negative 
value; instead of converging, the levels diverge. The unusual form of the potential 
curves for these states has been attributed by Mulliken (1936) to interaction 
with the potential curve arising from ionic products. Although a linear extra- 
polation is clearly impossible, the curve of AG(v+4) against v appears to be 
smooth, rising to a flat maximum and then falling fairly sharply, so that a graphical 
extrapolation to the dissociation limit is not necessarily ruled out ; a good value can 
be obtained in this way for KH, and the method can also be used for NaH. For. 
some of the excited states of the alkali-metal molecules (K,, Lig, etc.) a similar, 
but less marked, effect is observed. The vibrational energy levels at first show a 
slow rate of convergence which slowly increases, but at high vibrational quantum 
numbers the rate of convergence increases sharply and the graphical extrapolation 
appears to indicate a dissociation energy for these states of only about half that 
given by linear extrapolation. 

The most extreme anomaly is for the excited state of AgH. From the analysis 
by Bengtsson and Olsson (1931) it appears that the vibrational levels at first show 
rapid and increasing convergence, but the curve of AG(wv + >) against v flattens out 


* By lowest excited state is meant the lowest main term other than the ground one, ignoring 
multiplet splitting; for Al, etc., the ground state is *P and the lowest excited state 29. 


abruptly at higher v, and the rate of convergence must be practically zero for a) 
number of levels before increasing again (see Herzberg, 1939; Gaydon, 1946). 

It has been shown (Herzberg and Mundie, 1940) from the form of the limiting | 
curve of dissociation, as obtained from predissociations, that the excited "II states | 
of AIH and BH possess potential maxima. Mulliken (1937) has shown that some |} 
of the potential curves of the excited states of the silver halides, and probably the} 
copper halides, possess quite marked maxima. In these cases even a graphical | 
Birge-Sponer extrapolation would be of no avail, as it would lead to the energy of | 
the potential maximum instead of to the dissociation limit. The author is} 
inclined to think that the potential curves of the excited states of many other | 
halides, including probably those of Al, B, Ga and In, may possess similar 
maxima. | 

Although it is clear that there are many instances in which the extrapolation | 
(especially the linear form) is likely to fail when applied to excited molecular | 
states, there are still many for which it seems to work satisfactorily. As examples | 
may be quoted the upper states of the main systems of O,, P, and the various | 
halogen molecules; it is true that at one time Birge (1929) pointed out certain 
irregularities in the curves of AG(v+4) against v for some of these, but these 
irregularities have been largely smoothed away by later experimental work, and 
were in any case not of a magnitude seriously to affect the accuracy of the extra- | 
polation. Many of the lower excited states of N, and CO give satisfactory extra- | 
polations to the expected limits when the author’s (Gaydon and Penney, 1945; | 
Gaydon, 1946) latest values for the dissociation energy are used. | 

In general it seems that where a reasonable extrapolation from ground state | 
data is available it is to be preferred to one for an excited state, but for molecules 
with atomic (non-ionic) binding, especially homonuclear molecules and molecules 
formed from non-metallic elements, extrapolations from the lower excited states 
are permissible. For metallic halides, and probably hydrides and oxides as well, |}} 
it is unsafe to place reliance on extrapolations for excited states. 


§6. THE EXTRAPOLATION FOR IONIZED MOLECULES 


So far we have considered the extrapolations for neutral molecules which 
dissociate to two neutral atoms. For singly ionized molecules, however, dissoci- 
ation will take place to one neutral atom and onecharged atom. The force between | 
a charged and a neutral atom at large distances is likely to be much greater than | 


Table 4. Ionized molecules 


Molecule Do Dr, Do/Dys, 
H,> 2-649 2°50 1:06 
HClt 4-5 4-0) WD 
O,t 6-48 6-48 1-60 
N,* 8°73 9-2 0-95 
Cor 8-2 9-8 0-84 


that between two neutral atoms. Thus the potential-energy curve at large 
internuclear distances is likely to take a rather different form for ionized molecules, 
and it would seem best to make a separate comparison of the true and linearly 
extrapolated values of the dissociation energy for this type of molecule. The few 
available values are collected in table 4. 
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For O,+ the value of D, is derived from a good graphical extrapolation for the 
excited state; it is probably of quite high accuracy, but does, of course, assume a 
normal trend for the extrapolation. The values for N,+ and CO+ are derived 
from the values of D(N,) and D(CO) given by Gaydon and Penney (1945) com- 
bined with the ionization potentials of N, and CO; for CO the value 14-1 ev. for 
(CO) has been used, but it is just possible that the value is 14-55 ev. 

From the values collected in table 4 there appears to be little evidence for 
any systematic deviation from the linear value. It would therefore seem that 
ionized molecules are rather more suitable for making a linear extrapolation than 
neutral molecules, for which the true value, as we have seen, tends to be around 
20°%% below that of the linear extrapolation. | However, all the molecules 
included in table 4 are relatively light. For BeH*+, MgH+ and CdH? accurate 
values of Dy are not available, but the rate of convergence of the ground state 
vibrational energy levels does increase slightly with v, and the linear extrapolation 
is probably a little high. | Most probably the linear extrapolation is fairly good 
for light ionized molecules and a little high for the heavier ones, although not to the 
same extent as for the corresponding neutral molecules. 


§7. THE RELATION BETWEEN THE FORM OF THE 
EXTRAPOLATION AND THE POTENTIAL-ENERGY CURVE 

Clearly, the rate of convergence of the vibrational energy levels must be 
related to the shape of the potential-energy curve for the particular molecular 
state. Now for the lower vibrational levels the rate of convergence is frequently 
nearly a linear function of the vibrational quantum number, and we have been 
chiefly concerned with the linear extrapolation. One of the simplest functions 
for representing the potential-energy curve of a diatomic molecule is that due to 
Morse, which is of the form 


U=D, (1—ePe-ro)2, 


One of the main virtues of this function is that when inserted in the wave equation 
it leads to a series of vibrational energy levels which do converge linearly. Thus it 
seems desirable to make use of this simple function as a basis for comparing the 
form of the potential curve and departures from linear convergence of levels. 

For normal molecules, in which ionic binding forces are unimportant, the 
lower vibrational levels usually show a very nearly linear rate of convergence, 
put, as we have seen, extrapolation of this convergence leads to a dissociation 
snergy of the order 20% too high. We may take the ground state of O, as 
1 typical example. For this, some 21 vibrational levels are known up to over 
3-4 ey., and over this range the rate of convergence is nearly linear, showing only 
1 slight sign of an increased rate at the highest observed values of v. Figure 1 
shows a curve (dashed line) obtained from a Morse function (which uses the 
yalues of D,, w, andr,, which are derived from the three spectroscopic constants 
Do, », and B,) and a curve (full line) plotted from a modified function proposed by 
Hulburt and Hirschfelder (1941) which uses the five spectroscopic constants 
Do; Mey XH», B, and «,; the five-constant function can be relied upon in this case 
‘© come very close to the true potential curve. It will be seen from figure 1 that 
tt r>r, the true potential curve lies above that given by the Morse function. The 
rue curve tends to approach the dissociation limit more rapidly than the Morse 
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curve, and this more rapid approach to the limit is associated with a less rapid} 
rate of convergence of the lower vibrational levels and an increasing rate of con-}| 


vergence of the higher ones. - a. 
For molecules dissociating to atoms of which one is in a 1S) state, the vibrationa 
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Figure 1. Potential-energy curves for the ground state Oy. 


levels converge with unusual rapidity; even for the lower levels there is usually a 
marked departure from a linear rate of convergence, and the linear extrapolation, 
for which there is really little justification, comes much too high. Clearly, ir 
these cases the potential curves must depart considerably from the Morse form | 
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Figure 2. Potential curves for the ground state of CdH. 


Owing to this big difference it is not easy to plot the potential curve, even whe 
experimental data are adequate, as for many of the alkaline-earth hydrides. 
Figure 2 shows the potential curve for the ground state of CdH drawn to a Morsey) 
function (dashed line) and also a curve drawn to a five-constant Hulburt- 
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Hirschfelder function (full line, with upper part dotted). This latter curve 
orobably approaches the true curve quite well around the minimum, but may be 
subject to appreciable error as r increases, and the potential maximum to which it 
eads should not be taken seriously. However, it serves to show that the strong 
leparture from a linear rate of convergence of the vibrational levels is again 
issociated with a much more rapid rise of the potential curve towards the dissoci- 
ition limit than that given by the Morse function. 

The differences between the actual potential curve and the Morse form so far 
vonsidered have been in the region where 7 is greater than 7, but is still not very 
arge (~2r,); this is probably the region of most importance for the extrapolation, 
ind here the tendency is for the rate of convergence to be initially less than that 
siven by the Morse function, but to increase and exceed it at high wv, and for the 
rue potential-energy curve to lie above that of the Morse function. At larger 
nternuclear distances there may be a tendency for the reverse effect. The 
Morse function corresponds to an exponential falling off of the potential energy at 
very large r. ‘The actual potential energy at very large distances is, however, 
ikely, on account of interaction of two atoms, to fall off much less rapidly 
han that given by the exponential law. The weak van der Waals forces between 
itoms are known (see, for example, Margenau and Watson, 1936) to fall off as 
/r’, corresponding to a potential energy varying as 1/r6. There may be other 
orces, which vary as a still lower power of (1/r): thus the resonance interaction 
yetween like atoms in different states gives a force proportional to 1/r*; polarization 
orces also are likely to be important in many cases, and these will again give a 
yotential energy which, at large distances, falls off less rapidly than it would if it 
beyed the exponential law. The force between two atoms in P states may 
ary as l/r’. ‘Thus the true potential curve should cross the Morse curve some- 
vhere at large internuclear distances and then approach the dissociation limit less 
apidly. It would not, therefore, be surprising if towards the limit the rate of 
onvergence of the vibrational levels were to decrease slightly. 

There are not many examples for which we have adequate data on the vibra- 
ional levels very near to the dissociation limit, but there is some evidence 
0 support this conclusion. For the upper state of O, and for the ground state of 
1, Birge (1929) assumed an inflection in the curve of AG(v+ 4) against v at high v, 
ut later experimental work has not provided much support for any marked 
hange, although a slight trend towards a slower rate of convergence may exist. 
‘or the upper states of the alkali metals (Ky, Lig) there is no evidence for a decreased 
ate of convergence—indeed, the reverse is true. For the upper states of the halogens 
nd mixed halogen molecules, however, there is definite evidence of a decrease in 
he rate of convergence of the vibrational levels as the limit is approached. This 
ccurs for the upper states of the main systems of all the halogen molecules ‘for 
vhich data exist. 

The effect is well illustrated by ICl; the curve of AG(v+ 4) against v for the 
pper state is shown in figure 3, in which data are from Darbyshire (1932). The 
raph at first shows negative curvature corresponding to a slow but steady increase 
1 the rate of convergence of the levels; between about v= 20 and v=25 the graph 
hows an inflection and the ‘“‘ toe” turns up so that the rate of convergence 
.e. the slope) decreases to a value below that atv=0. ‘This tendency fora more 
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positive curvature to the graph in the region of the “ toe” would tend, if it occurred] 
for other molecules, to bring the dissociation energy a little higher than that which 
would be obtained by a graphical Birge-Sponer extrapolation using only point) 
before the “toe”? was reached. The effect must, however, be very small. Th d| 
laws of force which we have discussed will only apply at large values of r where tha} 
potential energy is already small, and the effect on the dissociation energy will} 
be almost negligible. Even in the case of ICI, which is one of the most extremal 
examples, the “toe” represents only about 0-04 ev., and the difference in dissocia 
tion energy obtained by using the true curve and by using an extrapolation fro | 
before the inflection (i.e. from around v= 20) would amount to only about 0:02 ev | 
There is little reason to think that this effect is important for any neutral molecule} 
As we have already pointed out, ionized molecules dissociate to one charged 
atom and one neutral atom. ‘The force between these atoms at large 
distances is likely to be that between a charge and an induced dipole, which is 
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Figure 3. Graph of A G(v+4) against v for the upper state of ICI. 


proportional to 1/r4. It is, therefore, especially likely that for ionized molecules} 
the curve of AG(v +4) against v will show an inflection towards a slightly highe \ 
limit. ‘There are no available data to check this supposition, but the fact that fo 4 
the lighter ionized molecules the linear extrapolation does not show any systema an | 
tendency to be too high may be connected with the existence of “‘toes”’ to the 
curves. If this is so, the dissociation energy of O,+, which depends on a good} 
graphical extrapolation for the excited state, may be a little higher than is generall 
assumed, although any considerable change would lead to a discrepancy betwenal 
the spectroscopic and the electron-impact values for the ionization potential of Oy. 
For molecules in which ionic bonding is important there is no case in which the} 
data are adequate for the drawing of an accurate potential curve. In many cases, 
however, the low value of the linear extrapolation seems likely to be associated|} 
with the fact that the true potential curve lies below the Morse function curve in 
the region where r is a little greater than r,. 


§8. SUMMARY 


. Although the remarks contained here may appear to be a severe criticism of they, 
Birge-Sponer extrapolation, they should really lead to an increase in its usefulness}, 
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by giving a fuller understanding of the limitations of the method and the likely 
deviations from the value given by linear extrapolation. 

Where data are adequate for a good graphical extrapolation for a ground state, 
the result should be satisfactory. Where, for a non-ionic molecule, the lower 
evels appear to converge linearly, the true dissociation energy tends to be around 
30 % of that given by linear extrapolation, being perhaps a little nearer this 
sxtrapolated value for light molecules and rather lower for heavy molecules. 
When dissociation leads to an atom in a 1S state the levels do not usually converge 
inearly ; the true dissociation energy is then often not much ovet half the linearly 
xtrapolated value. Extrapolations for excited states should be treated with 
aution, especially when ionic forces are likely to be important ; satisfactory 
ralues may be obtained for the lower excited states of molecules composed of 
lon-metals, especially if a graphical method is used . 

For molecules with largely ionic bonding the linear Birge-Sponer extra- 
olation does not appear satisfactory; in a number of examples it gives a value 
vhich is too low. For ionized molecules there appears little evidence to suggest 
ny failure of the linear extrapolation. 

_ For most types of molecule the above considerations should assist in the 
stimation of the dissociation energy, taking the linearly extrapolated value as a 
asis, and should also help in the assessment of the reliability of the value. 

There are two groups of molecules for which data are still inadequate to give 
auch guidance. For metallic oxides there will usually be some tendency towards 
mic bonding. ‘The electron affinity of atomic oxygen is not quite certain, but 
efinitely exceeds 2ev. and may be nearer 3; the ionic bonding forces are pro- 
ably not very important, provided the metallic ionization potential exceeds about 
ev., and the linear extrapolation can then be expected to be a little too high, as 
sual; for the oxides of the more electropositive metals any figure based on linear 
xtrapolation should be taken as only provisional. ‘The halides of the elements of 
sroup II are another unsolved problem, as we have the complexities due both to 
mic bonding and also to the presence of an atom with 1S, ground state; in many 
f them, too, the value of x,w, derived from head measurements may be inaccurate 
wing to the systems being weakly degraded. The evidence for CaCl and CaF 
iggests that the linear extrapolation may be appreciably too high, as in the hydrides 
f Group II elements, but this may not always be so. 
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ABSTRACT. A new ultra-violet band-system, E-x, of each of the molecules Sift 
SiSe and SiTe has been observed in absorption at 800-1000° c., and vibrational analysd#i 
have been done. The D-x systems, already known in emission, have also been photographe} 
in absorption. These observations have extended the D—x systems of SiSe and of Si 
and have led to somewhat improved values of the vibrational constants. The results at 
summarized in the following table :— 


D 


Molecule State Vo We huni Vee (By, 
SiS E a PABA ~ 402-3. * = ~—1-74 0-025 1:2:e¥ 
pti 35028°8 512-0 2:35 = 0):045 2°3 
eS 0 PADD SERS — 6°6 
SiSe E 38505°9 308°8 1-95 ==0°032 1:3 
D 32450°3 399-8 SNES: — 2:0 
x (43) 0 580:0 == ies == 6:0 
SiTe E 33991 242 (3:63) (+0-13) (1-4) 
D 28661°8 338-6 = er) == PA 
x (12) 0 481:2 E30 — 5-6 


§1. INTRODUCTION 
S: SiSe and SiTe each give rise to one band-system (D-x) in t 


near ultra-violet region, which appears readily in emission in posttigl| 
column discharges. _ For several reasons, including the fact that two ned 
ultra-violet absorption systems of GeS (which is isosteric with SiSe) were know ' 
we thought it likely that the absorption spectra of the silicon compounds woul 
contain systems lying at shorter wave-lengths than the D-x systems already know | 
This has proved so : SiS, SiSe and SiT’e each show banded absorption in t | 
ultra-violet region, and the bands of each molecule may be assigned to two syste 
D-X and E-x. 
The D-x system of SiS was originally developed rather well in emission, and} 
the absorption measurements have added nothing to the vibrational analysis of 
Barrow and Jevons (1938). Development of the p-x systems of SiSe and of SiT'l} 
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in emission was, however, somewhat incomplete (Barrow, 1939), and the 
absorption work has supplemented the original vibrational analyses and has led 
to improved values of the constants. 

Vibrational analyses of the three new E-x systems have been carried out, and a 
preliminary account of the work has been published (Vago and Barrow, 1946). 
The detailed results are given and discussed in the present paper. 


§2. EXPERIMENTAL METHODS AND RESULTS 


The absorption spectra were studied at temperatures in the range 800- 
1000° c. by methods already described (Barrow and Vago, 1944). An account of 
the preparation of SiS has been given elsewhere (Barrow and Jevons, 1938). 
SiSe and SiTe were made via the magnesium compounds. A mixture of finely- 
divided Mg and Se or Te was heated gently, using silica as a diluent to moderate 
the violence of the reaction. ‘The products were mainly MgSe or MgTe mixed 
with silica: they were heated zm vacuo at 1100° c. for about three hours to form the 
silicon compounds. ‘l'hese were found to sublime only very slowly, and since the 
other components of the reaction mixture were not volatile, the crude products 
were used without purification. SiSe and SiTe proved objectionable substances, 
as, like SiS, they hydrolyze rapidly in moist air. 

A hydrogen tube was used as source of continuous radiation. ‘The spectra 
were photographed in (i) a Hilger medium quartz prism instrument (E.498), 
(ii) a first order of a concave grating spectrograph with a linear dispersion of about 
7-4.a./ram., and (111) in a large quartz Littrow spectrograph with a dispersion of 
about 2:4.a./mm. at 2500 a. Ilford Ordinary, Process and Q.i plates were used : 
they were developed in Process developer to secure maximum contrast. 

The absorption spectra of all three molecules showed two distinct band-systems, 
but, whereas not much difficulty was experienced in getting good spectrograms of 
the three D—-x systems and of the E-x system of SiS, strong contrast was never 
obtained on plates of the E-x system of SiSe; and the E-x system of SiT'e was even 
harder to photograph. These difficulties seemed to be due to the intrinsic weakness 
of the systems: they could probably have been overcome by the use of longer 
columns of vapour, but in fact the spectrograms of the new systems of SiSe and 
SiTe proved adequate for their vibrational analyses. A reproduction of the 
3-x system of SiS is given in plate 1: the other E-x systems are generally similar in 
yppearance. 

The measurements of the E—x systems are set out in the Delandres schemes in 
ables 1-3. Additional or amended measurements of band-heads of the D-—x 
systems of SiSe and SiTe are given in table 4: these confirm the original v‘, v" 
assignments. ‘There can be few doubts about the correctness of the analyses given 
or the E-x systems, since the ground-state differences are. also known from the 
-x systems. Unfortunately, however, the distribution of intensities in the E-x 
ystems and the populations in the initial levels (determined by the Boltzmann law) 
oombine to make the v’=0 progressions very difficult to observe. Indeed, we 
ailed to measure any band-heads of this progression on the first plates of the SiS 
ystem, and they were only found on later plates of more contrast of the absorption 
t somewhat higher temperatures. It is therefore possible, but, we think, 
nlikely, that we have failed to observe v’ =0 progressions in the E-x systems of 


iSe and SiTe. 
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Table 1. Measurements of band-heads of the 
E<x system of SiS 
Italic numerals denote wave-lengths in air (1.A.); 


Large Roman, wave-numbers im vacuo (cm=+); 
Small Roman, wave-number differences. 


15 | 47664 
209735 
320 
14 | 47337-9 
2111-80 
330°9 
13 | 47001-0 
2126-94 
343°3 
12 | 46657-7 
2142-59 
349°4 
11 | 46308-3 
2158-76 
3555 
10 | 459528 
2175-46 
360-6 
9 | 45592-2 
2192-67 
305°2 
8 | 45227-0 
2210-38 
369°4 > 
7 | 44857-6 744°5 44113-1 
| 2228-58 2266-20 
1) aya 375°6 
6 | 44482-8 745°3  43737-5 740°2 4299733 
2247-36 2285-66 2325-01 
378°7 376°3 375°5 
5 | 44104-1 742°9 43361-2 739°4 42621°8 
2266-66 2305-50 2345-50 
384°6 384°2 385:0 
4 | 43719-5 742°5 42977-0 740'2 42236:8 735°9  41500-9 
2286-60 2326-11 2366-88 2408-85 
384'8 386°8 387-2 385°8 
3 | 43334-7 744°5  42590-2 740'°6 41849-6 734°5 41115-1 
2306-91 2347-24 2388-78 2431-46 
391°4 390°8 39176 
D 421988 740'°0 = 414588 735°3  40723°5 
2369-01 2411-30 2454-84 
394°4 394°7 393 
1 41804:4 7403. 41064-1733 40331 
2391-36 2434-48 2478-7 
395 399 
0 40669 737 39932 
t 2458-1 2503-5 
VD 
a EE, 
v’>0 1 D 3 
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The electronic and vibrational constants for the three known states of thes} 
molecules are collected together in table 5. The new data on the D-x and E-]] 
systems of SiSe and of SiTe have led to slight changes in the values of the constant} 
for the p and x states of these molecules. The ground-state constants for Si i 
determined from the extensive D-xK systems were found to be in good agreemeni} 
with the figures found from the analysis of the E-x system. ‘The k-state data fo} 


Table 4. Additional or amended measurements of band-heads of the 
D-x systems of SiSe and SiTe 


SiSe SiTe 
v', v”’ v v, vu” v 
8,0 35431-6 cm.-? 7,0  30867-8 cm.7? 
7,0 055°1 8,1 700-1 
8,1 34854-0 6,0 549°3 
6,0 679-1 Tl 385:0* 
ph 476-0 ; 5,0 232-3* 
5,0 302-2* 6,1 070:8* 
6,1 101-1* 4,0 29910-4* 
Sil 33727°-1 Sal 752-4 
4,1 346-1 30 586:3* 
4,1 431-6 © 
2,0 257-4* 
Sri 106:9 


* Band-head also measured in emission. 


Table 5. Summary of electronic and vibrational constants 


Molecule State Ve We Ay. Veo (SA 

SiS E 42322-1 402:3 —1-74 —()-025 1-42 x 10% 
D 35028°8 512-0 —2:38 —0-045 Zo 
x 0 749°5 — 2°56 — 4-94 

SiSe E 38505 -9 308-8 —1-95 — 0-032 1-16 
D 32450°3 399-8 —1:93 — 1:96 
x 0 580-0 —1°78 _— 4:4 

SiTe E 33991 242 (—3-63) (+0-13) 0:79 
D 28661-8 338-6 =1-70 _— Ty55 
ye 0 481-2 130) — 3:13 


Note.—-The units are cm.~1 (in vacuo), except for «,, the force constant, which is givew} 
in dynes/em. The signs of the vibrational constants are those applicable when th 
vibrational energy levels are represented by 


[w(v+3)+x0(v+3)+y,0(v+h)3+ oe]. 
SiS and SiSe call for no particular remark. The expression given for the vibra} 
tional energy in the F state of SiTe is, however, inherently rather unlikely, as i 
includes unexpectedly large anharmonic and cubic terms. Errors in these terml} 
might be attributed to unsuspected perturbations, or to the incorrect location ¢ 
some of the band-heads (the Te isotope effect might be partly responsible} 
However, the expression given represents satisfactorily the experimental 
measurements, and at the worst the value of w, is not likely to be much in error. 


Vi 


§3. DISCUSSION | 
It was largely the existence of two absorption systems of GeS which led us td, 


look for second absorption systems of the present molecules: in particular} 
the expectation was that SiSe, being isosteric with GeS (48 electrons), woulc 


be spectroscopically rather similar. How close the constants are for thes¢| 
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ympounds may be judged from the following figures :— 


GeS SiSe 
State Ve Ke ye Ke 
E 4:82ev. 1:27K10® 4:77 1:16%105 dyne/em. 
D 4-08 1:85 a 2a" 1:96 
x 0 4-36 0 4-11 


he isosteric group SiTe, GeSe and SnS show analogous behaviour in regard 
) their electronic and vibrational constants. The rotational constants of such 
roups of molecules may also be expected to lie rather close together. 

Lacking rotational analyses, we can only speculate about the molecular 
ypes of the E states of these molecules. One possibility—and it is no more than 
1at—may, however, be suggested. The lowest known excited states of SiO 
ad of SiS are ‘II, and the presumption is therefore strong that these states (and 
robably the similar states of SiSe and SiTe) correspond to the lowest excited 
nglet state of CO, alll. Now, although the bond in unexcited CO is much 
ronger than, for example, that in unexcited SiS, it might perhaps be expected 
vat the relative change in bond strength in going from the unexcited ! level 
) the *II level should be similar for the two molecules. ‘To test this possibility 
e set the force constants for the ground states at 100, and calculate those for the 
xcited states as percentages of the ground-state force constants. The figures 
re given in table 6. 


Table 6. Force constants, relative to those for the ground states 


co SiO SiS SiSe SiTe 
AHI 48-8 11 47-0 pil 46:8 Dp 47-7 D 49:5 
a’ 8X 29-7 ae E 28-8 E 28:2 E 25 
a*tl 64:3 = as = = 
x 1+ 100 1y 100 x1X 100 x 100 x 100 


It appears that for the three molecules CO, SiO and SiS the force constant 
wr the !II level is about 48°% of that in the ground state: the values for the D levels 
f SiSe and SiTe suspected to be !II are also close to 48°%. The figures for the 
levels are much lower and agree fairly well with that for the a’3% state of CO. 
‘his measure of agreement may be quite fortuitous—it is noteworthy that the 
35) level lies somewhat below the a1lI level in CO, whereas the E states for the 
licon compounds lie a little above the corresponding D states; but for what it is 
orth it suggests that the E states of SiS, SiSe and SiTe are of the type °B, 
valogous to the a’ 3X state of CO. 

The energies of dissociation of the silicon compounds are of interest. We 
ave, it is true, no very precise information, but only figures derived from long 
<trapolations of the vibrational levels: nevertheless some tentative conclusions 
in be reached. Let us consider SiS (see also Barrow and Jevons, 1938). The 
rergy of dissociation for the ground state, Dx, calculated by the Birge-Sponer 
‘trapolation, is 6-72 ev. The values for the D and E states, calculated in the same 
ay, but including the cubic terms, are 1-75 and 1-36, giving (vy. +D),=6-09 
id (vy.+D),=6:61. Now it may be presumed that the ground state of SiS 
formed from Si and S atoms in their ground states, *P. It is also most probable 
at the p1II state dissociates into two °?P atoms, since the next combination 
"atomic states that could give a singlet molecular state is, by the Wigner-Witmer 
les, 1D +1D, which lies at about Zev. above *P +P. If the combination 
)+1D were involved, we should have, using Dy = 6-7, (v.+ D)p = 8-7 or Dy = 4-4, 


36-2 
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which seems hopelessly out of agreement with the extrapolated value of 1- H 
Thus it is most likely that the D2II state of SiS also dissociates into Si(@®P)+ SCF 

The evidence about the E state is much less conclusive. It may be a trip 
state, in which case it could arise either from *P+*P or from 3P+1D ato 
states. However, the figure for (v,+D),, 6-6 ev. is close to the mean || 
(v,+ D)p and Dx, 6-4, so that the dissociation SiS(z)>Si(?P) + S(?P) is feasibll| 
and the data for SiSe (which molecule may be presumed to behave similar 
to SiS) suggest that this is the most reasonable hypothesis. | 

The extrapolated values of D for SiO, SiS, SiSe and SiTe are given in table 


Table 7. Energies of dissociation (ev.) 
Molecule: SiO SiS 


State Dextrapolated Vet D_ Dyean Dextrapolated vet D Dyean 


B = Stay 1:36 6-61 1:22 

D ae 27 (2505 1:75 6-09 2-13 

x 8-26 8-26 8:26 6-72 6-72 6:47 
6-47 


Molecule: SiSe SiTe 


State Dextrapolated vetD Dymean Dextrapolated ve+D Dryean 


E 0:89 5°67 (e247 — — 1:37 

D 2 OM 6°59 2°02 2-09 5:64 2-03 

x 5°86 5°86 6:04 Sy BeEYS  SoShe: 
5 6:04 5°58 


On the basis of the arguments outlined above, we have assumed that each of the} 
and E states of all these molecules gives a pair of unexcited atoms on dissociatio 
and have then determined mean values of D,. These mean values have be 
used to calculate the improved values of D, and D, given in the colum: 
‘““Drmean- For the SiO(pD) and SiTe(E) states, large cubic terms make extr 
polation impossible: the figure given for SiO(E) comes from a rough predicti | 
of that level (see below). ‘There may be considerable errors in these figures, bi 
they probably indicate correctly the general trends. Comparison with the fo 
constants given in table 5 tends to show—reasonably enough—that in a given sta 
the most tightly-bound molecule in the equilibrium position often requires all] 

| 


the greatest amount of energy to dissociate it into isolated atoms. Similar 
for any single molecule, the x state is more stable in both senses than the D stat 
which is in turn more stable than the E state. 
Finally, it would seem to be fairly safe to predict the existence of an E-X syste] 
of SiO at shorter wave-lengths than the II—1X system. Crude extrapolati | 
puts the E level at about 6-6 ev. above the ground state, so that the 0,0 band 
this system should lie at about 1850 a. The vibration frequency in the E stall 
should be about 670 cm-* and the energy of dissociation ca. 1-2 ev. 
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JRDER ABERRATIONS OF OPTICAL SYSTEMS 


TANGENTIAL ABERRATIONS OF A SYSTEM OF COAXIAL SPHERICAL 
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BSTRACT. ‘The behaviour of an optical system is known if the position of the point 
‘intersection of an arbitrary ray passing through the system with any chosen plane of 
ference can be determined. In the present paper tangential aberrations of any order 
e defined in terms of the distance of this point from an ideal image point, using elementary 
ethods only. It is shown how these aberrations may be computed in practice without 
ie use of any trigonometrical traces. 

As a practical example, the computations necessary for the determination of the 
mmplete primary and secondary aberrations and of the tertiary spherical aberration are 
ven in detail for the case of a Cooke Triplet. The information so obtained is compared 
ith that yield« d by strict trigonometrical tracing, illustrating the good agreement between 
em and the simplicity of the computations involved, 


§1. INDEX OF SYMBOLS 


HE following list gives most of the symbols employed in the text. Each 
[ nme in brackets refers to the paragraph in which the symbol first 


occurs. 

Symbol Meaning of symbol 
= coefficient of y?-”u” in the expansion of v’. (11) 
a (suffix) extreme ray of pencil. (3) 
m a quantity linearly related to y and u. (4) 

8 
A 2a, pi (12) 

v=1 

ledge A: 
A. habe (13) 
b, coefficient of y®-u” in the expansion of v’, (11) 
b (suffix) extreme ray of pencil. (3) 

5 
B DD Dinka 2) 

v=1 

ides 
B 0 (13) 
eC Nir. (4) 
ot Wir. (14) 
oe coma. (13) 
yy terms of I’, and T,. (8) 
a’, terms of the expansion of AA*. (5) 
d distance of diaphragm from preceding surface. (10) 


d’ distance between successive surfaces. (15) 


546 


H. A. Buchdahl 


SA) | 
v=1 

SAAS) 
xy piso) 


- ‘ : { 
coefficient of z2¥+1 in the expansion of A. (7) 
increments in y, win the equations giving sin Uj, etc. (5) 
increments in y, w in the equation giving sin Uj, etc. (5) ] 


increments in y,... in the equations giving tan U;,... ete. (14)) 
A(sinU+sinI). (4) 
X;—X;. (1) 

ee 


radius of aperture of diaphragm. (10) 

final surface. (5) 

posterior focal length. (10) 

paraxial focus (point conjugate to axial point of object). (3) 
coefficient of ¢” in the expansion of »/(1—t)(1—xt). (6) 
terms of the expansion of A. (5) 

height of object. (11) 

paraxial image height. (3) 

(Ii, Lj) tan Uj. (3) 

paraxial angle of incidence. (1) 

angle of incidence. (1) 


CS (9) = 
N/N’. (7) 


paraxial axial intersection length. (1) 

axial intersection length. (1) 

invariant of Lagrange type. (4) 

N(-—L)usinU. (4) 

N(i—L)utanU. (5) 

refractiveindex. (1) 

quantities refering to the ray (1,0). (4) 

coefficients in the expansion of y of principal rays. (10) 
quantities referring to the ray (0,1). (4) 

radius of surface. (1) 

arbitrary conjugate rays. (3) 
associated ray of R. (4) 

coefhicient of up? ’u” in the expansion of tan Us! (13) 
current symbol in computing scheme. (15) 
paraxial inclination of ray to axis. (1) 
=n 16) | 
inclination of ray to axis. (1) 
Hohe «(3) | 
semi-aperture of incident pencil of parallel rays. (13) | 
(L—njr. (7) 

kx. (7) 

eG) 
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y =. (5) 

= sin. 1(7) 

O(sufhx) quantities referring to focal paraxial ray. (5) 

; denotes any quantity behind surface corresponding to the ‘un- 

dashed” quantity before it. (1) 

(a, 8) a paraxial ray determined by y,=«, u,=f. (4) 

[x, B] a ray determined by L, sin U,=«, sinU,=f. (4) 

The quantities L, U, I, N, 1, d’ are the same as those used by Conrady (1929), 
so that the equations determining the passage of a ray through a surface are: 


ySiiea(i7e 7 sity Uae Pe eS (151) 
CN SID ae ey og ok (12) 
INCU sD Ve ie Se Se ease (153) 
7 Sindy ( L497) sindicdy |) gli hl ieee aise (1:4) 
ads) 0) Thies oe el neste (1:5) 
and similar paraxial equations, i.e. 

Fi Ent Ee Ag F as (1.6) 

ACN] es Oya ye Be ee Pe 88h oo (1:7) 
INGEt) = Os NP ON ie ees (1.8) 

Te (Lee ae att a aia’ dice (1.9) 
AShea— id ee” eer 1 ing Lares (1.10) 


It should be noted that H’ is positive if the point of intersection of R’ with 
the plane of reference lies below the axis of the system. ‘The phrase “ correct 
to the mth order ’’, as applied to any expression, is understood to mean that any 
errors involved in it are of at least the (+ 1)th order. 


§2. INTRODUCTION 


It has been stated occasionally (Conrady, 1929, chap. 6 ; Hardy and Perrin, 
1932, chap. 6) that the algebraic determination of the aberrations of lens systems 
ther than the primary involves computations of such magnitude as to make 
such analysis useless from a practical point of view. ‘Thus Conrady (p. 255) 
tates: “‘...due consideration of the disturbing effects... leads to the vast 
complication of the problem of determining the true secondary aberrations 
which up to the present has rendered the complete evaluation of the latter a task 
yf almost hopeless complexity’. Confining ourselves to the monochromatic 
angential aberrations of a system of coaxial spherical refracting surfaces 
ve shall show in the present paper (after suitably defining the “ aberration of 
yrder m”’ of an arbitrary ray passing through the system) that the computation 
»f the secondary aberrations represents relatively little labour in this case (see § 15, 
able 1); and that even the computation of the tertiary aberrations would not 
n practice be a task of insuperable magnitude. The methods used are essentially 
lementary in character; being entirely algebraic, the computations require 
10 trigonometrical tracing, and are peculiarly adapted to the use of calculating 
nachines. The simplicity of the method rests partly upon the choice of suitable 
variables, and partly upon a procedure which avoids the necessity for the explicit 
lerivation of any summation theorems. 
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§3. DEFINITION OF THE ABERRATIONS 


(a) Let an arbitrary ray R passing through the system from a given objec t 
point be defined by the values of Lsin U and sin U at the first surface (see § 4 (d)).4) 
A paraxial ray from the axial point of the object (the latter being perpendicula 
to the axis) intersects the axis in the conjugate point F’ in the image space. The 
plane through F’ perpendicular to the axis (i.e. the paraxial focal plane) is inter 
sected by the emerging ray R’, conjugate to R, at a perpendicular distance H" 
below the axis. If the system were perfect, all possible rays R’ from the given 
object point would have to pass through the ideal image point lying in the paraxial} 
focal plane and a certain distance h’ below F’, where h’ is proportional to thet 
height of the object. Let 


v= oh (3.1) 


v’ may simply be called the aberration of the ray R, v' being regarded as a function 
of the pair of variables L, sin U, and sin U,. (Regarding the height of the objea 
as variable, these two variables are independent.) | Now clearly, if the sign of Uy) 
be reversed, the signs of sin U,, L, sin U,, and of v’ are all reversed. Accordingly} 
v’ may be expanded in a series of homogeneous polynomials of the third, fifth,... . | 
(2n+1)th degree (n=1, 2, 3,...) in the independent variables L, sin U,, sin Uj, 
We define the aberration of order n to be the polynomial of degree (2n + 1) in this} 
expansion. ‘Accordingly the primary aberrations involve four independent} | 
coefficients, the secondary six, and so on. Hence the determination of the} 
aberrations up to and including the mth order requires the evaluation off] 


x (2v+2)=n(n+3) independent coefficients. 
vy=1 


It may be useful slightly to extend the definition of v’. Instead of considering 
the intersection of R’ with the paraxial focal plane we may consider the inter-| 
section, H’, with some arbitrary plane of reference (which is, of course, perpen-! 
dicular to the axis). ‘Then H’ may also be expressed as a series of homogeneous} 
polynomials of degree (2n+ 1), (n=0, 1, 2,...). We define v’ in this case as H'| 
minus the polynomial of the first degree. Our previous definition is contained} 
in this as a special case. 

(>) The four primary coefficients are, of course, intimately related to the} 
primary aberrations of von Seidel. The actual relationships between them may 
be established by defining the spherical aberration, distortion, tangential coma} 
and curvature of field purely geometrically and then expressing these in terms| 
of v’. At the same time these definitions may, if desired, be retained for practical} 


purposes on proceeding to higher orders. Accordingly we take as a measure} 
of the | 
(1) spherical aberration the distance below the axis at which a ray from the} 

axial point of the object strikes the paraxial focal plane; 
(ii) distortion the value of v' corresponding to principal rays, i.e. rays} 
passing through the centre of the diaphragm. 
If the extreme rays of a pencil of rays be distinguished by the letters a and bt 
respectively (the extreme rays being determined by the rim of the diaphragm), then} 
(iii) coma is measured by the length of the line drawn from the intersection} 

of the a and b rays to the principal ray, parallel to the focal plane; whilst the 
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(iv) curvature of field is the perpendicular distance from the intersection 
of the a and b rays to the focal plane. 

In the course of the evaluation of v’ to the nth order, the direction of the 
ray in the image space is in any case determined to the (n—1)th order, so that 
all the data necessary for the calculations (i) to (iv) are available. By far the 
most important of these in practice is thé coma, partly because of its connection 
with the question of unevenness of illumination in the image, and partly because 
it is not so easily found by inspection of the values of v’ corresponding to different 
rays of a pencil (see § 5), 

It should be observed that in the definition of v’, as well as its computation, 
the position of the diaphragm is not involved. The variation of the aberrations 
with shift of diaphragm may therefore be investigated after the main computations 
have been carried out. 


§4. ASSOCIATED RAYS 


(a) With any given ray R=[L,sin U,, sin U,] we may associate a paraxial 
ray Ro=(L,sin Uj, sin U,), i.e. a paraxial ray the nominal initial data of which 
are the same as the actual data of R. We shall say that R and R, are associated 
rays. (The particular choice of the data determining R, is dictated by the 
equations of §4 (d)). 

Let « be any quantity connected with the paraxial ray Ry which is linearly 
related to y, and w,. Let A be a quantity connected with the associated ray R 
which reduces to « in the (paraxial) limit. Then clearly the difference A—« 
will be small, and it is easily seen that it must be of at least the third degree in 
y,,u4,. (A is not thus uniquely determined; for instance, taking <= u, A might 
represent U, or sin U, etc.) 

Now one of the main difficulties of higher-order theory is the fact that we 
cannot simply replace any particular A by the corresponding « of the associated 
ray, as is done in primary theory, since such a procedure would make all but the 
primary terms meaningless. Accordingly it is desirable to find, if possible, 
expressions for the differences of the type A—«a, as explained above, for sub- 
sequent use, at the same time choosing the most convenient explicit form of A 
sorresponding to the particular «. 

In the following we define a surface invariant to be any quantity X; obeying 
the equation 


Mi Ne ANP ESO A Meciedeeeaef jue 1 weeny (4.1) 


dj é) 
X, shall be called a transfer invariant if 
Xj41—-Xj =A*X,=0; Claes f lye Ae eet (4.2) 


J 


whilst if X; satisfies both (4.1) and (4.2) simultaneously it shall be called a (true) 
yptical invariant, for which evidently 
Ajp=A,=Xy, for all7,7,k. = nae (4.3) 
Moreover, it will be convenient to use the term quasi-invariant for any quantity 
which reduces to an optical invariant in the (paraxial) limit. 
Now if we consider two arbitrary paraxial rays, say (y,, #,) and (y,, 4), then 
he expression 


NECN (ea ees oe (eee Vibert (4.4) 
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is an optical invariant. For A*A,;=0, obviously; and, at any surface, 


AA=AN{(I—r)—(—1) ui = ANr (tu — i), by (1.6). eee (4.5) 
*, since ANri= Ac =0, 
ar=cd(a-74) =ca(usi-F*7) -0, Subse (4.6) 


since A(u +2) =A(i/i) =0, by (1.7), (1-8). 
Now let A be a quasi-invariant which reduces to A in the limit. Then the 
fact that the expression 


if 
Dei SAA a): 1 eee (4.7) 


v=1 
does not vanish identically is implicitly connected with what may be called the} 
imperfection of the optical system. A heuristic examination of such quasi 


invariants may therefore be expected to yield a method whereby the desired 
differences, A —« (of which the difference A —A is itself merely a characteristi 
example), and thus the aberrations of the system, can be determined with reason 
able ease. 


(b) The equation (1.1) suggests the following explicit form of A: 
ASIN Lusi 0,  -) ye ay eee (4.8) 
where (y,,#,) and [L, sin U,, sin U,] are arbitrary rays. ‘Then at any surface 
AA =cA(sin U+sinI)=NiA(Lsin U), 8 ee (4.9) 
since AA =AN{(I—r)—(L—r)}usin U 
=ANr(zsin U—usin I), by (1.1), (1.6) 


=cA(sin U-— ssin1), by(heZ) 


=cA(sin U+sinI— pat sin 1), 


which is (4.9), since A(u +7) = A(sin I/7) =0. 

The great importance of equation (4.9) lies in its simplicity, in the sense} 
that AA is split up into two factors, one of which depends upon y, and mm) 
only, whilst the other factor is independent of these. Since the expressio 
A(sin J+ sin U) recurs frequently below, we shall simply write A for it, whenever} 
confusion is not likely thereby to arise. Thus 


AN =A. 30 ee ree eee (4.10) 


Equation (4.10) suggests that the adopted expression for A, (4.8), is the most| 
convenient for our purpose. 


Since, by definition, « depends linearly upon y, and u,, we may write 


Chg rg Va te Clay ee teas (4.11) 


where the «,;, % ; are constants of the system. These constants are easily 
obtained in any given case. It is only necessary to trace the two (paraxial) | 


rays (1,0) and (0,1) through the system. Then the actual value of a; in the forme 
1S %pj;, Whilst the value of «, in the latter is aj 
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We have, by (4.10), since A*A,=0, 


j—1 

Ni —lasin USN enon Ue eA. 

or eae 
. . . . 7a 

N,y; sin U; — Nju,L, sin U; = Nyy, sin U, — N,u,L, sin U, re Ne 
v=1 


Reenter (4.12) 
Hence, since y, and u, are independent variables, (4.12) splits up, by virtue of 
(4.11), into the two equations 


: : ' at 
NjV pj Sin. U; — Njuy;L,sin U,=N, sin U, + ~ Co»An 


Ngigg Sin U — Nou 


Oa 


j—1 
.L;sin U;= —N,L,sin U, + © egy. 
vy=1 


ree (4.13) 
These equations may be solved for L,sin U; and sin U;.. Remembering that 
VC TANIA EDENS EE (4.14) 
Vai Mail 5 


which follows at once by considering the optical invariant (4.4) for the special 
case of the two rays (1,0) and (0,1), we obtain the simple result 


: ; ee I 
sin U;=u,,L, sin U,+u,,;sin Uy + ee Ny (Cpvttas — Covlly,)A,, 
v= 1 


and analogous equations for L,sin U;, sinJ;, sinJ,’, the only change being that 
Up; is replaced by y,,, Zy;, tp; tespectively, and similarly for u,;. Equation (4.15) 
may be rewritten more simply in terms of the associated ray of [L, sin U,, sin Uj], 
i.e. 
: et Se 
sin U;=u,+ & Wy (Covas— Cavtns) Avs Are (4.16) 
y=15V] 
and analogous equations for L,sin U;, sinJ;, sinJ,’. Notice that the equations 
for sin U;’ and L,’ sin U,’ are similar, except that the summations are extended 
from v= 1 to y=] instead of toj—1. 
We have thus found the desired equations giving the differences between 
corresponding quantities of associated rays. 


§5. PRINCIPLES OF THE PRESENT METHOD 
(a) Writing A* =A sec U we have, since A*A* =O, 
' f 
INFN Pte INNS oi) Wie cowie (at) 
v=1 


Adopting the restricted definition of v’ from § 3 it follows, on taking the paraxial 
ray contained in A* to pass through the axial point of the object, that 


Aft’ = N'uog(log— L’) tan Up = NyuopH” wn ee ee (5:2) 
and Nee Nii Ld UN gilt = UN cUoflt 5 els oes (5.3) 


where / is the height of the object, so that 
ve = 
Njugsv' = & AAF =D*, say. Perit Grey) 
r=1 
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"st 


5 


(We have distinguished quantities referring to this particular paraxial ray by 
the suffix 0, and we shall refer to it in future as the “ focal paraxial ray ”’.) 
Accordingly we may take the “contribution by the jth surface to the final 
aberration * to be A.A3* Nyuos - 

Now suppose (dropping the suffixes referring to the particular surface) 
AA* to be expanded in ascending powers of sin U (see § 8), say 


AAF=TP,+I',+0s+ -.., (5.5) 


where I’, is the term in sin?"*! U. I’, may be rewritten by absorbing sin?”** U in 
the factor that multiplies it. We have then [’,, in the form of the sum of a number 
of terms each of which consists of the product of 27+1 sines of one or more of 
the angles U, J, J’, as well as of certain quantities referring to the focal paraxial ray, 
Replacing the sines formally by the corresponding quantities of the associated 
ray, the latter being regarded as functions of y, and w,, we are led to the series 


Dy(y, 4) + Py, u) + Psy, 4) + «5 Beer Gh 


where I’,,(y,) is a homogeneous polynomial of degree 2n+1 in y, u. (We have 
simply written y, u for y,, %,, respectively, and shall do so below whenever this 
procedure is not likely to cause confusion.) 

Now clearly the errors committed in this process cause T,, T3,... to lose 
their definite meaning for our purpose since the terms neglected in I’, are of the 
fifth and higher degree, i.e. of the same (and higher) degree as T',, and so on. 
However, returning to equation (4.16), we see that it may be written 


: - 1 j =< 1 ee 
sin U;= tas a >> cae) + mC + = a) a eo 


Hence, putting 


and $y,’, du,’ similarly with the summations extended from v=1 to 7 instead of 
j—1, we have 
A;=%p(V + dy;) + &gi(u +8u;), Se fo 


where 4 stands for ay one of the quantities Z sin U, sin U, sin J, sinJ’, and 
each a corresponds to the particular 4. When A stands for L’sin U’ or 
sin U’, Sy and 8u must be replaced by dy’ and dw’ respectively in (5.9). 

The important feature of (5.9) is that the value of any of the quantities which 
A represents is given by an equation closely analogous to the paraxial linear 
equation (4.11), i.e. by giving y and w on the right-hand side of the latter the 
same increments Sy and du, irrespective of which of the four alternative expressions 
A happens to represent. Hence, assuming dy; and du; to be expressed as functions 
of vy and w correct to the (~—1)th order (see below), we need only substitute 
y+8y; and u+5u; for y and w respectively in (5.6), rejecting all terms of degree 
greater than 2n+1, in order to make that series represent AA# exactly 
to the mth order. The terms of the resulting series may then be grouped 
together so as to give rise to a series of homogeneous polynomials in y, u of degree 
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3, 5, .... 2a+1. Then the polynomial of degree 2m+1 is the contribution 
(except for the constant factor 1/N’ui) by the jth surface to the aberration of 
order m. 

(6) We have yet to show how the dy, and 8w, are to be determined as functions 
of y and u. The simplest way of doing this, both in theory and in practice, is 
by means of a step-by-step method, following closely the procedure outlined 
above. For this purpose we imagine A to be expanded (at any surface) in a 
series exactly analogous to (5.5) (see §7), i.e. 


GEA Cars Cia C7nc a a a) wre (5.10) 
which again leads to a derived series, similar to (5.6), i.e. 
G,(y, 4) + Go(y, u)+Gely,u)+ ice. enews (5.11) 


Rejecting all but G,, dy; and du; may be calculated correct to the first order (i.e. 
3rd degree) by means of (5.8); thus 


igo 

ve eee (5.12) 
Sy —— Dy Un). 
Su, + ay 2 Cov Cwl > ¥) J 


Hence from (5.11) we can write down, using the expressions (5.12) at every 
surface, the terms 
Gij(y + 5y,u+8u)+Go;(y,U), nee eee (5.13) 

rejecting all terms of degree higher than the fifth. (5.13) is then the expansion 
of A; correct to the second order. This allows us immediately to calculate 
dy; and du; correct to the second order by substituting in (5.8). Hence with 
these values of dy; and 5u,, 

G1;(y +89; u + du,) + Goj(y + dy;,u+5u,) + Ga(y,u)  ...... (5.14) 
(rejecting all terms of degree greater than the seventh) represents A, correct to 
the third order. This process may be obviously repeated until we reach the 
expansion of A; correct to the (z—1)th order. (Notice that, of course, the only 
new terms to be calculated on passing from order m to order m+1 are those of 
order m+1, all terms of lower order remaining unchanged.) Then, finally, 


n 
2 Dy (YF OY, UO) awe (5.15) 
val 


(rejecting all terms of degree higher than the (21+ 1)th) is the desired expansion 
of AA# correct to the nth order. 

Alternatively we might calculate dy; and du correct to the mth order and 
then directly determine Af’. Though this method might be somewhat shorter 
in practice (particularly for large m), it would not yield directly the separate 
contributions by each surface to the final aberration. 


§6. -POLYNOMIALS 
The explicit expansion of A is greatly facilitated by the use of certain simple 
polynomials g,(«) defined by 


a/(1—2)(1 — xt) = 3 gn(x) eee (6.1) 
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where x, ¢ are independent variables whose range is suitably restricted so as to | 
ensure the convergence of (6.1). 
Replacing ¢ by t/x in (6.1) 
V(1-2]x)(1—2) = Eg, (xan & g,(1/x)e 
n=0 n=0 
SX) = #"S,(1/4). 


Hence ‘he, (4 V\y=x etx): | ORE ee (6.2) 


is a polynomial symmetrical in x and y. 
Expanding the square roots on the left-hand side of (6.1), it follows easily 


that eoraht F : 
(Gnice: Bilge Yond Jee biohe (6.3) 
gua) =(~1)" (2)F (=m, Bs 8-15 9), 


nN 


where F(«, 8; y; x) denotes the hypergeometric function of Gauss). 
Setting x= + 1 in (6.1) we have at once 


£,(+1)=0, (wT) 
Similarly £,(—1)=9, (n odd). 
The recurrence formula satisfied by the g,,(x) is 


(2+ 1) Sn+a(*) = (@ — 3)(% + 1) gn (*) — (4-2) aon a(%). nae (6.5) 


In particular 


go(x) =1, £1(%) = — 3(x + 1), 

Sox) = —3(e—1)?, go(x) = —79(x — 1)2(x+ 1), 

g4(2t) = — elo — 195002 + 60-45), pene 
galt) = — gh — Ma + 1)(7a2-+ 20 +7). } 


: §7. EXPANSION OF aA 
We have 


A=sin U’+sin/’—sin U—sin J, 
which, by (1.3), 
=cos UsinI cos lI’ —cos UcosI sin I’ +: sin U cos I cos J’ 
+sin Usin/ sin/’ —sinJ+sinI’ —sinU. 


che nee 7.1) 
For convenience, in this paragraph, ‘let a 
(L—nr)ir=x«; ~kx= Nx WN =; esin Ue: 
so that 
; Li (Md ip ly LGR” » LA eer ee (72) 
sinl=xz; sinl’=yz. 


Therefore, since cos /=(1—sin?J)}, etc., 


A=a{xy/(1 —2°)(1 —y*2?) —y/(1 — 2°)(1 — x22?) 
+ /(1 —x?2?)(1 —y?2")+ xy22?-x+y—-1},  ....., (7.3) 
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or, in terms of the polynomials of §7, 
A=a2{xyz*-x+y—-1+ 2 P8200") — Ve (a) + xo (yx) 2 toe ee (7.4) 
A=SscH (7.5) 
where ale 
8 = 2(%—y)(x+ I(y—1), 
82=3(%— (e+ Iy—- I -xt+y+a?+ay+y%), (7.6) 
83= 76% —)(x+1I(y—-1)l—x+y+s?+y?— 23 +3 + ot + xy 
+ x*y* + xy? + y4), 
and generally 
8n =a XE n(y") — VL, (x?) + 20”, nlR*), nZ2. 
Writing 6,=x«—y+1, we have in particular 
bn= 40,32 +y27+1)—B5% hee (7.7) 
It is evident that every 6,(n>1) has 6, as a factor, since 6,, vanishes when x = —1, 
y= +1, and k=1, in virtue of (6.4). 
Absorbing yn +1 in 6, we obtain the explicit form of the G,, G,, ... of 
(5.10); thus 
G,=3(sin/+ sin U)(sin J—sin /’)(sin I’ —sin U), 
G,=}G,(sin? U—sin UsinJ+sin Usin I’ +sin?/ 
+sinIsinl’+sin?J’),  { s*+++: (7.8) 


., etc. 


To obtain the G,(y,u) of (5.11) we neéd only replace sin U, sinJ, sinl’ by the 
values of u, 2, 2’ of the associated ray, giving 


Gy(y, 4) = 2 +u)(t—1)(r’ —u), 
G,(y, u) =4Gy(y, u{3(7? + u? +7’2) —u’?}, by (7.7), 
G,(y, u) =4G,(y, u)(u4 — 108 +10? + 70? +1? — But+ir tute 


4H! 4242493472) bo... (7.9) 
= WAG (9,0) + 4Gy(9, 0) uli? —P) — ia + =} 
., etc. 


where 7, u, ... are given by (4.11). Thus, for instance, 


Gi(y, u) a (e+ U) V1 a5 Us U) My A(t ae U)pV1 
+ (t—7') guy {(t —U)pVit (U — Ugly}. sooo (7.10) 
The case of a plane surface is easily disposed of. It is only necessary to 
calculate 7G, in place of G, (so that instead of the vanishing factor (7+) in G, 
we now have 7(i+u)=y, which remains finite), and to allow for this by sub- 
stituting Ni, for c, and Ni, for c, at this surface in (5.8), etc. 


§8. EXPANSION OF AA* 
We have 


AA* =A(A sec U) =sec UAA + AA sec U+(AA)(Asec U), 
AA*=secU’AA+AAsecU, = veser (8.1) 
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If we desire to carry out the computations to a high order, much labour will} 
evidently be saved if we express the I’, as far as possible in terms of the G, since; 
the latter are to be calculated in any case. This may be done as follows: 

By (7.5), remembering that 5)5=«—y +1, 

sin U3) 824s re (8.2) 
v=0 
also sec U=(1—sin? U)? 
sec U=1+42?+324+ 42°+... 
and 
sec U’ = 1 449,22? + (60,-+285)2* 2 ee (8.3) 
+ ($81? +8932 + $5971 + 7950°)2° +. - «J 
° AA* =cof1 + 485282 + (8951 + 350 Jatt... } {8122 + 8,25 +6,27+...$ +AA sec U 
= Cy{5132 + (Sy + 459251) + (53 + $5 75o + 5994? + 35945,)27 +... +AA sec U. 
hake (8.4) 

Combining the various powers of z with the factors which multiply them, and 
replacing the resulting sines of the angles by the corresponding quantities of the! 
associated ray, we obtain 


TD, =¢oG, + HAp(w’? — uv?) ] 
Ty =¢9(Go t+ 4u’2G,) +Ao(w’G, + 3u'4 — ut) 
T's = ¢p(G + $u"2G, +u'G,? +3u4G,) + 9(u'G,+ 4G? t sees (8.5) 
+3u'3G, + geu’® —u8) 
., etc. 


In practice [, and T, are.particularly important. Let us consider these in 
particular. We have by (8.5), (7.9) 
Dy = got + u)e—e')o" —u) + 32-7 )(ut wry, 
Also 
G+ $u?G, = G,{4u'? + 3(u? +7? +72) — Lu} 
= $G,(v? +7? +74 u’?) 

and 

Gy + 3(u'*— ut) = Bt —7')4u'(t+ uw)’ —u) + 3(utu'\(w+u)}, 
which may be rearranged, giving 


wGy+ 3(u'4—ut) = 30-1 ){3(ut+u')(w +2472 +u) + (¢+u)i—7’)(u' —3u)}. 
Hence if we write 
V1 = 2Co(2 + u)(t—2')(e’ — wv), 
Y2=3(0—1' (ut wu’ Ay = 3(u’ —u?)Ao, 
yee Rt A+ eye ares (8.6) — 
Ya= 4t(u' — 3u)Ap/co, 


P=" +72; } (8.7) 


we have finally 


P.=(1+Ye)¥s ty 


which is the form most convenient for computation. (Notice that at every 
surface Ay may be replaced by A), at once.) || 
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Correct to the second order, the contribution by the jth surface to Nyug,v' is 
then given by 


(¥15+ 723) 


10(y15 +2)? L O15 +24)? Cao 
ae + Yo, )¥e, +714 pg ee ele BOI LAR 8 Ld eee ee Op ae 
| Y1j V4). 5 Vij 43 N, Ou, cae Viv N, Oy; oe Y1 


§9. CHANGE OF POSITION OF OBJECT OR OF FOCAL 
PLANE 

Since the ray RR’ (§3) is quite arbitrary, it follows that the position of the 
object enters into the I’, of §8 only through the coefficients cy and Ay. In fact, 
using the extended definition of wv’ of §3, we are in every case only interested 
in the perpendicular distance from the principal axis of the point in which R’ 
intersects a chosen focal plane, the position of which is defined by Jj,, where 
the sufhix 0 is now used to refer to a paraxial ray having the required J} as its 
final axial intersection length. That is, if we are interested in the performance 
of the optical system for widely differing positions of the object (or, more generally, 
for widely differing positions of the plane of reference), we need only compute 
v, and wv, separately (i.e. v’ corresponding to the rays (¥91, %;)= (1,0) and (0,1) 
respectively). For then, by (5.4), 


if uf 
Pies DA + Hoy 2 BD. Oe ep itihs (9.1) 
j= j= 


Here y,, and u,, are connected by the condition that yo/ujr=[r. Hence writing 
(9.1) in the form 

, UpgU =YVolnfUp t+ Ug llaflg, ~— sn venee (9.2) 
and applying (4.11) and (4.14), it follows that 


/ 
tf 


Ny f / / / / x fF Hy 
O = Ry Ulos Has — Vas Mos Uy — (los Uns — Yt Wa ase veeeee(9.3) 
According to this equation, if the reference plane be moved through a distance 
8it7, v’ changes by the amount 


/ 


N, 7 
= Nt ig (Og —Wy olgp2 aw s'ee (9.4) 


oOo 
Q 


Notice that the actual magnitude of 6/; is not restricted to be small. 

The amount of labour involved in calculating both vj, and wv, is not much 
greater than that involved in calculating v' alone. We shall indicate how this is 
to be done most easily, considering only primary and secondary terms. (Higher 
orders are not essentially different.) Thus, from (8.6) we first determine in the 
ordinary way the quantities y,, y2, y; and y, with (V1, %)= (1,0), i.e. we deter- 
mine the four quantities y,, Y2,, 7p» Yap and hence the three products y,,7s, 
Yeps and 1,7 sp- 

Let c,/c,=9. Then, since 


¥ A,= ons VyAy/M ee rine tons 
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we have 
a> OV1p> 
Vag —ViY2p/M15 | 
VieVa= Vigan (9.6) 


V29g¥3 = —V1Y2p7s/ U4, 

V1¢qV4q = —V1Y 1pYap Uy. 
(Notice that the dy, already involve y,,.) Thus, having determined Uy, it Wi 
be seen (remembering that the dy; and 6u are unchanged) that the only ney 


multiplications required are those of the coefficients of y,;,y3 by 9; apart from 
these only a number of additions is required. | 


§10. THE DIAPHRAGM 


A pencil of rays is frequently specified by the position of the object poirj 

1T1 1 “ce ” | 

and the position and aperture of the diaphragm. ‘The “extreme rays”’ of t 
| 


pencil (assuming that there is no vignetting) are the rays grazing the rim of thy 
diaphragm. (The aperture of the latter is supposed to be circular.) For a give} 
position of the diaphragm there exists therefore a definite relationship betwee 
the values of y and wu (i.e. y, and u,) specifying principal rays. To determin 
this let the diaphragm be situated a distance d behind the surface immediate] 
preceding it, say the jth. Then, paraxially, the desired relationship is given b} 


jad = Lott Iait 


eb PI) 1 A 10.1). 
Upjyy + Ug U4 
Ug gy 
SS u, Say. eae tw ielie ls! 10.2 
ne ae Up, Pi y ( ) 


to the first surface.) However, we have strictly 


Lad VAYBYI TI gu 
u' + Up, du’ + uy,du' ’ 
which may be written 
y = pu — jy,’ + p10u,;'- eceene (10.3) 
‘To the first order the dy,’ and du,’ may be calculated by setting y=p,u, (8, a ij 
du,’ have already been computed for general values of y, and u,), the error 4 
committed being only of the second order. Hence, correct to the second orde 


Y=pietpy, = ae (10.4) 


where p; is a polynomial of the forth degree in p,. ‘This process may be continue 
to any desired order. Writing (10.4) in the form 


Yi =Py= (Lori Pr Pale eee (10.5) 


we see that the right-hand side of (10.5) gives the variation of the location wit 
respect to the first surface of the centre of the entrance pupil, due to the impe 
fections of the optical system, as a function of the inclination to the axis of tit 
incident principal ray. 
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Similarly the relation between y and u for the extreme rays, corresponding to 
10.2), is 


y=put Dee (1006) 


e 
Nab 
vhere ¢ is the radius of the aperture of the diaphragm. Higher approximations, 
orresponding to (10.5), may of course be determined. 

In the present notation, the stop-number of the system, defined as the ratio 
f the posterior focal length f’(=1/w,,) to the diameter of the entrance pupil, 
3 given by 

stop number=-222— “oi 
QU, fe 


§11. U, AND SPHERICAL ABERRATION 


(a) Let us assume that v’ has been computed (§§ 5-8) and is available in the 
nal form 
3 U' =ay? + ayy*ut agyu? + azu? + byy® + by y4u + by y?u2 
+ b,y?u? + byut+bsu'+higher terms.  ...... (Listy 


] 


11.1) here refers to the last surface of the optical system considered (which may 
f course be any intermediate surface of an actual system), the coefficients being 
valuated for a definite position of the object. We shall take the plane of reference 
) be the paraxial focal plane. 

To gain information about the system, which can be easily visualized, it will 
requently be desirable (though by no means essential) to plot v’ as a function of y 
ndwu. Ifthe object is situated at infinity it is advisable to plot v’ as a function of 
for various fixed values of u, whilst if the object, of height A, is situated at /,,, one 
aay plot v’ against u for various fixed values of h; y and w in this case are related by 
te condition 

(Jy, — L,) tan U, =const. =h, 


r in terms of the associated ray, 


By =) ua) (Man a (11.2) 
The position of the diaphragm may often be decided upon by inspection of 
1e graphs so obtained. 
 (b) If the object is situated at infinity (as is frequently the case with photo- 
raphic objectives) the spherical aberration is given by (11.1) on making u=0, i.e. 


Sper abe (Uoyee, = gly ucek | eae ee ere (11.3) 


‘his equation gives the radius of the disc of confusion projected by the system 
n a screen placed in the position of the focal plane. In practice we may also 
1easure the spherical aberration by the value of the expression 


[Op Ea sin Usleco, | hea 0 2 reese (11.4) 
hich is more convenient from the point of view of the computation of higher 


rders. 
Thus for the tertiary spherical aberration we require essentially the calculation 
f G,(y,u) omitting the term in u3; of G,(y, 1) rejecting all terms except those in 


37-2 
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y® and y4u; and of the term in y’ in G,(y,w). (Since the dyy and du; are then 
any case known to the second order, sin Uj’, and hence (v’),.9 may, if desired, 
calculated subsequently directly from (11.4). The difference between theg; 
results is, however, in general, quite negligible.) 

The great simplicity apparent in the calculation of the higher-order spheric) 
aberration ab initio, as indicated above, is there due to the fact that one of t 
independent variables, namely u, is zero. Hence to obtain the same simplici 
when the object is not at infinity but at /,,, say, we must proceed somewhat! 
differently, viz. by introducing the independent variable 


=i — lal ee eee (11.5) 
in place of y,. Equation (5.7) then becomes 


; Sag ee | j-1] 
sin U;= m6 = eo N, Gear le»)4») + (os tla) ( oe ay N, cso), 


We have therefore again a set of equations of the type (5.8), (5.9), etc., excep 
that where previously any quantity referring to the ray (0,1) appeared we no 
have in its place the corresponding quantity referring to the ray (J),,1). Accord 
ingly the computations now begin with the tracing of the two rays (/),,1) and (1,0 
And since one of the independent variables is now zero, the calculations are simild 
in every detail to those in the case where the object point is situated at infinity 
If, as will often be the case, the primary and secondary terms are already know’ 
it may nevertheless be worth while to calculate them again by the method abov 
For in this way a good check is obtained on equation (11.1), from which tbh 
spherical aberration is obtained by setting y=/)u, so that all the coefficien}s| 
appearing in it are checked in this way. | 
(c) In practice one will generally limit oneself to the computation of the exad#! 
primary and secondary aberrations; partly because the information so obtaineff 
will be sufficient in many cases, and partly because the evaluation of the genial 
tertiary (and higher) terms is rather lengthy. As we have seen, the computatio 
of the tertiary spherical aberration is, however, relatively simple. Consequenth 
it will frequently be advisable to carry it out. For obviously (taking the objed 
to be at infinity) when the inclination of the incident pencil (i.e. u,) is sufficientll 
small, the difference between the actual aberration and the sum of the primar 
and secondary terms is in general dominated by the next higher term in (11.1) n 
involving w; the latter being the tertiary spherical aberration. Similarly whejij 
io isottuite: . 
If the plane of reference is other than the paraxial focal plane, the considet 
ations above are not essentially altered. Thus, if the system were perfect, thi 
radius of the disc of confusion would be given by a simple term involving y « { 
u linearly. In the actual system the radius of the disc of confusion is, therefor | 
the sum of this linear term and v’ (extended definition). | 


i 


| 
; 


§12. DISTORTION 


From the definition of $3 it will be seen that the distortion is represented bi? 
v' when y is replaced by its value in terms of w as given by (10.4). We shal 
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estrict ourselves in this and the next paragraph to the first and second orders, 
he extension to all orders being obvious. Thus 


Distortion = (a p,° + a,p,? + agp, + ag)u3 
+ {(Bopr? + bypy4 + bop? + bap,” + bypy + b5) + Pa(3a9p12 + 2a,p, + ay) 0. 


Seaaete (1221) 
Therefore, if 
3 
A= ayp 
beer bie a EE ee ea: (12.2) 
nd B= Xi byp*”, 
v=1 J 
his may be written more compactly 
Distortion = Au? + (z + page). eee (12.3) 
1 


As might be expected, the term in p, will in practice often be found to be negligible 
‘compared with the other terms, in which case we have simply 


Distortion= Ane Bue (12.4) 


§13. COMA AND CURVATURE OF FIELD 


_ Inthis paragraph we take for the sake of simplicity an object situated at infinity. 
[his implies no restriction in principle. (For objects finitely situated we only 
equire the additional application of equation (11.2): see (13.5)). In equations 
n which several rays of a pencil are involved, quantities without suffix refer to the 
yrincipal ray, whilst, as before, the suffixes a and 6 refer to the extreme rays. 
It may be easily verified that according to the definition of §3, the coma, C’, is 
riven by 
Ce ("= Cpe) of b Va — Vy " tan U’ —(tan U,’ + tan “2 
Z 2 tan U,’ —tan U,’ 4 


vhilst the curvature of field, W’, is given by 


/ / 
ean A Sis 
tan U,' —tan U;,” 


v 


13.1) may therefore be written 
C’ ={h(u,' + v,')-—v'} — W’{R(tan U,' + tan U,’)—tan U’}. ...... (13.3) 
In this equation, since the lowest terms of the factor of W’ are of the third 
legree, W’ need in general only be evaluated to the (m — 1)th order if C’ is required 
o the nth order. Thus, when n=2, since dy,’ and dy’ are known to the first 
rder, the calculation of sin U;’ and hence of tan U;’ to this order presents no 
ifficulty. Therefore (omitting the sufhx /), let 
tan U! =up'ytUg Ut Uy (oy> +hyrut teyu?+tyu®). sss. (13.4) 
Ve have assumed the incident pencil to consist of parallel rays; 


~ ) 
Va Vu, i Se (13.5) 
Very — OB; 
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where w is the semi-aperture of the incident pencil. (13.5) applies becaug} 
y—Yal =(La1— 11) sin U,] is the perpendicular distance between the a ray 
the principal ray; and similarly for the 6 ray. (The fact that, if the pencil | 
actually limited physically by the diaphragm, the principal ray will not in gener : 
lie precisely halfway between the extreme rays at the first surface is of no practical 
consequence. ) 


Now to the required order, by (13.2), (13.4), (13.5), and (10.2), 


aes 2 i{(3aop1? + 2a, py + ay)U + ap} sven (13.6) 


and (tan U,’ + tan U,’)—tan U’ =u, '(3tpp, +t )uw. ....-. (13:79 
Also, by (11.1) and (10,4), 


}(vq' + Up’) —V' = (3app + ay )uw? + {(10b9p, + 6b, py? + 362p; + bs) 
+ 3appstu®w2+ (5byp,+b,)uwt, = sane (13.8) 


C’ = (3app, + a )uw? + {(10b p13 + 6b, p17 + 3bop, + b3) 
+ 3app3 — (3app1? + 2a, p; + ay)(3top, + ty) }uPw? 
+ §(5Dop1 + By) —g(BtoPy + ty)}ut#. we. (13.9) 


As in the case of distortion, it will frequently be found in practice that the 
terms in ps3, ty, t; may be neglected. If so, (13.9) may be written, in terms of 
(i222) 


C= Anus" + (Baw? By); eee (13.10) 
where 
A = pe \ 
ee a eee (13.11) 
1 a"B 
and Bae nl dp > 


By means of (13.2), a formula analogous to (13.9) may be written dows 
giving W’ correct to the required order. From our present point of view this is off 
no great interest (see end of §3)._ Accordingly we content ourselves with observing 
that if the terms containing pg, fo, t,, t can again be neglected, the equation for 


W’ analogous to (13.10) becomes 
Uy’ W’ = (A,u? + Aw?) + (But + Byu?w? + Bswt). ...... (13.12) 


§14. PRINCIPLES OF AN ALTERNATIVE METHOD 


It is interesting from a theoretical point of view that whereas previouslyj | 
we made use of the two quasi-invariants AA and AA*, each of which splits up intoy, 
two (quasi-invariant) parts, i.e. 


AA pc A} 
aCe \ ree (14.1) 
AA,=c,A, 
and 
AAS = ¢,(sec U’A + sin UA(sec U)) —c¥ sin IA(sec U), i 
Ber ac 14.2 
AA? =c,(sec U'A + sin UA(sec U)) —c¥ sin IA(sec U), | oar 


whe tN in = ae (14.3) |} 
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(together with their corresponding sums D,, , D, and D#, D# respectively), we 
can also proceed in a way which avoids the use ue the forte? tesaites 

If we now take the associated ray of the arbitrary ray R to be (£, tan U,, 
tan U,) we can write down a set of equations, by (5.1), similar to (4.13), but with 
the sines of the angles replaced by their tangents, and with the AA replaced by 
theAA*. Write 


1 i=t 
by;= — N. » (NINE 

v1 

ee ( tae (14.4) 
Sy= tap D AAS, j 


and similarly dy,*’ and du,*’, as before. ‘Then in place of (5.9) we have now 


tan U; = u,,(y + dy,*) + ug;(u + du,*), 
L,tan U; = Vpj(y + 89;*) + Yyj(u + 8u;*), 
sec U;sin J; = tj(y + 89;*) +24;(u + du,*), 
Js ys oer Petes 
AA} and AA} are now to be regarded as expanded in series of ascending powers 
of tanU. If, as before, we absorb, in the mth term of these series, tan2”+! U in 
the factor that multiplies it, then the mth term will be composed solely of the 
products of quantities given by (14.5). Hence we may proceed as before. 
Although AA is no longer explicitly contained in these equations, it is doubtful 
whether their practical application would be more convenient than that of the 
formulae developed earlier. 


§15. EXAMPLE OF COMPUTATIONS 


(a) The system used for the purpose of illustrating the practicability of the 
methods above is a slightly altered design of a corrected Cooke Triplet as given 
by Hardy and Perrin (1932, Chap. 21). The specifications of the system are 


Surface 1 @ 3 4 5 6 


r +20:05 |—128-3}—58-8 +18-9 + 311-3) —66-2 
Pm eOlG2 metals 72547 A, | 161G2/ Ci Bh (15.1) 
d' 3-896 1:63-0:93 13-42 | 3-03 - 
(f’ =96-729) ; d=4-0 (behind 4th surface) 
Max. stop-number of system =/"/5-6. 
The object is assumed to be situated at infinity, so that we are essentially computing 


vp. Moreover the factor Nyw'y, may be replaced by 1/f’ wherever it occurs. 
where f’ is the focal length of the system. 


(b) Table 1 

This gives the complete detailed computation for v’ (primary and secondary). 
Little comment is needed, as the computing scheme is arranged so as to explain 
itself. (A lowered number after a nought following the decimal point gives the 
total number of noughts following the latter.) A current symbol, ¢,, is attached 
to such quantities as recur at a later stage of the computation. The entries are 
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reduced to six significant figures (the original calculations were generally 
carried out to seven). In the case of a corrected system it is important that the 
number of significant figures carried should exceed by at least two the number of 
significant figures desired in the final result. Since no trigonometrical tables are 
required, this creates no practical difficulty. It may, however, be convenient at 
times to “modify” the dimensions of the system by dividing all the given radii 
and distances by the focal length (or a similar convenient quantity). This avoids 
the continual occurrence of a great number of noughts after the decimal point. 

Altogether it wili be seen that the computations can be rapidly carried out 
if a high-speed calculating machine is available. For each surface there are 
113 entries consisting entirely of simple additions and multiplications. Many 
of these can be further speeded up by observing that some of the coefficients 
repeatedly appear multiplied by various other coefficients. ‘Thus it will be seen 
that in the last part of table 1, t.,, for instance, appears 11 times, i.e. multiplied 
by ta; Ear tign taastans Laeshtazee tanta eC ele 

The first-order contributions by the different surfaces appear in the rows 
to1) tz) tog, tog; Whilst the second-order contributions appear in the rows fgg, ..., tyg 
(disregarding the factor f’ by which they should be multiplied to give the actual 
contributions as defined). Hence the final result is 


v’ = +0-0,145259y3 — 0-0,52286y7u + 0-18468yu? — 3-4308u3 
— 0-05;103856y° — 0-0,70536y4u + 0-0,93953y3u? 
+ 0-061374y2u*—1°2694yut—2-52210 2 es eee (45.2) 
(c) Table 2 


This gives the additional computations required to determine the secondary 
and tertiary spherical aberrations, measured by (l,¢— Ly’) sin U;’ (see (11.4)). 
There are only 33 entries per surface; fg) and ¢,, are the contributions by the 
individual surfaces to the secondary and tertiary aberrations respectively (except 
again for the factor f’). The final result is 


Sph. ab. = + 0-0314526y3 — 0-0;10463y> — 0-0,56284y". ...... (15:3) 
i f 
Since buy’ = (= tn) 9+ (z ta) * 
1 1 
oh ii 
and — dyf = (= ta) y+ (z tx) ye: 
i 1 


we can determine sin Us’ at once, by (5.9), correct to the second order. Hence 
the spherical aberration measured by (v’),_. may be calculated. Thus 
sec Uj’ = 1 + 0-0,53439y? + 0-0,20476y4. 
Sph. ab. = +0-0,14526y3 — 0-0;10386y> — 0-0,56813y7....... (15.4) 
The difference between the results (15.3) and (15.4) is negligible in practice, 
as will generally be the case. 


(d) Figures 1-5 


These are drawn, as explained in §11, on the basis of (15.2) and (15.4), i.e. 
v is plotted against y for fixed values of u. (The units correspond to f =96-729; 
i.e. they represent approximately percentages of the focal length.) In this way the 
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Figure 1. Sph. 


Figure 2, 3°. 
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values of v’ corresponding to a whole pencil of rays, incident parallel at the first 
surface, are graphically represented for the case of pencils inclined at 0°, 3°} 
6°, 12°, 18°, respectively to the axis. Roughly the whole field for which the 
objective is designed is thereby covered. 


Figure 4. 12°. 


Figure 5. 18°. 


The meaning of the various curves is as follows: 


Curve Tage cee ner , primary aberration. 

Curve II, ------ , (primary + secondary) aberration. : 
| F 1-5. 

Curve III, —— , curve obtained from a number of (Figures (oa 


strict trigonometrical traces. 


(primary + secondary + tertiary) sph. abn. 
(Figure 1.) 


(strict primy. + secndy.) abn. + terty. sph. abn. 
(Figures 2-5.) 
The values of y marked pr, a, b refer to the principal and extreme rays corres- 
ponding to the diaphragm with the given aperture and position. (As has been} 
pointed out, these are not involved in the calculation of Os) 


Curve IV, -—------ 
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(e) Distortion and coma 


The coefficients a o,..., 6; being known, the calculation of distortion, coma, 
etc., 1s so simple that its inclusion here is hardly necessary. Accordingly we shall 
content ourselves with quoting certain results. Thus, with the given diaphragm 
position, 


PLO OSZIW Mh ie ws Soe ay tiene (15.5) 
Putting y=p,u in the expressions for —Sy,’, +5u,4’, (given by (tg9)5,----, (ta7)s)s 
we have 
pg= 264445 8 eee (15.6) 


Substituting in (12.1), 
Distortion = —1-8444u?—8-6839u®, =... se. (1532) 


Notice the large value of the secondary term. 


DIST. 


20 


Figure 6. 


As has been mentioned, the term in p; is often negligible in practice. Actually 
n this case 


Dea = 2032370) ee eS (15.8) 
“or the coma we obtain, by (13.9), 
C’ = —0-0345580uw? + (0-027238u5w? — 0-0312949uw*). 


If the terms in Pz, fo, t, be neglected, we find 
C’ = —0-0545580uw? + (0-028985u3w? —0-0312741uw4), ...... (15.10) 
vhich, again, is nearly the same as (15.9).) 
We may note here that if a pair of extreme rays is defined by a particular 
alue of w, then the corresponding effective stop-number of the system is given by 
stop-number=48-4/e. - ss wegen (5.01) 
(f) Figure 6. ‘This is a graphical representation of (15.7). ea 
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(g) Figures 7-10 


~ 06 =°05. - O48 


Figure 9, 42°. Figure 10. 18°. 


Secondary coma. Thus for fixed u (at any rate up to U,;=12°), the coma is 
virtually given by a single term Proportional to the Sourth power of the aperture of} 
the pencil, whereas the Primary term alone would predict a residual coma pro- 
portional to the Square of the aperture. 

If a direct measure of the contribution to the coma by any individual surface 
is required, this may be obtained from table 1 by means of a formula analogous to 
(13.9) or (13.10), with the as teplaced by Ze1» ta7, tog, tog and the bs replaced by 
25g) 0++, tgs, Tespectively, 
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§16. CONCLUSION 


It has been shown how, by means of relatively short and straightforward 
computations, the tangential aberrations of higher order than the first can be 
obtained for a given system. In principle, the properties of the system are 
completely characterized for arbitrary positions of the object and the plane of 
reference (and of the diaphragm) by two expressions D* and D¥, or one such 
expression D* if the object and plane of reference have fixed positions. For 
the calculations of D¥ and D# it is convenient, though not essential, to make 
use of two expressions D,, and D,, similarly defined, the latter being needed cor- 
rect only to the order next lower than that to which the D¥ and D# are required. 
All four expressions are sums of certain expressions, of which one is contributed 
by each surface of the system, and all of them .vanish identically in the 
paraxial region. 

It is hoped to extend the present analysis to cover the case of arbitrary skew 
rays at a later date. 
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ABSTRACT. In viscometers of the Ostwald type used for the determination of 
kinematic viscosities relative to water, surface tension causes a reduction of the head 
available and increases the time of flow. The correction is greater for water than for 
organic liquids, but is not proportional to the surface tension. Earlier estimates of the 
correction are shown to be suspect, and the classical tables of Bashforth and Adams are 
applied to show how the capillary rise in a measuring bulb of specified size (3-85 ml.) 
and form varies during the discharge. ‘The calculated corrections for two liquids with 
surface tensions differing rather more than those of water and benzol diverge by 0:5% 
when a mean head of 10 cm. is assumed, and a simple experiment demonstrates that this 
estimate is probably low owing to effects of drainage on the shape of the meniscus. ‘Two 
methods are proposed for eliminating the correction during the calibration of viscometers. 


SMe tN RO DUC TON 


HE possibility that errors in measurement of relative viscosity may occur 
when liquids of different surface tension are tested in instruments of 
Ostwald (or similar) type has frequently been discussed (Barr, 1931) and 
much more frequently overlooked. ‘The calculation of the necessary correction 
for viscometers having bulbs of a given shape and size should be feasible when 
the surface tension is known, but several attempts have been based on inadequate 
theory. Since the blown bulbs of a U-tube instrument are not likely to be 
geometrically perfect, the method used by Bénard (1907), in which the effect of 
38-2 
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surface tension is measured at successive levels of the liquid in the bulb, would 
be more suitable for general use; this method is not, however, readily applicable 
with great accuracy to a completed viscometer, and could hardly be used with oils 
and opaque liquids that wet glass. 

The desirability of either avoiding or correcting for the effects of surface 
tension has increased in recent years with the demand for higher precision in the 
measurement of kinematic viscosity, especialiy of oils. As absolute determina- 
tions have been considered impracticable for routine work, there is fairly general 
agreement that the viscosity of water at 20° c. be adopted as a standard, for which 
the value 1-007 cs. is provisionally accepted, and that all viscosities be relative to 
this standard. The writer is in accord with this tendency, but it must be 
recognized that two disadyantages attach to the selection :— 

(a) If water gives a flow time of 200 seconds in a viscometer that may be 
calibrated by its means to +0-1%, an oil with a viscosity of only 1000cs. will 
require a flow time of over 55 hours. ‘This direct comparison may, of course, be 
avoided by calibrating a series of viscometers by a step-up procedure, but with 
three instruments, allowing a maximum time of 2000 seconds for the calibrating 
liquids of intermediate viscosity, the cumulative errors, for two determinations 
of the same accuracy in each tube, will amount to + V6x (071) 0-247 

(b) Water has an exceptionally high surface tension (72-7 dynes/cm. at 20°), 
while the value for most mineral oils is about 30 dynes/cm. If the surface-tension 
correction with water in a viscometer amounts to only 0-5°% and the correction is 
proportional to surface tension, the error in a comparison between water and oil | 
will be some 0-3°%; actually it will be shown later that the correction falls more 
rapidly, so that the error will be greater. 


§2. RECENT CALCULATIONS OF AN APPROXIMATE CORRECTION 
FOR SURFACE TENSION 

Mention may be made of three attempts later than those cited by ‘Barr (1931) 
to obtain an estimate of the magnitude of the effect of surface tension on the 
hydrostatic head in a viscometer. In the first (Barr, 1933) the “capillary” rises, 
h,, at various levels in a bulb of the shape specified for B.S. viscometers were 
deduced from a table due to Sugden (1921) and were subtracted from the 
hydrostatic head h operating during the outflow of successive elements of volume 
dV. Summation of dV/h and of 5V/(h—h,) over the whole discharge gave 
quantities proportional to the total time of efflux without and with correction for 
surface tension (a). For bulbs of 1 ml. and 4-7 ml. capacity the corrections were 
calculated as 21% and 0-8% respectively, for water fillings. Tis calculation is 
invalidated by the omission of any attempt to allow for the effect of the slope of 
the walls of the 45° cones on the capillary rise. 

In a long paper dealing with the use of his “ suspended level” in viscometry, 
Ubbelohde (1936) analysed the effect of surface tension in a spherical bulb on the 
assumption that, when the radius at the line of contact of the meniscus with the 
wall is r, the capillary rise is simply 20 cos 6/(rpg), where @ is the deviation of the 
wall from the vertical, g is the acceleration due to gravity, and p is the density of 
the liquid (more strictly p is the density of the liquid minus that of air). For the 
whole bulb he thus finds that the operative head is reduced by 20/(Rpg), where 
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R is the radius of the bulb. Relying on this value, he deduces a radius for a 
hemispherical “ suspended level” at the lower end of the capillary such that the 
correction should be eliminated. The calculation ignores the variations in shape 
of the meniscus during its passage down the bulb. The experimental 
comparisons of viscosities of benzol and water, determined (i) by applying the 
indicated correction to results obtained with a plane ‘“ suspended level”’, and (ii) 
by using the hemispherical level, differ by 0-7°% (or by nearly 2°% if kinetic-energy 
corrections are similar in the two series). 

Jones and Fornwalt (1938) used a method of summation similar to that 
indicated by Barr (1933) for the effect of surface tension in two viscometers : one (7) 
with a cylindrical length of 3-8 cm., radius 0-90 cm., and conical ends, and one (F) 
with a cylindrical length of 5-6 cm., radius 0-85 cm., and hemispherical ends: the 
_discharge in the first case was into a lower reservoir of radius 3-6 cm., and in the 
latter into one with a radius equal to that of the bulb. They assumed again that 
capillary rises are inversely proportional to radius, but they allowed for variations 
of slope of the wall by postulating that the rise was also proportional to 
(1—sin@)/cos#@; no explanation of this postulate is given. Since their 
calculations indicated corrections of 0-44°% and of 0-02% for comparisons of the 
viscosity of water and methanol in instruments T and F, respectively, and 
experiment gave a difference of 0:4% in the relative times of flow, they concluded 
_ that the correction was practically eliminated by the substitution of hemispherical 
for conical ends in the measuring bulb and the provision of a recipient bulb of the 
same diameter. 


§3. APPLICATION OF THE TABLES OF BASHFORTH AND ADAMS 


Gross errors may be inherent in any estimation of the effect of surface tension 
on viscometry when it is assumed that the capillary rise is proportional to surface 
tension and inversely proportional to density and to radius. Thus for water 
(c=72:7 dynes/cm., p=0-9970 gm./ml.) and benzol (o=28-9 dynes/cm., 
p=0-8774 gm./ml.) at 20°c. the capillary rises in tubes of three different radii 
rcm. may be obtained from Sugden’s tables as— 


0-683 cm. 0:273 cm. 0:067 cm. 0:0046 cm. 


0-280 cm. 0-090 cm. 0-010 cm. 0:0002 cm. 
Ratio B/W 0-41 0-33 0:15 0-04 


The last row gives the ratio of the rises for benzol and water in the same tube; the 
ratio of the values of o/p is 0-45. Now Sugden’s tables were based mainly on the 
much more comprehensive work of Bashforth and Adams (1883), who showed 
how the form of the meniscus dividing a liquid from air could be calculated in 
detail. The contour of a meniscus is defined by the value of a parameter 
B=bog/o, where b is the radius of curvature at the vertex (O in figure 1); tables 
(No. Il) are given for values of 8 up to 100, showing the ratio x/b and 2/6 all round 
a sessile drop, x being the radius at a vertical distance z from the vertex. ive 
and p are known, each value of f fixes a value for b and hence a series of values 
for x and z. 
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The points of which the co-ordinates « and z are thus listed are at intervals | 
of 5°, from 0° to 180°, of the angle ¢ between the vertical through the axis of | 
revolution of the drop or meniscus and the normal to the surface. Sugden’s| 
table refers to the special case  =90°, which applies to the meniscus of a liquid 
in a vertical cylinder—when, as is assumed throughout this paper, the contact 
angle is zero—but the complete work of Bashforth and Adams allows the 
capillary rise to be found for any small vessel of which the internal surface | 
constitutes a surface of revolution. The simplest form of bulb for calculation is 
the biconical, but the form specified for B.S. viscometers, viz., a cylinder with | 
conical ends, may be treated by finding the shape of the meniscus at successively | 
lower leveis in the bulb when the angle ¢ is first 135° (for the upper cone), then 


Figure 1. Diagram to illustrate definitions of x, z, r and ¢. The contour is that of Bashforth 
and Adams’ smallest drop (of mercury resting on the horizontal plane, shown dotted) ; 
it has been inverted so that the part below ¢=90° represents the meniscus of a liquid 
ae the wall of a vertical cylinder and touching a cone, of 45° semi-vertical angle, 
below it. 


90° (for the cylindrical part) and finally 45° (for the lower cone). For the small 
elements of volume enclosed between the fiducial marks on the tubing above and 
below the bulb, the angle ¢ is 90°, and the capillary rise may be found by either 
method. A constant a is defined by the relation a? =2o/(pg), where p is again 
strictly the difference in density between the fluids on the inner and outer sides 
of the meniscus, 1.e., in our case, the “apparent” density of the liquid in the 
viscometer. ‘The capillary rise h, is given by h,=a?/b; Sugden’s table gives 
r/b for various values of r/a, where r is the radius of a vertical tube containing the 
liquid, while we shall use the x of Bashforth and Adams’ tables to refer to the 
radius of the circle of contact between the meniscus and the bulb. 


Method of calculation for a spectfied bulb 


‘The surface-tension correction has been calculated for a bulb consisting of a 
cylindrical portion with a length equal to its radius, 0-896 cm., connected above 
and below by 45° (semi-vertical angle) cones to cylindrical tubing of radius 
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0-21 cm., with timing marks 0-42 cm. beyond the junctions of the tubing with the 
cones. ‘I'he form of the bulb is nearly the same as that which was shown by Barr 
(1933) to give the minimum ratio of wall surface to volume for a cylinder with 
conical ends; the length of the cylinder is assumed equal to the radius (instead 
of to 0-8 of the radius) since in practice it is convenient to lengthen it even more. 
The actual radii result from arbitrary selection of certain round values of 8 (1 and 
100) as minimum and maximum and of a value a?=0)-16 cm? corresponding with a 
surface tension somewhat higher than that of water (for which a2 =0-149 cm? at 
20° c.). The total volume of the bulb was then subdivided by finding the volumes 
discharged when f assumed the values 8 = 1, 4, 10, 40 and 100 as the meniscus fell 
in the upper cone; while the cylindrical part was discharging, 8 remained 
constant at 100 and in the lower cone it fell again, to B=1-5, before the junction 
with the tubing below. In the upper cone, particularly, the lowest point of the 
meniscus is so much below the circle of contact that a considerable addition must 
be made to the volume delivered from the frustrum above the plane of this circle; 
Bashforth and Adams provide a table (No. III) which gives the value of this 
addition, in terms of 6, for all values of 8 and of 6. For each element of volume 
discharged, AV, the capillary rise h, was assumed to be the mean between those 
characterized by the values of f at the beginning and end of the discharge of that 
element. The product AV ./, was finally summed for the whole bulb volume V 
and divided by V to give a mean f,. Assuming that the recipient bulb has the 
same diameter as the measuring bulb, the capillary rise in this arm will have the 
same value as that already calculated for the cylindrical part of the measuring 
bulb; the excess of h, over this value represents the correction to be applied to 
the mean head causing the flow. ‘The method of calculation may be illustrated 
by detailing the procedure for the first two elements of volume considered :— 

(a) The tubing above the upper cone has radius 0-21 cm. and the graduation 
is 0-42 cm. above the junction, assumed sharp. ‘The meniscus starts with its 
lowest part in the plane of the graduation (see figure 2, position I) and its shape 
will remain unchanged until the circle of contact reaches the joint (position IT). 
With a=0-40 cm. we have r/a=0-525, whence from Sugden’s table r/b =0-9165 
and b=0-229. The capillary rise in the tube is, therefore, h, =a?/b=0-698 cm. 
We require to know the distance of the vertex of the meniscus below the plane of 
the joint and also, for later use, the “ bubble”? volume v bounded by this plane 
and the surface of the meniscus. We find first the value of B=2b?/a?=0-655. 
Interpolation in Bashforth and Adams’ table II for B=0-655 and ¢=90° gives 
x/b=0-917, z/b=0-825, whence x=0-210 and z=0-189 cm. The volume AV 
discharged while A, remains at 0-698 cm. is, therefore, 7(0-21)?(0-42+0-19) 
=0-082 ml. Hence AV .f,=0-0590. Bashforth and Adams’ table IIi gives the 
volume of the ‘‘ bubble” as 1:47b? =0-0176 ml. 

(b) We pass next to the position where B=1; here b?=}a°B=0-08 and 
b =0-282, whence h, = a?/b=0-568 cm. For 6=135° the tables give r/b =0:7369 
and v=2-091b3, whence r=0-21 cm. and v=0-0473 ml. The meniscus has, 
therefore, become tangential to the cone at the junction with the tubing above it 
(position III) and the volume run out since the meniscus touched the joint at 
90° is 0-:0473 —0-0176=0-0297 ml. This volume is considered to have been 
discharged while the mean capillary rise was }(0-698 + 0-568) =0:633 cm. 
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The radius of the cylindrical part of the bulb was fixed at 0-896 cm. because}, 
this value of r results from the figure given by the table for B= 100 at 6=90°; the | 
meniscus touches the upper cone at the same time, when r=0-846 cm. for) 
$=135° (position VII). In this position the lowest point of the meniscus is || 
0-410 cm. below the circle of contact with the cylindrical wall and 0-536 cm. below) 
the circle of contact with the cone; a surprisingly large fraction of the total content| 


This minimum persists until the meniscus touches the lower cone (cf. figure 1), | 
after which there is a relatively rapid increase. The variations in capillary rise | | 
during discharge are plotted in figure 3 (upper curve). | 


Figure 2. Positions of the meniscus at four stages of the flow from a bulb shaped 
as in figure 3—not to scale. 


The total volume of the bulb between the timing marks, 3-85 ml., was sub- | 
divided into 13 elements in all. The sum X(AV.h,) divided by 3-85 gave | 
h,=0-1241 cm. If the discharge is supposed to be into a recipient cylinder of 
radius 0-896 cm., in which the capillary rise is 0-0566 cm., the mean amount by 
which the driving head is reduced is 0-1241—0-0566=0-0675 cm. Assuming |} 
that the average hydrostatic head is 10 cm., the surface-tension correction will | 
therefore amount to some 0-7°% when the liquid tested has the surface-tension 
postulated. ; 

As Bashforth and Adams’ tables do not extend beyond f = 100 it is not possible 
to complete a precise calculation by the above method for a liquid of lower surface 
tension in the postulated bulb. The positions of seven points on the lower curve 
of figure 3 were, however, calculated for a liquid for which a? =0-0625 (the value 
of a” for benzol at 20° c. is0-0671). These points refer to positions of the meniscus 
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in the upper and lower cones. The capillary rise in the cylindrical part of the 
bulb was found from Sugden’s table, and the extremities of the horizontal straight 
line were displaced in the direction of the hg axis by an amount arbitrarily 
estimated to allow for the flattening of the meniscus; the curve was then completed 
by the dotted portions. Measurements, with a planimeter, of the area between 
this curve and the ordinate h,=0-0057 and of that between the curve for 
a =()'16 and the ordinate h, =0-0566 gave a ratio 0:29 :1. The surface-tension 
correction when a*=0-0625 and the hydrostatic head is 10 cm. is, therefore, 
estimated as 0-29 x 0-68°%,=0:20%. The correction falls more rapidly than does 
a’, for which the ratio is 0:0625/0-16 =0-39. 

The correction for water in a larger bulb of the same shape is also not directly 
obtainable by the use of the tables. For a smaller bulb, of volume 2-3 ml. and 
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Figure 3. Calculation of surface-tension correction for viscometer bulb 
of volume: 3°85 ml. 


maximum radius 0-75 cm., with a recipient of the same radius, the correction is 
found to be 0-88% (instead of 0:68°%) for a hydrostatic head of 10 cm. when 
fe —0-16. 

Examination of the curves illustrates the correctness of the common assumption 
that an increase in the volume or in the maximum diameter of the bulb tends to 
decrease the effect of surface tension. They show, however, that an increase in 
the diameter of the tubing above and below the bulb will produce only minor 
improvement, the main portion of the effect being due to conditions in the upper 
cone. If adequate drainage occurs with a cone of greater semi-vertical angle than 
45°, it might be preferable to adopt an angle of 60°, say, for the upper cone, so 
that the meniscus more rapidly reached its maximum diameter. Observation 
shows that some improvement in the drainage from the lower cone is desirable, 
and the curves suggest that the capillary effect is here so small that a steeper slope, 
say a 30° cone, might be adopted without undue increase in the surface-t2nsion 
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correction. Such a change would require that the maximum diameter be} 
increased and the parallel portion practically eliminated so as to prevent} 
inconvenient increase in the overall length. A 60° cone and a 30° cone, base toy 
base, would need to have a maximum radius 29%, greater than that of the form) 
assumed in the above calculations, in order to have the same content, and would 


have almost exactly the same length. 


$4) CRIULLGISM OR Whe Vir OD 


A stricter method of calculating the correction for surface tension requires) 
that the height of the measuring bulb above the recipient be known, so that for) 
each element of volume, AV, discharged, a quantity AV/(A—h,) may be deduced} 
which is proportional to the time required for the discharge of that volume. ‘The} 
sum SAV /(h—h,) for the whole bulb is then proportional to the flow time} 
observed in an experiment. This method, which was used by Barr (1933) and | 
by Jones and Fornwalt (1938), takes account of the fact that the flow time is less; 
sensitive to variations in the driving head in the earlier than in the later stages of! 
the discharge. For a viscometer of normal design the simpler procedure now) 
adopted will produce only a small error in the correction, and the summation of 
AV .h, or the plotting of A, against volume allows examination of the behaviour} 
of different bulbs without reference to one specified mean head. 

A more serious defect is common to all calculations of this type. It has been} 
demonstrated in several papers (e.g. Jones and Stauffer, 1937) that the “ drainage’ 
error’, due to variation in the volume of liquid left on the walls of the bulb, may 
be serious when externally applied pressures are used to shorten the flow time, but} 
becomes negligibly small for viscometers of the type here considered, in which the} 
times allowed for drainage are, of necessity, proportional to the viscosity of the} 
liquid. Nevertheless, as the meniscus moves down the wall it must be distorted 


by the flow of liquid from the still wet surfaces above it and the distortion will Pe 


most severe precisely in those regions where it will cause the greatest effect, viz.,| 
where diameters are small. 


§5. EXPERIMENT TO ILLUSTRATE THE EFFECT OF SURFACE TENSION] 


A simple experiment has been made to demonstrate the necessity for a surface-| 
tension correction in the use of a No. 0 B.S. viscometer. After determining the| 
time of flow (1348 sec.) for water at 25°c. in the usual manner, the water was} 
blown up until the level was about 1 cm. above the upper timing mark and al} 
volume of 0-02 ml. of a solution of “‘ Wettol”’ (a wetting agent for photographic use| 
supplied by Messrs. Hford), of density 1-1010 at 25° c., was added at the surface. 
Flow was then allowed to proceed and the time observed was 1334 sec.,| 
approximataly 1%; less than that for water alone. During the flow the Wettol, 
tended to remain at the air-water interface, and probably none of it reached the: 
capillary during the timed interval. When equilibrium had been reached, thel| 
contamination of the water surface was shown by the greatly reduced rise of the: 
meniscus in the capillary above the level of the water in the other arm; the exalt 
rise observed, 2-5 cm., is not significant as the extent of dilution is unknown, but) 


pure water gave a rise of about 7-5 cm. (when the levels in the U-tube were lowered | 
to allow the full rise to occur). ; 
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The capillary rise for water in the 2-cm. diameter recipient bulb of a No. 0 
viscometer is calculated from Sugden’s table as 0-034 cm., and would have been 
reduced to less than 0-002 cm. if Wettol had been added to this surface also. The 
surface-tension correction with pure water is thus some 0-7°, greater than with a 
liquid having about 1/3 the surface tension of water. 

The correction calculated for water in a bulb of only 3-85 ml. capacity was less 
han 0-7°% ; for the 6-5 ml. bulb of a No. 0 viscometer it would be appreciably less, 
and it cannot have been completely eliminated by the addition of Wettol. Hence 
some factor, such as the distortion—suggested above—of the meniscus by liquid 
draining from the wall, must have been operative to produce the observed 
lifference of flow times. An attempt is often made, especially in viscometers, 
such as that of Bingham, which are operated by external pressure, to eliminate 
he effect of surface tension by constructing the recipient bulb to be of the same 
shape and size as the timing bulb; if the falling and rising menisci differ in shape 
he object cannot be completely achieved. 


§6. DETERMINATION OR ELIMINATION OF THE 
SURFACE-TENSION CORRECTION 

From the absence of any simple relation between a@?, 7 and h, it will be obvious 
hat the need for a surface-tension correction cannot be avoided by selecting such 
1 diameter for a cylindrical lower bulb as will give a capillary rise equal to the mean 
sapillary rise in the upper bulb: the matching can be correct only for one value 
of a®. ‘The same objection applies to Ubbelohde’s proposal (1936) to eliminate 
he correction, in his ‘‘ suspended level”’ viscometer, by causing the lower end of 
he capillary to expand into a hemispherical surface of calculated radius instead 
»f into a horizontal plane. 

‘Two alternatives appear to remain :— 

(1) 'To calibrate viscometers of the usual forms with two or more liquids of 
cnown viscosity and known surface tension, using such times of flow as will avoid 
rrors due to kinetic energy, and to state the different constants occurring for 
lifferent values of a2; and 

(2) 'To design viscometers so that the maximum surface-tension correction 
vill be negligibly small. 

(1) The former procedure presupposes that the viscosities required can be 
iccurately determined by other means. ‘The dynamic viscosity (poises) may be 
neasured in an instrument of the Bingham type, in which external pressures are 
ised that are 30 or more times the variation in hydrostatic pressure in the 
iscometer during flow. When all the numerous precautions and corrections are 
aken (including consideration of the ‘“ drainage error” which may be marked in | 
his type) accurate results are obtainable, especially if it is necessary only to make 
hem relative to the viscosity of water. ‘The higher pressures used necessitate 
maller diameters for the capillary, to minimize kinetic-energy corrections, than 
n the Ostwald type, and troubles due to suspended matter in the sample and to 
bstruction of the capillary are exaggerated. For the calibration of an 
)stwald-type instrument a further determination, that of density, is required to 
onvert poises to stokes. 

(2) If an instrument for comparison of kinematic viscosities can be produced, 
uch that times of flow are proportional to kinematic viscosity to the precision 
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required, this viscometer, though possibly inconvenient for routine determinations 
might be used to determine the relative viscosities required in method (1}} 
Cannon and Fenske (1936), and more recently Cannon (1944), have describeg| 
‘master’? viscometers for this purpose, in which the capillary length and meas | 
head are increased to some 45 cm. on the assumption that errors due to surfac} 
tension and to kinetic energy will become negligible. Their estimate of the orde! 
of the surface-tension correction (0-7 mm. difference in loss of head as betwee 
runs with water and oil) suggests that the error in the master instrument shoul 
be some 0-15°%. But in their fusiform 3 ml. bulb the effect is probably as muc 
as in the 3:85 ml. bulb, for which calculations on the lines indicated above woul 
show an error of more than 0-2°% ; and the experiment with Wettol points to thij 
being an underestimate. The error would be less if the large diameter of th} 
recipient bulb (3cm.) were reduced to equal that of the bulb, ss in the for 
described by Cannon for low viscosities. We propose to try the effect of usin; 
heads and capillary lengths up to nearly a metre: the correction should them 
become progressively less, and this will be tested by the variation in the relative 
times of flow of water and an organic liquid of lower surface tension as the lengt 
is increased. 
We have found, however, that it is possible to time with good precision thi} 
fall of a meniscus in a uniform tube of 1 cm. diameter on which two fine markt 
have been etched 5 cm. apart. With the same diameter for the other arm of thif 
U-tube the surface-tersion correction should be zero, apart from effects olf 
drainage on the shape of the meniscus. The importance of such effects will 
obviously be much less than in Cannon’s “‘ semi-micro”’ form, which has a bor 
of 2-5 mm. for the measuring and receiving tubes. ‘The occurrence or absence o 
sensible correction should again be demonstrable by increasing the length. Since 
the ratio of initial to final head is greater in this type of instrument than when wide 
bulbs of similar capacity are provided, it may be desirable, at least in the shorted} 
models, to make use of two recipient tubes in parallel, of the same diameter as th | 
measuring tube. ‘This arrangement, besides reducing the variation in head by : 
quarter and halving the effect of errors in loading, can be made to avoid the mai 
error due to maladjustment to the vertical of a U-tube viscometer; if the tube: 
are equidistant from the vertical capillary they are equivalent, as far as concerns 
the effect of tilt, to a recipient vertically below the upper bulb, as in Cannon and 
Fenske’s “routine” form. 
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ABSTRACT. In studying the absorption of liquids by rigid gels, certain thermodynamic 
elations are required in terms of stress components and generalized co-ordinates. These 
ave been derived using the analytical methods of Gibbs as developed by Guggenheim, 
nd the conditions under which they hold are clearly defined. The analysis includes a 
lew derivation of Kirchhoff’s relation for the heat of absorption, in which it is shown to 
ye an exact relation between the isopiestic and isothermal heat of absorption and the 
sopiestic and isosteric variation of equilibrium humidity with temperature, where the 
yressures concerned are the external pressures on the gel-liquid phase. The effect of 
nternal stresses existing in a homogeneous gel-liquid phase due to the absorption of the 
iquid is obtained by assuming that the energy of expansion of the gel is the same for 
‘xpansion due to absorption as for expansion due to external forces. The effect of external 
hear stresses On vapour pressure is considered, although it is not required inthe case of 
Lomogeneous gels, since the analysis is similar to that for the internal energy of expansion, 
ind the effect of such stresses is required when dealing with more complex structures such 
is wool cells. 


§1. INTRODUCTION 


ARKAS’ work on wood (1945), and in particular on the effect of natural 
restraints on the wood-water isotherm, has shown the importance of 
structural rigidity for the vapour pressure isotherms of gels. Wood 

s an aggregate of cylindrical cells with rigid sheaths which are partly filled with 
rel substance, and Barkas has been mainly concerned with the effect of the rigid 
heath on the water absorption of the gel substance; the restraints on the gel 
substance are, in this case, external restraints imposed by the sheath. ‘The 
estraints need not, however, be external restraints, for any energy of swelling 
yf a gel must be reflected in the free energy of the gel-absorbate system, and the 
sotherm will therefore depend to some extent on the swelling energy. Cassie 
1945) has already attempted to estimate the effect of this swelling energy on the 
water-vapour isotherm of keratin, but his work has been considerably restricted 
yy lack of adequate elastic data for the gel, and of appropriate relationships 
yetween the thermodynamic functions and the elastic parameters. 

The purpose of this paper is to establish the necessary relations of the thermo- 
jynamic functions and the elastic parameters, using Gibbs’ method. The 
irst requirement is a transformation from the load variables of thermodynamics 
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to the stress variables of elasticity, and from this the relations of chemical potentia} 
and entropy changes to changes of external stress are readily derived. Tha 
problem of the relation of chemical potential to the energy of expansion fo! 
internal stresses is, however, best approached through the intrinsic energy] 
term, and this and the heat of expansion are derived in terms of the swelling of 
the gel. ‘| | 

Gibbs’ method has the great advantage over the cycle method used by Barkag| 
that it avoids reference to swelling pressures, which are not easily visualizeq 
for the internal stresses in the gel itself due to its swelling. 

The analysis is only concerned with cases where the stress can be reduced 
to pressures along the three principal directions and with a fixed amount 0} 
gel with varying quantities of absorbed material; the state of the gel is thereford 
completely specified by the temperature 7, the number of moles absorbed, ny 
and the loads in the three directions X, Y, and Z, or alternatively the stressed 
X, Y, and Z, and where shear stresses are involved a rectangular block is com 
sidered. 


: | 
§2, RELATION BETWEEN VARIATIONS AT CONSTANT STRESS 
AND A'T CONSTANT LOAD 


From the thermodynamic relations expressed in terms of generalize | 
co-ordinates, expressions of the type 


(SB) rn = (0x/0A)xyz5X + (0y/0A)xyz oY aF (dz/0A)xyz OL ey. eee (1) 


can be obtained (e.g., equations (11) and (12) below), where x, y, are linear di 
mensions, A is either 7 or J and B either the chemical potential, u, of the absorbed 
substance or the entropy, S. 

It is required to express this in terms of (0x/0A)xyzand 8X etc., the fifth para} 
meter, determining the state of the gel, remaining constant throughout. Ti 
the stresses are small compared with the moduli, one can write 6X=yz6X and 
(0x/0A)xyz = (0x/0A)xyz without appreciable error: but in the general case 


Ox Ox Ox Ox 
ox = | = 5A ae 
g (a (se) «2X : eee. (2) ad” 
Bee @a) PD nl Dan Bal 
“"\@A/]xyz \OA/xyz \OX/4yz\0A te) Ge 


(35) ee(St ee? oe (2) 


but X=yz2X, 
J“. dX=y2dX + X(ydz+ 28y). 
(5) __¥ 1 (dz 
Saha Talo) eat @) 
also (5 ics 
OX ) ayz | bie? ' 


ss 
(=) ~Cou E> oi aa 


Some thermodynamic relations of rigid hygroscopic gels 587 


where £,, is Young’s Modulus in the x direction, o,,, is the ratio of the contraction 
in the x direction to the extension in the y direction for a stress in the y direction. 
Since a positive value of X denotes a pressure, whereas in defining £, tensions 


are regarded as positive, the signs in equations (4) are interchanged with their 
usual values. 


Substituting from (3) and (4) in (2), 
dx\ _ (ox x ysl (oz 1 *) } 
ee (Fa). ie re (55) cr si y (i XYZ 
ee 
BOWE NEN OA Jere 2 NOA) aus 
eer fox 1 (dy 
Ags 24 (5 Feet ts x) ae 5) 


or, rearranging and multiplying by yz6X, 


Ox Ox 
»(5a) Oke =ysdX(1+ 2 Ey et aie 


a 
Ties 1 oy Soh Oay gh N55) E 
+enbx(F22-7-x)(34) ayox( B45 xX aA) 
similarly 


oy is 1 ) =) 
aa =z ire & 
an() OY y, sy (74 Zs E, SEN oe 


sexd¥(14 5 + eet 2)\() ae (ee eee (Fee 
and 
© (Fa eas? Y 2) 6a a 
8 aoe ER ) (Pi) ove tay s(t : xt ae Baal. 
Adding the three equations gives 
ya (5) ove OX + 3x (3) bY ty (3) ve 
= (3) fax E, (SY —0,,8X — 42,52) — 702 —6,,0X — oy8Y)} 
Cl Se ee 
(5) ve DY ose. Ae sae (6) 


But 6X =yzbX+Xyz (= + *2) , where 5z and dy are due to arbitary changes 
= iy 


in stress at constant A, and 
bz 1 Fy Cr ? 
sh eee eye 2 SV = (02 —o,,0Y —0,,04)/E,, +-+0-5> (7) 
PSA hee E ee 


since o,,/E,=%y,/Ez ete. (Love, 1920, p. 105). 
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8X = ya {SX —X(8YV — omy SX — 0492) [Ey —X(8Z — oy, 8Y — on¢5X)/E,}; 


Varden (8) 
and similarly for 5Y and 6Z, so that if X=Y=Z 
ss (0x/0A)xyz 6X + 2x(0y/0A)xyvz Op ese xy(02/0A)xyz bZ 
= (0x/0A)xyz ox ap (dy/0A)xyz oY SF (02/0A)xyz OZ5 thee (9) 


Hence when the initial stress is a pure hydrostatic pressure, (1) gives 


= xyz Ol hen 6Z 
C2) De : (Faye X +5 (34 XYZ + 2\a4 XYZ 


1 (dx a) : (3) } 
=V4— oe feel Wee OZ fa) enue 10 
are Gees C XYZ a x \0A) xyz ce 


When the stress is a pure hydrostatic pressure we can express the variations }\p 
in B at constant temperature and composition in terms of variation of linear 
dimensions at constant stress and independent variations in the stress components. }} 
When 8X =8Y=5Z=5P the relation can be obtained more simply from the}} 
ordinary thermodynamic relation for fluids (e.g. Guggenheim (1933), equationay 


(151) and (153)) by the relation 
(1/V)(@V/0A)p = (1/x)(0x/0A)p + (1/y)(Cy/0A)p + (1/2)(02/0A)p. 


When the initial stress involves pure shear stress the error involved is that of If 
neglecting (X — Y)z5y+(X—Z)yéz in comparison with 6X, which is small if} 
(X — Y) etc. are small compared with the moduli. 


§3. DERIVATION OF THE GENERAL RELATIONS IN TERMS 
OF GENERALIZED CO-ORDINATES 
For a rectangular block of gel of the type considered, the general equation 
is given by Fowler and Guggenheim (1939), p. 60, asd F= —SST—=UXbx+ psn, 
where F is the Helmholtz free-energy, and if we define G as F+ X Xx by analogyiy 
with the definition of Gibbs’ function G*, 


6G = — ST + xXxdX+ psn, 


hence 

S= —(0G/OT xyz, *=(0G/OX)yryz, etc., and w=(0G/n)pxyz; 
hence (0S/OX)nryz = —(0?G/OT OX) = —(0x/0T )nxyz 
and (Ou/OX)nrvz = (0x/On)rxyz. 


For variations at constant temperature and composition 
(OS) ur = (0S/0X)nryz5X + (0S/OY)nrzx SY + (AS/0Z)yrxy 5Z, 
and 


(04)nr = (0x/0n)rxyz 5X + (Ay/On)oxyzdY + (Oz/On)pxyz5Z. oe eee (12) 


* G is not identical with G since in the case of a pure hydrostatic pressure Y Xx=yz Xx }}| 
+2xVYy+xyZ2=3PV. 
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When the initial stress is a pure hydrostatic pressure these can be transformed 
by the relations obtained in the previous section: 


1 (ox 1 (oy 1 /dz 
( ) : t5(5P) wxes®X +5 (33) cp, 2¥ 42S), 52h aaa. (itd) 


1/ Ox 1 2 | 1 /dz 
Of np =V4—\ = OX ES Vira Fl) Ban LO Lae eee ne 
ae A tle sere ® a5 ene p er) 


but (1/x)(@x/0T),xyz is the coefficient of thermal expansion in the x direction, 
a, and (1/x)(0x/dn)rxyz is the isothermal swelling per mole in the x direction 
at constant pressure =s,,; hence 


(8S)er= —V(og3X+0,S¥+0,8Z) eee, (11.2) 
and (Sar =V(s,3X+5,8V¥+5,8Z). see (12.2) 


Equation (12.2) can be used to calculate the effect of shear stress on vapour 
pressure, and the variation of vapour pressure of a hygroscopic gel initially 
subject to a pure hydrostatic pressure (or when shear stresses are small compared 
with the rigidity) with any variation in its state. Equation (11.2) can be used 
to calculate the contribution of the energy of expansion to the heat of absorption 
of the moisture by the gel. These will be considered separately. 


wan PEPRECT OF CHANGES IN STATE OF A-.GEL- INITIALLY SUBJECT 
TOPPA URE Sry ORO SWAT CePRES SURE ONS LE EQUILIBRIUM 
VAPOUR PRESSURE 


Equation (12.2) gives the variation in w at constant temperature and compos- 
sition with changes in the stress components; the net variation when temperature 
and composition may also change is therefore given by 

du =(0u/0n)prdn+ (Ou/0T)pndT + V(s,5X +5,5Y +5,6Z), ...... (3) 
since constant X, Y and Z when X = Y =Z is identical with constant pressure. 
If jy is the chemical potential of the vapour in equilibrium with the absorbed 


component in the gel, wy =p at all times. 
The only parameters on which wy depends are the vapour pressure p and the 
temperature 7. Equilibrium can only be maintained if T is the same for both 


phases. 
J. Ou =Sup =(Ony/Op)rdpt+(Our/OT)pST. ss ae (14) 


Combining (13) and (14), 
(Ou/On)pr hn + (On/0T pn dT + V(s,5X +5,5Y +5,82) = (Ony/Op)rdp + (Quy /OT )pdT. 
We have (Guggenheim, 1933, p. 35) 
(0u./0T pn = —(0S/On)pr, (On7/OP)r =2, (Quy/0T), = —Sy, 
where v= molar volume of the vapour. 


“. (Qu/8n)prdn—(0S/an)pr8T + V(s,8X +5,5Y +5,8Z) =vdp — SydT. 


ony (15) 
Considering first isothermal changes at constant stress: 
(O./0n) prin =v op 
. (0u/0n)pr = v(0p/0n)pr. 
39 
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The general equation for isothermal changes is 
V(s,5X +5,8Y +s,5Z)=vdp—v(0p/On)xyz7on. —... «ss (16) 


This is identical with equation (13) in the paper by Barkas (1945) on swelling}! 
stresses in rigid gels, but the derivation given here shows that it is only rigid}\lp 
if the initial stress is a pure hydrostatic pressure. ‘This can also be written 


V4 (sat sy ts.) (oe) (82 — Sy) (") 


SY 207 VAY. 0 
+(-9( 3 ) +650) 3 )=sop—0(Z) am 


thus separating the effect of changes in hydrostatic pressure and shear stresses. 

Equation (15) can also be used to determine the variation of vapour pressure} 
with temperature keeping the composition constant. This can be done even} 
if the absorbed liquid is in more than one form, as equilibrium must be maintained] 
between these forms, and, therefore, according to the phase rule, in discussing}! | 
equilibrium with external conditions the absorbed liquid can be treated as aj 
single component. In this case the only variation on the L.H.S. of equation (15) 
is in T’, whereas for the vapour phase the pressure must be allowed to alter tof 
maintain equilibrium. 


= (52) ,,97 =080—SyaF, 


5,2 (yea eye 
é On}prp \OT/ap p \OT np 


(if the vapour behaves as a pérfect gas) ; 


oS a] 
{Sy ‘ i = (symp). 


oS oH 
but re =| = — 
: De Le & Je a 


(Guggenheim, pp. 37 and 38), where H is the heat content of the gel ; 


0H 0 
Hy (pnt eer) =R TZ inp) ; 


nP 
but p=py, 


dH 3 | 
Tie (Fo RMF. Inp) tt ae (17) | 


. This may be combined with the relation obtained by Guggenheim (1933, p. 64),}) 
using a similar analysis, for the equilibrium vapour pressure of a pure liquid] 
maintained at constant hydrostatic pressure, which may be written 


rs) ai) 
Hy Hy=RT(snInpy) «ad deoe (ge . 


where Hy, is the molar total heat of the liquid and py the pressure of the vapour 
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in equilibrium with the pure liquid. Whence—since for a perfect gas Hy is 
independent of the pressure at constant temperature— 


0H 0 p 
Hi _ je =RP\ 55 a Oe ‘ 


The heat absorbed when 1 mole of liquid is absorbed by a large quantity of gel 
at constant temperature and pressure is given by the increase in total heat 


0 
(= AH pn 


(@AH|8n)pp =(8H/On)pp — Hy. 


where 


é O.4 p Wha) 
as SOS * Sa ar) (ain Pak i ee ae 
(; ai) ; RT (sn a = RT» Fn n) $i. adoiie (19) 


where / is the relative humidity of the vapour in equilibrium with the liquid 
absorbed by the gel. 

The above equation is identical with Kirchhoff’s relation for the heat of 
dilution of a dilute solution, but the above derivation shows that it is perfectly 
general and makes clear precisely what quantity is given by RT?(0/0T Inh)np. 
Previous derivations of the relation by Shorter (1924) and Hedges (1926) have 
either directly applied the relation for dilute solutions on general thermodynamic 
reasoning or have been based on the assumption that the volume swelling per 
mole was identical with the molar volume of the pure liquid, which does not 
hold for keratin at low regains. It should be noted that the heat of absorption Q 
as defined by Hedges and Shorter = —(0AH/0n)/M, where M is the molecular 
weight. 

Equation (17) has also been derived by earlier workers, notably Kruyt and 
Modderman (1930) and Brunauer (1943), who applied it to the heat of absorption 
of gases by charcoal, etc., and who give it as RT?(0/0T Inp),, =q, where q is variously 
defined, whereas in equation (17) the heat is accurately defined in terms of 
generally understood thermodynamic quantities. The relation of this to the 
observed heat of absorption in any experimental arrangement is a matter for 
thermodynamic analysis of the exact procedure employed, and will differ for 
different experimental arrangements. 

In the case of direct absorption by a gel from the pure liquid, the usual pro- 
cedure is to immerse a given quantity of the gel containing m moles of the absorbed 
liquid in a calorimeter containing the pure liquid and observe the heat evolved. 
Since in this case the pressure on the liquid and the gel remain constant through- 
out and the rise in temperature is negligibly small (although the quantity actually 
measured), the heat evolved is clearly, 


sat (OAH 
~ | ( ) dn, 
n On / pr 
where the integral covers the range from the initial content of absorbate to 


saturation. For an amount of gel containing 100M grams of the dry substance, 
where WM is the molecular weight of the absorbed liquid, m is the regain C 


39-2 
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(amount of absorbate per 100 g. of dry substance), and hence the amount of 
heat evolved per gram of the dry substance is 


~(1/100M) if (5) 


sat A 
=~ (1/1800) | ; (=), ,ac for watce) eee (20) 


PT. 


on 


In the case of direct absorption of the pure liquid by the gel, the heat evolved || 
is thus directly related to the quantity given by Kirchhoff’s relation, which has || 
been shown to be an exact relation between the heat absorbed when 1 mole of | 
absorbed substance is absorbed by a large quantity of gel material at constant 
temperature and pressure on the one hand, and the isosteric variation of equi- 
librium humidity at constant hydrostatic pressure on the gel, on the other hand. 


§5. EFFECT OF GENERAL EXPANSION OF GEL ON THE INTRINSIC 
ENERGY AND THE CONTRIBUTION OF THE ENERGY OF} 
EXPANSION TO THE TOTAL HEAT OF ABSORPTION 

Equation (11.2) expresses the effect on the entropy of varying the stress com-| 
ponents on a rigid gel initially under a pure hydrostatic pressure. This can 
be used to calculate the change in intrinsic energy 6U, when the gel is expanded} 
by the application of external forces by the same amount as when an additional 
mole of liquid is absorbed at constant temperature and pressure. Experiment | 

does not directly give the changes in stress required to do this, and equation (11.2) 

must be transformed by expressing the changes in stress in terms of changes in| 

the geometric parameters x, y, and z. This can be done by means of the ordina | 
elastic equations (for rhombic symmetry) : | 


one a6 éY 6Z : 

= phan ot Beg = ee 
Sy SX Va. Oo 

= y vena he ellen eaeeoe (21) 
bz Ox OY O72) 


gi Gen Bom Sie 
The negative sign appears on the L.H.S., as the conventional definition of £ 
assumes that tensions are positives, wiiereaes in the thermodynamic convention 
tensions are regarded as negative pressures. 

These three equations can be transformed to another set expressing the} 
stresses in terms of the resulting strains. Writing 4, for 'X/E,, this becomes 


4 
| 


— DA, = (1 — ogyoyz) 8x/" + (Czy + FxeFzy) BY/Y + (Coe + CxyFye) 92/2; 

— DA, = (1 — Gg) 8Y/Y + (Gye + SyaF ne) 82/2 + (Gin ye NON) oa) Meee (22) 
— DA, = (1 = 6 gySyz,) 83/8 + (Oge + zy Fyq) 8%/% + (Szy + FeaxF cy) BY/Y, 
where D is a determinant of which the value is 


D=1-—o,,,c,,,—6,,0 


awyl yx 2eP za — Fyz% zy ~ Fay y2 2a — Sy a %zyF xz: 
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A similar transformation for colour equations is given by J. Guild (1924). These 
values can be substituted in equation (12.2), giving 


(6S) nx = ( V/D)[E, a{(1 ais THe) dx/x a (Ge) te Trg) éy/y 
sl fees as Gi ye) 33/3} ag Eye, ((1 > Cue) oy/y zee } a EA 30S oatG }). 


If U is the intrinsic energy, 5|U = 75S — X Sx at constant composition ; or, when 
the initial stress is pure hydrostatic pressure, 


S8U=T3S— P8V; 
(8U)nr =(TV/D)[E pial... ty ao eee ay eae REN 


If we assume that the internal energy of differential expansion is the same whether 

the expansion is due to external forces or absorption of liquid,* the contribution 

of the expansion to the total heat or, more strictly, the total heat of absorption, 

when there is no heat of reaction between the absorbed liquid and the gel when’ 
absorption takes place at constant temperature and pressure, is given by 


Bord exnancion= (L074) Ee Ont mie: Fael Cyoe ioe « ct ce Ocean \] — PdV + PSV, 
since at constant pressure 6(PV)=P6V. Expressing this per mole, 
(OH /0n) expansion = (TV /D)[E,%q{(1 — GeyFy2)(1/%)(Ox/On)pr+ ...... | A erates ] 
= (LV {DYE 04 (b = gy Gys)S e+ (Cet CaP ey Sy to oes es RA oo (F 
Serie (24) 


writing s,, for (Ox/dn)/x, etc., in each of the terms in the equation. 
If there is symmetry about the 2 axis, 


ay YX Lx? 
Cy Oxz_) 
Oyz = Ox7z) 
and Sy=Sys Ch op = yy and £,=E,,. 


The expression for D then reduces to 


D= (1 =e Fru)(1 BO a 265532) 


and 
(0H, / ON) expansion = (T V/ D \[E a am O xeF2)S Pap (Cre fF Ono ¢)S x 
+ (Cig + Fn xF ne)Sat +E phat. sss. + EyaJ{(1 —0%.,)s, + 2(Oz¢+ Coe Fide Say l 
= [TV /a a8 orl = Ceca 20 eFax) [2E %94(1 + Orn) ae Oz( 1 te Le): 
a8 E,aA(1 ar On d PP Cae) s 25 205,(1 a Fre)Sah] 
a {[TV/a —~9g9— 20822) [2B (Sip ci @ p83) ty Ef Ql zy5 x we (1 ea O pap)Se > 
Jatin (25) 

but , EeGag = E Gee) 


* This assumption is the simplest one, and would appear to be justified by the fact that it gives heats 
of the correct order when this is the principal source of change in intrinsic energy on absorption of 


water. 
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therefore an alternative form is 
(0H/dn \excpaseton = {TV/( Si 26 nxn) LZ AeaS a 
ar Ops S7 ar O28) ate Bett im ee || ee Seen (26) 


For most hygroscopic gels s, is only about 3% of s,, and it is probable that the 
same is true of «, and «,, (as it certainly is for wood). In this case, the last term 
in the square brackets can be omitted and we have 


(0H/ON)expansion =2 TIVE ROS + Oy, (% 25, ai aS) }/(1 er cae ZG pT pave 


fy 


It should be emphasized that the above analysis assumes that the gel 1s homo- 
geneous in the limited sense that external stresses can produce the same dis- 
tortions at every point as absorption of the liquid. 


§6. APPLICATION TO THE STUDY OF THE ABSORPTION 
OF WATER BY TEXTILE FIBRES 

Urquhart and Williams (1924) determined the variation in the equilibrium 
vapour pressure of cotton with temperature and regain and give the results in 
the form of a series of iso-humidity curves. From these it is seen that at high 
relative humidities the equilibrium regain passes through a’ minimum as the 
temperature increases, the temperature corresponding to this minimum in- 
creasing as the humidity decreases. At these minima the relative humidity 
in equilibrium with the corresponding regain is a maximum, ie. (0/0TInh),p is 
positive for lower temperatures and negative for higher temperatures. Hence, 
from equation (19), (@AH/dn)pp for the particular high regain concerned is 
negative for lower temperatures and positive for high temperatures. 

A somewhat similar result was obtained by Speakman, Stott and Chang (1933) 
for wool. They found that the saturation regain and that for 97:-5°% R.H. 
showed a minimum at temperatures between 40 and 50° c. and also, as might 
be expected, that the swelling in water was a minimum at about the same tempera- 
ture. ‘The change in sign of the value of (@AH/én)pp with increase in temperature 
is most easily explained by the assumption that (@AH/dn)pr consists of two com- 
ponents of opposite sign, and that either the value of each changes with temperature 
or the contribution of each to the total value of (@AH/dn)pp changes. Such an 
assumption is made by Cassie (1945), who attributes the value of (@AH/dn)pr as 
observed by Hedges to two components—a negative component changing 
rapidly with regain, due to the heat of reaction between the gel and the water, 
and a positive component due to the heat required to expand the gel. Combining 
this assumption with the theory of Peirce (1929) and others, that the absorbed 
water consists of two portions, one of which has no heat of reaction, it is seen 
that (0AH/0n)pr for this portion is positive, while for the remainder it is negative. 
If equation (19) is applied separately to these two portions it is seen that the 
relative humidity in equilibrium with one portion decreases with increase in 
temperature, and for the other portion increases, if the proportions of the two 
components remain the same. Equilibrium can in fact only be maintained by | 
taking into account the first term on the L.H.S. of equation (15), which will have 
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two components, and adding the relation that the sum of the changes in 8n for 
the two portions is zero. When this is done, equation (19) is found to hold 
for the water as a whole, as deduced from the phase rule, but the proportion 
of the component with a positive value of (@AH/0n)pr will increase with rise in 
temperature, thus tending to make (QAH/dn)pz, for the water as a whole, become 
less negative as the temperature rises, in accordance with the experimental 
results quoted. 

At high regains and temperatures, when practically all the water is absorbed 
without heat of reaction, (@AH/0n)pp should approach the heat of expansion of the 
gel as given by equation (27). However, the comparison requires a knowledge 
of the transverse value of Young’s Modulus, £,,, as well as the three principal 
values of Poisson’s ratio for a material showing radial symmetry. Unfortunately 
it is impossible to obtain these constants for textile fibres, and an accurate com- 
parison can only be made by the use of a material like horn keratin, which is 
obtainable in bulk, and shows similar absorption properties to wool keratin. 
Until full experimental data for horn keratin are available, the best comparison 
that can be made is between (dAH/0n)pz, calculated from Urquhart and 
Williams’ isostere for cotton at 19-9% regain, which corresponds to about 98 % 
R.H., and for which full experimental values are given, and (0H/02)expansion 
calculated on certain assumptions for the value of the elastic constants, etc., of 
cotton. The assumptions are: 


(1) £, for cotton is about 104 bars, which is approximately the value for horn 
keratin and the gel material of wood at high regains (W.I.R.A. and 
Barkas, 1945, p. 55). 


(2) «, is that for beech across the grain, =6 x 10~° (Kaye and Laby). 
(3) 2V{s, +o,,(s, +5,%,/«,,)} = molar volume of water = 18 c.c. 

(4) One = Cn, =, =t (value of o,,, for keratin (W.I.R.A.)). 
Substituting in equation (27) for T=360° K.=87° c., 


10 =F 
(GEL On \estit = os x 24 x 10-8 cals/mole 
Sein S 


= 160 cal/mole. 


Urquhart and Williams’ data give the equilibrium vapour pressure of cotton 
at 19-9% regain when it is absorbing water at three temperatures, 70, 80, 
and 90° c. These can be converted into equilibrium humidities with the aid 
of saturation vapour pressures at these temperatures as given by the Smithsonian 
Meteorological Tables and the logarithm of h plotted against the reciprocal 
of the temperature. (@\AH/dn) at any temperature is then given by the slope 
of the curve. ‘The value found at 90° c. is 64 cal/mole, which is of the same 
order of magnitude as the heat of expansion obtained from elastic constants. 

Urquhart and Williams also give the absorption and desorption isosteres 
at 12% regain in the form of curves of the logarithm of the vapour pressure 
against the reciprocal of the temperature, obtaining a greater (negative) slope 
for desorption than absorption, showing that (0H/dn)pr is less positive for de- 
sorption than absorption. 
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This is consistent with the assumption of Barkas (1942) and Cassie (1945) 
that the hysteresis of the absorption isotherm is a reflection of the elastic hysteresis 
of the gel, as on this assumption the heat of expansion will be less for desorption || 
than absorption. 

The absorption isotherm gives the change in the thermodynamic potential 
or Gibbs’ function when one mole of the absorbate is transferred from the liquid 
in bulk to a large quantity of the gel at constant temperature and pressure. The 
general expression for the change in this function is given as 


AG=AH-TAS; 


differentiating with respect to gives the change in chemical potential under the 
conditions mentioned above, i.e., 


Aw =(@AG/On)pp = (OAH /On)pr — T(OAS/On)pr. 


The study of the isotherm is thus resolved into a study of the two components 
of the change in Gibbs’ function. The analysis of the entropy term involves 
statistical considerations and is outside the scope of this paper. ‘The energy |) 
term can be obtained experimentally either from the heat of wetting or from the 
isostere, and can therefore be analysed independently of the entropy term. 
The change in total heat AH can be written 


AH =AU+APV 


and the change in intrinsic energy can be split into two portions: AU, the change 
in intrinsic energy due to the reaction between the liquid and the gel, and AU, | 
the change due to the expansion of the gel. The APV term involves the change 
in the value of PV for the gel, and also for the liquid, so that, differentiating, 


oAH ce dAU, oAU, r) 
( On \z = On Ny +( an Ne +P(S eee Spee (28) 


where (OV /0n)pr refers to the gel and Vz is the molar volume of the liquid in bulk. 
The second term is that obtained from the elastic properties of the gel, and hence 
the second and third terms are equivalent to (@H/0n)expansion. If the absorbed 
liquid takes two forms with different energies of reaction, the differentiation 
takes place under conditions of equilibrium between the two forms. It is 
generally assumed that the liquid in one of these forms is localized, and in 
the other form relatively mobile. If A is the number of moles of liquid in the |} 
mobile form and B the number in the localized (bound) form, 


dAU, dAU, aA dAU, OB 
( an je 5 ( aA are( 35) es ( aB eled iter ae 
(9 Ne @ U,\ (aA @AU,\ (eB oH 
On PT a 0A Nae (¥ ) u ( oB ) seal eon iu ho. 3 fVa 
ee (30) 


At low contents of the absorbate it is generally assumed that all the liquid is 
absorbed in the bound form, i.e. when n =(; 


0A OB 
Gye Smee (F) 1 ‘ 
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ar) 7 i. =| ia EY. 31 
on PP n=0 os OB APT a0 on expansion oe Te secede ( ) 


Hence the value of (@AU,/0B) can be determined at low contents. If it is further 
assumed, as has been done by Cassie (1945), following Pierce (1929), that (@QAU,/0A) 
is zero and that (@AU,/0B) is a constant, substitution of the value of the latter 
in equation (30) enables 0B/dn to be found for all values of x. The values of B so 
obtained can be used in analysing the entropy term on statistical and other theories. 
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OPTICAL PROBLEMS OF THE ROTATING- 
PRISM CINEMATOGRAPH PROJECTOR 


By: J KUDAR 
Communicated by Professor A. F. C. Pollard. MS. received 4 March 1946 


ABSTRACT. 'The most characteristic optical aberration in the rotating-prism non- | 
intermittent cinematograph projector, the non-linearity of the image shift, can be eliminated 
by a suitable periodic variation of the direction of the illuminating beam. In addition | 
to a description of this unusual method of optical correction, other aberrations of the 
image formation by the rotating prism are discussed, such as the prismatic astigmatism 
and the tangential and sagittal prismatic coma. 


§1. INTRODUCTION 


device is used for producing a stationary projected image from the uniformly | 

moving film. This may be an arrangement of rotating or oscillating mirrors, 
or a rotating polygonal prism. A ray suffers parallel displacement by refraction | 
in a plane parallel plate, represented by opposite sides of the rotating prism | 
(cf. figure 1); a polygonal prism of suitable dimensions and speed of rotation can, | 
therefore, produce a stationary image. However, the parallel shift is not a linear 
function of the angle of incidence, and this gives rise to considerable optical 
aberrations. A theory of these was worked out by H. Dennis Taylor (1937) and 
applied especially for high-speed cinematographic cameras with rotating prism. 
The use of such cameras in practice is made possible by the circumstance that 
the light is cut off by a rotating shutter between two consecutive images in those 
positions of the prism where the angles of incidence, and thus the aberrations, 
would be too large. “This limitation of the angle of incidence is impracticable in | 
the case of projectors, because it would involve a frequency of interruption equal 
to the number of images per second. With the usual number of 24 images/sec., 
this would cause troublesome flicker* ; a dark period in the middle of the projection 
of an image, applied for doubling the flicker frequency, would cut out most of the 
time available for the projection of the image if its length were made comparable 
with that dark period between two images which would result from the exclusion 
of large angles. In addition, it would result in the loss of the optically most 
perfect interval. Increasing the number of sides of the polygon would reduce 
the greatest angle of incidence, but the natural limitation to the circumference of | 
the rotating prism prevents sufficiently good correction by this means. It seems 
more promising, therefore, to investigate optical possibilities for correcting or 
compensating the non-linearity of the image shift. 

Detailed consideration shows that, for practical reasons, it would be very 
difficult to attempt correction by means of rotating or oscillating optical elements 
between the film and the projection objective or in front of the objective. Fortun- 


|: cinematograph projectors with non-intermittent film movement, an optical 


* With modern intermittent projectors, the number of interruptions is two or three times the 
number of images. 
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ately, however, the non-linearity of the image shift can be eliminated by a very 
unusual method of correction, without introducing any constructional parts 
between the film and the objective or in front of the objective. This can be 
understood by means of figures 1 and 2. Figure 1 shows small pencils of rays 
coming from the film F parallel to the optic axis A~A of the projection objective O 
and deviated by refraction in the prism (of which p and p’ are opposite parallel 
sides). 
§2. DETAILED THEORY 


Let D be the distance between the sides p and p’, » the refractive index of the 
prism, x the angle of incidence which is equal to the angle between the objective 
axis and the prism face normal and y the angle of refraction. he parallel shift 
h of the rays due to the refraction is 


D sin (« —y) COS x n—1 n+1 1 
h=———_= = eet se NaS — tee) oe 
a Dsin( 1 Visa 5) D “+ (55 :)e +f, 


where the terms in x° and higher powers of x can be neglected for practical 
purposes. Equation (1) shows that the non-linear aberration of the parallel 
shift, represented by the term with x, is positive. A small pencil from any point 
M of the film produces an astigmatic image M’—-M”’ which is nearer to the optic 
axis by the distance (1) than the original point M. 

Let now figure 2 represent another case in which narrow parallel pencils of the 
illuminating rays remain perpendicular, not to the film, but to the face p of the 
rotating prism. Let again x be the angle between the axis of the projection and 
the prism face normal. In this case there is no parallel shift, but, seen from the 
objective O, a point M of the film appears as an anastigmatic image M’ produced by 
the refraction at the prism. For pencils of small aperture the distance between 
M and M’ is 


Relatively to the point M, the image M’ is nearer to the objective axis by the 
distance 


=) 3 
hv =MM’ sinv=D°—=(«—-F +... ), Prsker (3) 


Equation (3) represents the compensation of the movement of the film in the case 
shown by figure 2, while equation (1) represents the same for the case of figure 1. 
Both cases involve non-linear aberrations, as shown by the presence of the terms 
with x?; the essential difference, however, is that the sign of the aberration is 
different. 

Figures 1 and 2 represent extreme cases in which the movement of the film is 
over-compensated or under-compensated. By a suitable choice of conditions, 
therefore, an intermediate case can be obtained in which the aberration term in x? 
disappears. The condition for this is provided by the following rule: the ratio 
of the angle of incidence of the principal ray illuminating a point of the film to 
the angle between the prism face normal and the objective axis must be constant 
during the rotation of the prism. This ratio is 1 in the case of figure 1, and 0 in the 
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case of figure 2; the non-linear aberration disappears at an intermediate value 
which depends on the refractive index of the prism. . | 

This optical condition for the elimination of the non-linear aberration of the | 
image shift can be fulfilled, for example, by introducing between the film and the 


Muminati 
Beam “9 


Figure 4. 


light source a rotating optical device that influences the direction of the illuminating || 

beam without taking part directly in the formation of the projected image by the 
objective. Figures 3 and 4 show a schematic arrangement of the projector witha 
suitable beam-deflecting device satisfying the optical condition just mentioned. |/ 
In figure 3, O is the Projection objective, P the rotating polygonal prism which || 
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produces a stationary image from the uniformly moving film F, R a rotating hollow 
polygonal prism which is the moving component of the beam-deflecting device. 
Inside R is the stationary component S with a cylindrical surface c concentric with 
“the cylindrical inner surface c’ of the rotating component R. The uniform 
movement of the film F and the uniform rotation of P and R (as shown by the 
arrows) are mechanically synchronized and the consecutive pictures 1, 2, 3, 4 on 
the film correspond to the prism faces 1, 2, 3,4 of Pand R. The illuminating 
beam crosses the plane face a of the fixed prism S, the diagonal face b of which is 
a reflecting plane surface. The beam thus reflected at right angles crosses the 
concentric surfaces c and c’ and leaves the moving faces (2, 3 in figure 3) of the 
rotating polygonal prism R for traversing the film F, the rotating prism P and the 
objective O. The illuminating beam is more or less deflected according to the 
varying angle of refraction at the prism faces of the rotating polygon R. It is 
easily seen that the ratio of the angle of incidence of the principal ray traversing 
a point of the film to the angle between the prism face normal of P and the objective 
axis is approximately constant and less than 1. If this ratio has a suitable value, 
the non-linearity of the image shift is eliminated. 


§3. FURTHER ABERRATIONS 


In what follows, other aberrations produced by the rotating prism P will be 
discussed which, although less striking than the non-linearity of the image shift, 
must nevertheless be considered in any practical solution of the problem. These 
aberrations can be classified in two groups: (a) those which are essentially inde- 
pendent of the angular aperture of the illuminating beam, and (6) those which are 
proportional to the first and second powers of this angular aperture. 

(a) The non-linearity of the image shift belongs to those aberrations which are 
essentially independent of the aperture; another is the chromatic variation of the 
image shift, which does not exceed the practically tolerable limit. A further 
aberration may arise as a consequence of film shrinkage, in which case the length 
of the film per image would become less than that on which the design of the prism 
isbased. An optical compensation of film shrinkage by an automatically respond-. 
ing device is of course desirable when using considerably shrunk films in non- 
intermittent projectors; however, it would be difficult to discuss this question 
before treating the optical problems of the image formation in the standard case of 
non-shrunk film. 

(6) The illuminating rays crossing the film have wide apertures in practical 
cases; consequently, the aberrations depending on the aperture need detailed 
investigation. ‘These aberrations, which determine the quality of the projection 
if linearity of the image shift is achieved, will be treated in the following discussion 
of the image formation by wide-angle pencils starting from any point of the film. 
Lens calculations, based on the approximation sin « =« — «3/6, give account of 
the following five aberrations (apart from the chromatic error): spherical aber- 
ration, coma, astigmatism, curvature of field and distortion. It is possible to 
design the projection objective and the prism behind it, in the position in which the 
optic axis is perpendicular to the prism faces, together as a well-corrected lens 
system. ‘Then the rotation of the prism produces periodic aberrations which 
_ have no axial symmetry and are to be related to the tangential and sagittal sections 
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| 
| 


of the rotating prism. These prismatic aberrations are obviously analogous to 
the five lens aberrations above mentioned. For instance, the non-linearity of the || 
image shift is an effect corresponding to the distortion in lens optics. There is a}| 
prismatic astigmatism, while the field curvature of lens optics corresponds here to | 
the varying positions of the plane of the astigmatic image. Finally, the periodically 
varying prismatic coma is also a very interesting and practically important aber- | 
ration. The analogy between these periodic aberrations due to the rotating 
prism and the axially symmetrical aberrations of lens optics becomes immediately 
obvious by considering the variation of the former during the rotation of the 
prism. ‘Thus the axially symmetrical aberrations due to the prism, when its two | 
faces are just perpendicular to the optic axis, will change continuously to the 
tangential and sagittal prismatic aberrations. 


| 
§4. TANGENTIAL RAY-TRACE 


We begin by tracing tangential rays through the prism. Thus figure 5 is 
crossing the prism axis at right angles. Let M bea point on the film; p, p’, two | 
parallel faces of the prism; x, the angle between the objective axis and the prism | 
face normal; e-e, the second principal plane of the objective; MO, a ray through | 
M parallel to the objective axis, crossing the principal plane e—e at O after traversing | 
the prism faces p, p’. Let further B, and B_ be points of intersection of the 
plane e-e with two other rays from M; +7 and —7, the angle between these rays 
and MO; M’, virtual intersection of the rays MB, and MB_; M1, virtual 


e 


Figure 5. 


intersection of the rays MO and MB,; M’, virtual intersection of the rays MO 
and MB_; t-t, a plane through M’ parallel to e-e and to the film plane. 
We calculate first M’O+ and MO, i.e. the distances of these virtual image 


points from the principal plane. Let m be the distance of M from e-e. The 
parallel shift of the ray MO is 


Dei ees 
sin( 1 Wir Becca (si 
and the parallel shift of MB, 


sin (x+7 ERTS OSCE aga 5 
Dsin(x+ | Ve=amteeaik oe uh) 
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wherein D is again the diameter of the inscribed cylinder in the prism. Then 
the distance OB. is 


: cos x ; cos (x +7) t 
t D Ws ey er |i WA Oe sine ae) 
mtan7r+l sin x( a/(n2 —sin2 5) Dsin (x Se “fit /(n? —sin2 (x aN >| coor. 


Introducing ¢, =OM/,, we obtain 


_OA+ _ pam*(1 cos x )-pe eta cos (% +7) 


t — WF > oS SS ea EE : a SW ee eee 
PS tans tanz\ 1/(n? —sin? x) sin t a/(n? — sin? (x +7)) |’ 


which, breaking off the series of the trigonometric functions after the second 
power, results in 


T15: 2 mya 


n—1 n+1 a 3nt1 ntl 
2 


t. =m— D— 
me n 


Equation (8) reveals several aberrations : 


The term in 7?, being independent of the rotation of the prism, is simply the 
spherical aberration, and is the same as in the position when the prism faces are at 
right angles to the optic axis. ‘The term in x? represents the position of the image 
formed by an infinitesimal pencil of rays in the tangential section through the 
prism. Finally, the term in x7 is a measure of the varying prismatic coma. It is 
of practical importance to find the intersection M’ of the rays MB, and MB_. 
The distance of M’ from e-e is 


n—1 Dek ON eal 
sb area 1+(5.45-3 )e Eel SO0CO6 (9) 
and the distance of M’ from the film plane 
n—1 oft hi n+1 
ares [1+(5. 5h Sa |. BESO (10) 


The term in 7? is of no importance, because the spherical aberration of the prism 
(in the perpendicular position to the optic axis) can be corrected by the projection 
objective. 

In the image plane t-t the tangential coma produces the vertical aberration 


MM’ ,M’'_ Sonal Ss al eta ae 
—7 + tant=D—— «5. =) wr tan tr D . 5 5 


§5. SAGITTAL RAY-TRACE 


The tracing of sagittal rays through the prism can be effected by reference 
to the perspective drawing figure 6. Taking M, a point on the film; MK, a ray 
parallel to the optic axis; MK’, another ray in a sagittal plane, i.e. in a plane 
parallel to the prism axis (and parallel to the objective axis, in the particular case 
of figure 6) ; and taking also K and K’, the points of incidence of these rays 
on the prism face p; V, the point of intersection of the perpendicular from M with 
the prism face, we get T for the point on this perpendicular for which MKLTK 
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and write x, the angle of incidence of the ray MK; «, the angle of incidence of the | 
ray MK’; o, the angle between MK and MK’; @, the angle between MK’ and || 
K’T; «, the angle between KT and K’T. We use the sign <¥ for a right angle. 


For the triangle KVK’ we have 
sin?% =sin? o+cos?o . sin? x. a IA eee (12) 


It may be noted that the refracted rays of MK and MK’—which are not shown 
in figure 6—have no vertical intersection, unless o is infinitely small. Con- | 
sequently, there is no Pythagorian equation for the corresponding angles of | 
refraction. 


Figure 6. 


We find for triangle KTK’ 


sing coso.tanx 


a ee enn eee (13) 
sin w cos w 
and for triangle MK’T 
sng cose. . i 
ain oe ae oe) Uh Re ee) re Pein aded (14) 
(13) and (14) give 


sin x sina. 15 
cos w = sin? eeeveece ( ) 


The ray MK’, after refraction at two parallel faces of the prism, has suffered 
a parallel shift of magnitude, 


é COS “ 
@=Dsino(1- Fea), oeecee( 16) 


in the plane of incidence. Produced backwards, the refracted ray intersects |} 
the film in a point whose distance from M is 


d/sin 3, oh tes alee 


since the film plane and the triangle KK’ are parallel. The displacement (17) | 
in the film plane has a vertical component 


dcosw/sind Pr os) 
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and a horizontal component 
ESI oy/Sin he 8 ee Nae (19) 
the prism axis being assumed horizontal). 


The vertical component (18) can be obtained from (15) and (16); if the 


power series for the trigonometrical functions are broken off after the second 
term, we have 


pcosw |. *. COS « ptz 1 n+1 1 n+1. 
sind EAN / (n? — sin? yh See 2n2 6 on cal 
Fears (20) 
The horizontal component (19) is given by (13) and (15) 
sin w COS % 
Sn p7D tang. coss(1- ene), re (21) 


Let us now consider in the plane MK’K of figure 6 another ray at an angle —o 
with MK, so that this new ray and the ray MK’ have a point of intersection the 
distance of which from the film plane is 


sinw 1 COS % n+1 n+1. 
iad tan 5=Deosx(1- Tera EP [14 (So ~ 5) 2n? 0] 


in (22), the term in o? is the spherical aberration of the prism, which does not 
shange during the rotation, while the term in x? gives the position of the image 
ormed by an infinitesimal pencil of rays in the sagittal section of the prism. 
The difference between the x?-terms of (10) and (22), 
n* — | 


n® 


D 


RS erA limes. 108 (23) 


s the distance between the tangential and the sagittal points of intersection, 
.e. the prismatic astigmatism. 

In (20) the term with x? is the well known non-linear aberration of the 
arallel shift, and the term in xo? is the sagittal coma. 


§6. CONCLUSION 

These aberrations have been calculated for the special example where the 
rincipal ray illuminating any point of the film is parallel to the objective axis. 
“he formulae are the same for any other direction of the rays; in the general 
ase, however, x does not mean the angle between the prism face normal and the 
bjective axis, but the angle of incidence of the principal ray from any point 
n the film. 7 and o are always the tangential and sagittal angular semi- 
pertures of the wide-angle pencil illuminating a point of the film. In the 
eneral case, these aberrations as vector quantities have to be related to the 
irection of the principal ray in question, not to the axis of the objective. 

A more detailed discussion of how the results given in the present paper can 
e used for the optical correction of the image formation in the rotating prism 
rojector will be published later. 
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ABSTRACT. A rotational analysis has been done of the 0,1, 0,2, 0,3, 0,4, and 1,5 bands 
of the near ultra-violet band-system (v, & 35030 cm-+) of SiS. The transition responsible 
for the production of the band-system is of the type p11 —x1l. The rotational constants 
have been derived from two sets of experimental observations : (i) of the bands in absorp- }| 
tion, by G. Herzberg, N. Olson and F. J. Wilson (University of Saskatchewan), and (11) of || 
the bands in emission, by the author. The more important constants are given by 


B,” =0-3036 —0-0014, (v”-+4), for the ground state, and 
B,,/ =0-2664—0-0020; (v’+-3) for the upper state. 


A-type doubling in the upper J] state is negligibly small, at least up to J=100. The values 
of the equilibrium internuclear distances are r,“=1-929 a. and r,’=2:059, a. A brief 
discussion of the results is given. 


§1. INTRODUCTION 
T: E production and vibrational analysis of a near ultra-violet band-system of 


silicon monosulphide (v,> 35030 cm71) were described some time ago | 
by Barrow and Jevons (1938). Since then, a preliminary account has been 

given of the rotational analysis of the 0,1, 0,2 and 0,3 bands of this system (Barrow, 
1944), which was shown to be due to a transition of the type 1II—1X. The} 
present paper gives the results of a more detailed analysis of these three bandsand ! 
of the 0,4 and 1,5 bands of the same system. 
Two distinct sets of experimental material have been used. When the pre-}} 
liminary analysis of the 0,1, 0,2 and 0,3 bands from my own plates, taken in 
emission, had been completed, I learned from Professor G. Herzberg that} 
Mr. N. Olson and Mr. F. J. Wilson, working in his laboratory at Saskatchewan on 
plates of the bands in absorption taken by Olson, had already carried out analyses}, 
of the 0,1, 0,2, 0,3, 0,4 and 0,5 bands. Their results proved to be in close agree- 
ment with mine. Professor Herzberg then most kindly sent me the detailed{ 
measurements and analyses: the results to be given are based, where possible, |) 
on the averaged measurements of the two independent sets of investigations. | 


§2. EXPERIMENTAL METHODS 


Olson photographed the bands in absorption. SiS, prepared by the method} 
described by Barrow and Jevons (1938), was sublimed into an electrically-heated i 
‘T-shaped quartz tube. ‘This consisted of an absorption tube, 15-20 cm. lone ! 
and 2:5 cm. in internal diameter, fitted with a central 'T-piece which acted as j 
reservoir for SiS and through which the whole system could be evacuated and 


f 


, 
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ultimately sealed. Independent heating of the reservoir and of the main tube 
minimized deposition of SiS on the windows of the latter. The tubes were 
operated at about 900°c. A hydrogen lamp was used as source of continuous 
radiation, and the bands were photographed in a second order of a 6-m. concave 
grating spectrograph on Eastman “ 40” plates. 

The experimental arrangements used in developing the SiS system in emission 
have already been described (Barrow and Jevons, 1938). The region 2800- 
3100 a. was photographed in a fourth order of the 10-ft. concave grating at Imperial 
College, London. ‘The grating itself largely absorbed the overlapping fifth order: 
the overlapping part of the third order was absorbed in a 2-cm. column of an 
aqueous solution of cobalt and nickel sulphates containing 4:15 g. CoSO,.7H,O 
and 14:5g. NiSO,.H,O dissolved in 100c.c. water (Bowen, 1942). Ilford 
“Zenith” plates were used, and exposure times were of the order of two hours. 


§3. ROTATIONAL ANALYSIS 


The 0,1, 0,2, and 0,3 bands were first selected for analysis as their structure 
is little overlaid by that of other bands. However, the 0,2 band, which appeared 
to be particularly simple, proved to suffer from considerable overlapping of the 
branches, and the 0,1 and 0,3 bands were therefore examined first. In no bands 
were the lines resolved close to the heads, but further away from the heads 
three branches were apparent, one of which—the Q branch—was considerably 
stronger than either of the others. 

By a trial and error study of the first differences, R(J)— Q(J) = OQ(J + 1)— 
P(J +1), it turned out to be an easy matter, first to identify the P and R branches, 
then to ascertain the absolute J-numbering, and finally, using a rough value of D, 
obtained from the relation D= —4B?/w, to get fairly accurate values of the con- 
stants B for the 0,1, 0,3 and 0,2 bands. 

At this stage I received the data of Olson on the 0,2 and 0,3 bands and of Wilson 
on the 0,1, 0,4, and 1,5 bands. Their assignments of the lines of the 0,1, 0,2 and 
0,3 bands agreed with mine in detail, and, after taking account of a small systematic 
difference in the wave-numbers of the lines (my values were on an average 0-35 
cem;+ lower), means of the wave-numbers of the lines of these three bands were 
taken. Only Wilson measured the 0,4 and 1,5 bands. ‘The wave-numbers are 
displayed in table 1. 

The procedure adopted for the finai evaluation of the constants was as follows. 
The constants which could be determined with the highest absolute accuracy 
were obviously those for the upper state with v’ = 0, for which there were measure- 
ments of four bands. Either A,F’(J) or A,F’(J) differences could have been 
used: eventually the former, 


A, F'(J)=RVJ)- OV) OJ +1)- PU +1), 


were chosen after it had been demonstrated that, within the limits of error of the 
measurements, the differences R(J) — O(J) were equal to QJ +1) — PJ +1), up 
to the highest J-values (J~100). Thus the A-type doubling of the upper "II 
state is very small—as would be expected by analogy with the results of the 
analysis of bands of the similar system of SiO (Saper, 1932). 
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It is a familiar result that 
A, F(J) 
(J +1) 


where B, D and F are the constants in the expression for the rotational energy 
levels, 


=2B+4D(J +1)? +6F(J+1)4+..., 9 vs... (1) 


LEON ISIC IN OEE Opa ed eC Maa ee ee ae (2) 


and generally, B> D> F. 
APY) 


Thus, a plot of against (J+1)? should—at least for moderate 


Lal 

values of J—yield a straight line of intercept 2B and of slope 4D. A selection of 

AF y'(J) 
J+1 


A,Fy'(J) values for v'=0 is given in table 2, and a plot of against 


2000 4000 6000 8000 70000 12,000 
(ile = 
: AF UT : 
Figure 1. SAD) as a function of (.J+1)?. 


J+1 
The two points © represent the mean values for J/=45-76, 77-111: their coordinates are 
(0°52626, 3868) and (0°52027, 9127), giving 2By’ =0°5306¢, 4Do’ =1°14 x 10-8. 


(J +1)? is shown in figure 1. From the latter it appeared safe—in the first 
instance—to neglect terms in /’, and the best straight line was therefore drawn 
through all the points with J>45, using the centre-of-gravity method. This 
gave By’ =0-2653, and Dy’ = —2-85 x 10-7 cm=1. 

It was considered that the measurements of the Q-branch lines were more 
reliable than those of either P- or R-branch lines, and the lower-state constants 
for the 0,1, 0,2, 0,3, and 0,4 bands were therefore determined from the Q line in 
the following way. ‘The wave-number of any such line is given by 


Oerenthoah) Ill) = Do Its Ply toc san (3) 


where v, is the band-origin. Provisional values of B” and of D’’, obtained from 
A,F” (J) and A,F’’(J) differences, were inserted in (3), and preliminary values of 
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which we may write 7, calculated for each line, thus 
i) = O(J) -- (B’— B”) .J(J +1)-(D" DD"). Jal ve...(3a) | 
where the horizontal bars indicate provisional values. In hignen oe Vo values i 
were not quite constant with J, since the (B’ — (B’— B”) and (D’ — D’) values were not || 
quite correct. If the differences between the true and the provisional values are) 
indicated by A’s—thus, Avg = vg — ¥9, — then 
Avy = —A(B'— BY”). J(J+1)—A(D’ =D") .J*{F +1)? —... -«-- 2. (3b) 
The 7, values were then plotted against J(J +1): the intercept gave vp, and the) 
slope at low values of J(J +1) the correction —A(B’—B”). New values of v9 


Vo, 


Table 2. Some values of A,Fo'(J) 


| 
RG) OO) QJ +1)—P(J+1) Mean 


d/ obs. Calc. 
0,1 Oo 0,3 0,4 0,1 0,2 0,3 0,4 


40 | 21-48* 21-56% 21:°77* 21-46% | 21-58% 21-57* 21°35*. 21:80* | 2t-S7 ioe 
41 1-97® 22-05% -2-45* - 1°96") 198% 2:45") 2-00Mee 2:35 eal pees Pinas 
42 2°57* 2-42¥° 29-01%) 2-587) "261s 2-7/5" 27 er OS Dols 
43 3-20® 2:87* 3:47* 3°14%| 3:14* 3:-47® 93-27 3:31") 3:24 DoS 
44 |. 3-95* 3-40% 3:98* 3°75*| 3:68* 3:93* 3-61 3:30") mai 3°78 
45 4:44* 4:05* 4:64* 4-54*] 4:16* 4:86* 4:04  4:24¥*| 4-33 4-30 


65 | 34:62 34:65* 34-71 34:74 | 34-61 34-68* 34:58 34:50* | 34:64 | 34-69 
66 DS Ons 2S  eo-4 ORL, B83 Sey sil Oakes’ || Dow S2i 
67 S00) 0)7 4 OZ -OS BY Se A MO” || Sez! So? 
68 6:33 6-3) AO 22 OS 6:29 0-557 ae Oc SOD ecm On 6°25 
69 6:76 6:90" (6:84 = =6:66 6:52" 55 6:37.07 / an" 7 Ona mOoO 6°76 
70 TEES A PSDP 6:99 a7 40% cS Seed aa ee) UY 


90 | 47-48* 47-20% — — AT-41 47-83% — = 47-47 | 47-42 
91 3:09 = 57/2O1 — 7:94 8:26®  — — oi, doz 
92 8:60* 8-21* — — 8-41 8-65 — —= |) ee) 8-42 
8 eilik — hey) — — 8:88 9:05 == — 8-99 8-92 
94 Dao Doss — — 9-81 9-25 = == 9-49 9-42 
hy || SORE? Os: — — 50:00 9-90 = —_— DY 9-92, 


* Asterisks indicate values involving the wave-numbers of unresolved lines. 


were then calculated, with the corrected (B’—B’’) terms, but with the same, 
(D’—D’’) terms. Plots of Avg from these 7, values against J?(J +1)? enabled | 
A(D’ — D’’) to be evaluated. | 
When this had been done, it was found that the values of vp calculated from) 
(3) with the now corrected values of (B’—B’’) and of (D’—D”) still showed some! 
systematic departure from constancy in the region of the highest J values (90-110). 
The magnitude of terms in (F’—F”).J3(J+1)? was therefore considered.) 
F’ and F” were calculated from the relation | 


| 


' 
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(see, for example, Jevons, 1932), The results were 
F’ = —1-7x 10-8 cm?z1, 
F" 210 ema, 
so that, approximately, 
(fF — Ff") = —-1-7 x 10-8, 


The addition of the terms —1-7 x 10-8 . J3(J +1)? to (3) was found to lead to 
values of vy which were then sensibly constant over the whole range of J values 


. for all these four bands. 


The constants derived in this way were finally checked by calculating the 
wave numbers of all the measured lines of the four bands, by first evaluating the 
Q branches and then using the calculated A,Fy'(J) values to get the P- and R- 


. branch lines. 


The lines of the 1,5 band were so often superimposed on other lines of the same 
band that accurate evaluation of the rotational constants from either A,F(J) or 


_A,F(J) differences was found to be impossible. Accordingly, values of B,’’ 
_ and of D,’’ were obtained by linear extrapolation, and the quantities (B,’—B,”’) 
_ and (D,’—D,"’) determined from the Q-branch lines in the way described above. 


§4. RESULTS AND DISCUSSION 


The vibrational and rotational constants are collected together in tables 3 
and 4: some comments follow. 


(i) The structure of the bands makes it quite certain that the electronic 
transition is of the type 1I]—1X. ‘The ground state of SiS is most probably 
1y'+ (see, for example, Herzberg, 1939). ‘The system is therefore exactly 
analogous to the “4th Positive” atI]—x1X* system of carbon monoxide, and 
to the similar system of SiO, as was, indeed, suspected before the rotational 
analysis. 

(ii) The energies and products of dissociation of SiS in these two states will 
be considered in a later paper on another ultra-violet system of SiS (see also 
Vago and Barrow, 1946). However, it may be noted that the fact that the p1II 
state is a singlet makes it most probable that it, like the ground state, dissociates 
into the ground state atoms, Si(?P)+ S(3P). 

(iii) The new values of w,’’ and of «,w,'’, determined from the band-origins, 
are very close to the values obtained from the original band-head measurements. 
As the analysis shows, this agreement arises from the circumstances that the 
band-heads are not very far from the origins (~2 cm-*), and that the separation, 
head-origin, changes only rather slowly from band to band. 

(iv) It is to be expected that SiS should have a similar electronic configuration 
in its ground state to the isosteric molecule P,, which also has 30 extra-nuclear 
electrons. As with other sets of such molecules, this may account for the rather 
striking similarity in the ground-state constants (table 5). However, no upper state 
of P, corresponding to the D 'II level of SiS is known. 

(v) The bond distance in the ground state of SiS, about 1-93 a., is some 
0-22 a. shorter than the uni-covalent distance (Schomaker and Stevenson, 1940). 


The bond in SiS is, of course, expected to be stronger, and thus shorter, than a 
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single covalent bond, and the observed degree of shortening is in reasonable: 
accord with the figures of 0-25 and 0-15 a. for the ground states of the similar 


molecules SiO and GeO. | 
Table 3. Constants of the p1II and x1X+ states of SiS 
State : x 1+ p 4 
v 0) cm71 35028-8 cm>! H 
: 749-65 
(749°5) H 512-0 H 
Xp —2:575 
(—2°:56) H —2-38 H 
Vee — —0:045 H 
B, 0-2969, — 
Bs 0-2984, = 
B, 0:2999, — 
Ba 00-3014, 0:2632, 
Bs 0-3028, 0:26533 
3 00-3036, 0-2663. 
ay —0-0014, —0-0020; 
D, —1:9,x 10-7 — 
DD: —1:9,x10-" es 
iD —-1:9,« 10-7 a= 
‘De —2:-0)x 10-7 == 2-95 10 a2 
iDY. = —2)-85 x 10-7 
ier — —1-7 x 10=8 
ie 1-931, A. 2:063, A. 
We 1-929, A. 2-059; A. 
: 4-1065 F 
Notes : (1) The values of r were obtained from r= Tage A., where pu is the reduced. 
pB 


mass in physical atomic weight units. js was calculated for 785i 38S, 


with 


(2) H signifies a constant derived from band-head data. 


Table 5. 

State 

SiS x lyr 
P, . x 1 


R. F. Barrow 


| 


Si= 27-987 and S=31-982. 


Table 4. Positions of band-origins 


Band Origin 
sit 34163-8, cms 
0,2 3424:54 
0.3 2690°3, 
0,4 1961-2, 
5 1744-6, 
Comparisons of the constants of SiS with those of 


the isosteric molecule P,* 


Dissociation 
We Force constant energy Be 72 
749-65 4:94 x 105 dyne/cm. 6:6 ev. 0:3036 1-929 


780°43 5-56 5-03 0°3031 1-895 |i] 


* Herzberg, Herzberg and Milne (1940). 
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VELOCITY OF SOUND IN MIXTURES OF 
ARGON, HELIUM AND HYDROGEN AT 
LOW TEMPERATURES 


By A. VAN ITTERBEEK anp W. VAN DONINCK, 


Physical Laboratory, University of Louvain, Belgium 
Communicated by Dr. E. G. Richardson, 26 March 1946 


ABSTRACT. Measurements are described of the velocity of ultrasonics in mixtures 
of argon, helium and hydrogen at a frequency of 524 kc./sec. as a function of pressure 
and at fixed temperatures near 90° and 20° absolute. From these data, information about the 
interaction between the gases in the mixtures is deduced, using the equation of state of 
Kamerlingh Onnes. 
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OME years ago one of us, with W. H. Keesom, developed an indirect method 
S of determining the second virial coefficient of gases at low temperatures 

(v. Itterbeek and Keesom, 1930). ‘This method is based on measurements 
of the velocity of sound as a function of pressure at different fixed temperatures. 
Since then this method has been employed in our laboratory using ultrasonics, 
and has been applied to the following gases: A, O, Hy, D,, CO (v. Itterbeek et al., 
1938, 1943, 1944). In the last research we extended this method to gas mixtures. 
In this paper we report on measurements carried out on A-He, A—-H, and H,—He 
respectively at liquid-oxygen and liquid-hydrogen temperatures. From these 
measurements it is possible to obtain some information about the interaction 
between the gases concerned. 

§2. METHOD 


The velocity of sound W can be expressed by means of the equation 
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Taking into account the equation of state of Kamerlingh Onnes, 


B 
po=RT(1 =n =) ceo eee (2) | 
equation (1) can be transformed into 

S | 
W= WA1+ R- eT) ae ty (3)} 

W, being the velocity of sound in the gas mixture in the ideal state (p =0). 
Liye re) Ree | 
dere Cy X+Cy(1—x) My. 4+ M(1—x) (4)} 


(M,, M, are the molecular masses, R the gas constant, Cy, Cp’, Cy and ¢,’ the speci 
heats for the ideal state). 
Further, 
S=B T dB fh d*B 
=X dT Dede re 
A is a numerical factor given by 
— #.ot+(1l—#) ey 
Mi we laa\M, 2 ee (6) 
For the case of He-A and He-H, respectively at liquid-oxygen and a 
hydrogen temperatures, we have A= 3/2. 


Writing 
B=B,.28+.84.(1=0)+ Bax) re Me) 
together with equation (5), we obtain 
S=S,.¢2+.cry(1—n)4Sy. (1x) ip (8) 
with 
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For the case where the two components do not possess the same specilf 
heats (A-H,), A depends on concentration (6). | Equation (8) must be changy 
because the factors S, and S, are no longer independent of concentration. |} 
that case S, and S, will be indicated in what follows by the notations S,' and | 
and computed by means of equation (5). Discussing our measurements we cou}, 
however, observe that |S —S,’.x?—.S,’.(1—«)?]/x.(1—x) does not depend | | 


concentration, which proves that the expression 
I te i Ge d*p 
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is small.* Thus we can put «=f approximately, and obtain in this manner sos | 
indication of the interaction between the components at low temperatures. 


* Admitting for 8 the following expression: B=a+ 2 + im » we can observe that | 
dp, 7 ap | 
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becomes small when A is approximately equal to 3/2. 
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We proceed in the following way :—We measure the velocity of sound in the 
/ gas mixture as a function of pressure at different fixed temperatures. From these 
measurements, we compute by the method of least squares the different values 
of S corresponding to the different concentrations. S, and S, are in most cases 
known from previous measurements carried out on the components, or S,(S;’) 
_and_S;(S.’) are computed from the temperature-variation of B given by measure- 
ments on isotherms. 

The gas mixtures are prepared within a large vessel of capacity 4 or 101. 
respectively for the measurements at liquid-oxygen or liquid-hydrogen tem- 
| peratures. The concentration is at first determined by means of the partial 
pressures. ‘The mixture remains in this vessel for several hours. After this the 
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Figure 1. Acoustic interferometer. 


gas is let into the acoustic interferometer represented in figure 1. The velocity 
of sound is measured in the mixture at room temperature and from these measure- 
ments the concentration is finally computed from equation (4). At low tempera- 
tures the concentration is checked again by using the value of W, obtained from the 
measurements represented by equation (3). The values of x found at room 
temperature and low temperatures agreed generally very well. 

For the measurements on the pressure dependency of the velocity of sound, 
the pressure is lowered in the interferometer by letting the gas stream into an 
evacuated vessel of capacity 51. The gas pressure within the interferometer 
is read on a closed mercury manometer by means of a cathetometer. 

The frequency of the piezo-quartz used is 523-78 kc./sec. For the details 
of the electrical system, we refer to our previous publications. 
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§3. MEASUREMENTS ON THE GAS MIXTURES A-He 


The values obtained for the velocity of sound as a function of pressure} 
temperature and concentration are given graphically in figure 2. 
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Figure 2. Velocity of sound in A—He as function of pressure. 


Using these experimental data we have computed the factor S and also 
corresponding to equation (3). W,.,,. is the value of W, computed by meanjf 
of equation (4), « being the concentration observed at room temperature. It wi 
be seen that the agreement between experimental values and theoretical on 


is very good. 
Table 1 


The values S,-A and S,—He are taken respectively from the measurement} 
of van Itterbeek and van Paemel (1938) on the velocity of sound in pure argoi| 
and from the measurements of Nijhoff (1928) on the isotherms of helium gas 
These values are indicated in table 2. | 

The values of S of table 1 are drawn graphically in figure 3, the full lin} 
corresponding with the theoretical curves computed by means of the equation | 


S=S,.x2+ Bu. (1—x)+S,.(1—x)2. 
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Table 2 
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Figure 3. S in mixtures of A and He. Figure 4. Velocity of sound as function 


of pressure at 20° k, 


The interaction coefficient f is indicated in the third column of table 2. 


§4. MEASUREMENTS ON THE GAS MIXTURES He-H, AT THE 
BOILING POINT OF LIQUID HYDROGEN 


The values obtained for the velocity of sound as a function of pressure at 
20°-0 K. are indicated in figure 4. 

In table 3 (figure 5) we give the values of S and W, computed by means of the 
experimental values of figure 4. We see again that there is a good agreement 


between the values of Wo,,,,.. and Wo, 
Tables3 
Concentration (%) 
iP. 

(7K.) Wo exp 10° ..S Wo theor 
H, He 

20:3 0-0 100 S755 —3-39 375°8 
78-6 21:4 341-1 —2-44 341:0 
60°8 B92 Sul Cg —1-65 317-9 
41-4 58-6 297-6 —1-06 297-9 
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We found again that S as a function of the concentration satisfies equation (8), 
with B= —2:11.10°°. 
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Figure 5. Comparison of theoretical and observed velocity of sound. 


§5. MEASUREMENTS ON A-H,; MIXTURES 


The experimental results for the A-H, mixtures are given graphically ir 
figure 6. 
The values of S, Wo .., and Wy 4... Computed by means of the experimental} 
data are given in table 4. 
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Figure 6. Velocity of sound in A—H, mixtures. 


The values of S,’ are computed by means of the formula 


: : 3 ; 5 | 
10°. B=0-8742 >> 818 5 1515 x 10 a4 0099 x 10 Gate) 


T iB ie 
published by Keesom and van Lammeren (1932). 
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Table 4 
Concentration (% 
a RS sag ° exp. 10°. S UA ee 
CX.) nN H, (m./sec.) (m./sec.) 

90 100-0 0) — —6-6 —— 
85-6 14-4 189-7 —5:21 189-7 
54-6 45-4 232-0 —2-38 232-0 
36:1 63-9 278°5 —1-36 278-2 
19-5 80:5 358°5 —0-36 385-2 
0 100-0 — +0-29 — 

85 100-0 0 —10°3 — 
85-6 14-4 184-5 — 8:02 184-3 
54-6 45-4 DISD) — 3-62 224-4 
Sol 63°9 270-4 — 1:98 270°2 
19°5 80°5 348-0 — 0:62 348-2 
0 100-0 — + 0:24 = 

80 100-0 0 —- —12-6 — 
85-6 14-4 178°8 — 9-85 178:9 
54-6 45-4 218-6 — 5:17 218-8 
36°1 63-9 262:°3 — 3:36 262:2 
19-5 80°5 837-6 — 1:47 337-8 
0 100-0 oa + 0-17 — 
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Figure 7. Velocity of sound in A-H, mixtures as functions of temperature. 
We computed f by means of equation (8), replacing Sj, S. by S;'and S,’. The 
values of S,’ are calculated by means of the equation 
10? 10° 10° 10° 
10?. By =0-129 + 0-5923 Vr. 75923 et S209 Ts 1-40030 ar Sar: (11) 
published by van Itterbeek and van Paemel (1938). The values of S,' are 
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computed by formula (10), the specific heats for hydrogen being taken from the | 


book of Farkas (1935). | 
The values obtained for f for the different concentrations are given in table 5. jij 


Table 5 
T x 10°. B 
103. 
(OK) (2) average 
90 80 dey 
70 —1-8 
60 Ae 
50 —1-8 ae 
40 =1:8 
30 —1-8 
85 80 —4-8 
70 —4-7 
60 —4.5 
50 43 ee 
40 — 4.3 
30 459 
80 80 —$8-5 
70 = 99 
60 Ley) 
50 ng pe 
40 — 3-4 
30 = 75 


We see, as remarked in § 2, that 6 does not dened upon pressure within the limits} 
of experimental accuracy. 


§6. MEASUREMENTS ON A-D, MIXTURES 


In order to compare the interaction of H, with that of its isotope D,, we havel} 
also made some measurements on A-—D, mixtures. ‘These measurements, | 
indicated in figure 8, are only made for a concentration at which the influence of the} 
interaction is largest. | 


Figure 8. Results for A-D, mixtures. 
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From the experimental data of figure 8, we have computed the values of S 
and W, a (see table 6). W, theor, Was calculated by means of the values of the 
specific heats of D, taken from Farkas. 


Table 6 

| 

eve Concentration (%) 

i a ee ee 

(°K.) iN De! Wo exp: 10° .S Wo theor. 

90 68-4 Bile O 200°5 —3-64 200-7 
85 68-4 31-6 OSA —5-34 195-0 
80 68-4 31-6 189-2 —7:34 189-3 


In table 7 we compared the values of S obtained for A-D, mixtures with 
those of A-H, mixtures. It will be seen that there is no difference within the 
limits of experimental error. 


‘Table 8 
Dizrarey | Oe Oa a 1027 Ba | Orns 
‘Table 7 , (OE) 

cL 102258 103. S A-H, 90 — 8-0 —(\ AG 
©.) A-D, A-H, 85 =eo°0) —0°3 Sas 
80 — 10-0 —0-6 fell 

90 — 3-64 e655 
gs ge oe A-He 90 — 8-0 0-5 aes 
80 —7:34 — FAG) 85 — 9-0 0-38 4.3 
80 — 10:0 0:28 75 
H,—He 20 — 6:56 0-10 = ei li| 


§7. CONCLUSIONS 


In table 8 we have summarized our results obtained for the interaction terms 
and compared them with the interaction betweeu the individual components. 
Thus we can observe that £ is generally of the same order as the B-value of the 
component with the higher boiling point. pe aoe when the temperature 
decreases, 6 increases very Sirota 

Further work is in hand on H,-N,, H,—-He, H,-O, and He—O, mixtures. 
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REVIEWS OF BOOKS 


Optique Physiologique : La Dioptrique de L’iil et sa Correction, by Y. Le GRAND} 
Pp. 356 and 105 diagrams. (Paris: Revue d’Optique, 1946.) | 


The book under review represents Part I of a treatise by M. Le Grand on Physiologica 
Optics. He is here dealing more particularly w‘th the optical system of the eye on 
Gaussian principles, and with the method of testing and correction of refractive error i 
The main contents of the book will be generally familiar to those who have studied visua} 
optics, but the material is presented in a very clear and orderly manner. A necessary} 
amount of algebra is included, but the text is not overburdened with it, neither is tha 
clarity of the diagrams obscured by too many rays of light or too many details of the optica 
systems. 

Since M. Le Grand has for many years been engaged in teaching physiological optics 
it is not surprising to find a list of exercises given at the end of the book. As the solution} 
of the problems, and not merely the answers, are also included, it is possible for thi 
conscientious student both to test and teach himself by this means. | 

A more unusual, but none the less welcome, feature in a book of this type is the inclusioaf} 
of two appendices dealing respectively with elementary statistics and with the calculation} 
of matrices in optical work. Whether the average optician will consider that he has an 
need to venture into these fields is perhaps a matter of doubt, but at least the material iif} 
provided in a readily assimilable form for those who have the inclination to attack it. 

It is interesting to note that in the list of original papers quoted in the bibliographjf 
at the end of the book, half of them appeared in the Tvansactions of the Optical Society— 
a small matter, perhaps, but significant that the Society in its time performed a usefu 
function and might have continued to do so if sufficient support had been forthcomin}h 
to keep the Society alive. Readers in this country will also be interested to see thaf| 
Fincham’s theory of accommodation is described as one of the most satisfactory. 

Quite obviously M. Le Grand regards his subject as one worthy of very serious; tudi 
by optical designers. As he remarks on p. 345, the modern correcting lens is a com} 
plicated optical instrument and the last word has not yet been said about the best desigufi 
of spectacle Jenses. M. Le Grand can at least claim to have given a clear account of thip 
basic principles which must govern any future developments in such lens design. 

W. D. WRIGHIp 


A Survey of General and Applied Rheology, by G. W. Scott Brair. Pp. 194) 
(London: Sir Jsaac Pitman and Sons, Ltd., 1946.) 18s. 6d. | 


Rheology is defined as the study of the deformation and flow of matter. Thus define 
it should obviously cover a large proportion of the field of physics. In practice thi 
subject is concerned almost exclusively with those types of matter which do not obe}+ 
the simple laws of deformation or flow to be found in the elementary textbooks of physics}! 

Dr. Scott Blair’s book is written primarily for the scientist in industry, who has t 
deal with matter as it is rather than with matter as the physicist would like it to be. Th 
first part of the book gives an account of the wide variety of rheological phenomen . 
encountered in industrially important materials, and discusses the various types of apparatul| 
by which these phenomena may be investigated. "Those wishing to pursue the subject inf 
greater detail will find the annotated bibliography particularly helpful. | | 

The second, and more controversial, part of the book deals with the interpretation a 
rheological data. The central question here discussed may be stated in the following way | 
What physical factors are responsible for the formation of a subjective judgment of a giver 
theological property on the basis of sensory impressions? Dr. Scott Blair has mad 4 
extensive investigations into this question, and his discussion of this and related proble 


makes interesting reading. L. R. G. TRELOABH. 
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